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Abstract 

In this thesis we investigate the stability of few body symmetrical dynamical 
systems which include four and five body symmetrical dynamical systems. 
We divide this thesis into three parts. In the first part we determine some 



QJ • special analytical solutions for restricted, coplanar, four body problem both 

^ \ with equal masses (Roy and Steves, 1998) and with two pairs of equal masses. 



The Lagrange solutions L\ and L4 are obtained in the two triangular solutions. 

The equilateral triangle and the isosceles triangle solutions are also obtained. 

We also provide a comprehensive literature review on the four and five body 

j> ■ problems to put our research on these problems in the wider context. 

(N 

' hi the second part we investigate more complicated and general four body 

§: — * - 

O ■ donian Symmetric Four Body Problem (CSFBP), a symmetrically restricted 



four body configuration first introduced by Roy and Steves (1998), by perturb- 
ing the position of one of the bodies and using the general four body equations. 
We show that the CSFBP is stable towards small perturbations and there is 
no significant change in the symmetry before and after the perturbations. 

In the third part we introduce a stationary mass to the centre of mass of 
the CSFBP, to derive analytical stability criterion for this five body system 
and to use it to discover the effect on the stability of the whole system by 
adding a central body. To do so we define a five body system in a similar 



fashion to the CSFBP which we call the Caledonian Symmetric Five body 
Problem (CS5BP). We determine the maximum value of the Szebehely con- 
stant, C = 0.659, for which the CS5BP system is hierarchically stable for all 
mass ratios. The CS5BP system has direct applications in Celestial Mechan- 
ics. The analytical stability criterion tells us for what value of C the system 
will be hierarchically stable but we would like to know what happens before 
this point. To understand this, we determine a numerical stability criterion for 
the CS5BP system which compares well with the analytical stability criterion 
derived earlier. We conclude this thesis with the generalization of the above 
analytical stability criterion to the n body symmetrical systems. This new 
system which we call the Caledonian Symmetric N Body Problem (CSNBP) 
has direct applications in Celestial Mechanics and Galactic dynamics. 

Research presented in this thesis includes the following original investiga- 
tions: determination of some analytical solutions of the four body problem; 
stability analysis of the near symmetric coplanar CSFBP ; derivation of the 
analytical stability criterion valid for all time for a special symmetric config- 
uration of the general five-body problem, the CS5BP, which exhibits many of 
the salient characteristics of the general five body problem; numerical investi- 
gation of the hierarchical stability of the CS5BP and derivation of the stability 
criterion for the CSNBP. 
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Chapter 1 
Introduction 



The motion of systems of three or more bodies under their mutual gravitational 
attraction is a fascinating topic that dates back to the studies of Isaac Newton. 
Because of the complicated nature of the solutions, few-body orbits in the 
most general cases could not be determined before the age of computers and 
the development of appropriate numerical tools. Today the few body problem 
is recognized as a standard tool in astronomy and astrophysics, from solar 
system dynamics to galactic dynamics (Murray and Dermot, 1999). 

During the past century, Celestial Mechanics has principally been devoted 
to the study of the three body problem. Due to the difficulty in handling the 
additional parameters in the four and five body problems very little analytical 
work has been completed for greater than three bodies. 

In our galaxy it is estimated that roughly two-thirds of all stars exist in 
binary systems. Furthermore it is estimated that about one-fifth of these sys- 
tems actually exist in triple systems, while a further one-fifth of these triple 
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systems are believed to exist in quadruple or larger multiple systems. There- 
fore, out of the 10 11 stars in the Galaxy, of the order of 2 x 10 9 of them exist in 
quadruple stellar systems (Steves and Roy, 2001). It is therefore of important 
to study the dynamical behavior of such systems. 

In this thesis we investigate the dynamics of small clusters of stars and 
planetary systems using n-body symmetrical dynamical systems where n > 4. 
We have generalized the investigations of Steves and Roy on symmetrical four 
body problems to symmetrical N-Body problems which gives us insight into 
the stability of symmetrical stellar clusters with planetary systems. 

The main thesis of research is divided into two parts. In the first section 
(chapters 3), we discuss the equilibrium configurations of four body problems 
as particular solutions of four body problems. In the second section (chapter 4 
to 7), the more complicated four and five body restricted symmetrical problems 
are studied. We particularly focus on the Caledonian Symmetric Four Body 
Problem (CSFBP) Steves and Roy (1998, 2000, 2001) and the Caledonian 
Symmetric 5 Body problem (CS5BP). 

The current state of research on the stability of the four and five body 
problems is reviewed in chapter 2. Because of the greater complexity of hav- 
ing a higher number of bodies, the main focus of the literature on the four 
and five body problems has been on the analytical study of their Equilibrium 
Configurations . An equilibrium configuration of four-bodies is a geometric 
configuration of four bodies in which the gravitational forces are balanced by 
the centrifugal force so that the configuration is maintained for all time. We 
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also review the few existing papers which address the analytical stability of 
four body problems in general: In particular the papers of Loks and Sergy- 
sels (1985) and Sergysels and Loks (1987) are considered. Publications on 
the equilibrium solutions and the stability of the Caledonian Symmetric Four 
Body Problem ( Steves and Roy (1998, 2000, 2001), Roy and Steves (1998) 
and Szell (2003)) form the main source of knowledge relevant to the research 
of this thesis. These papers are reviewed separately in detail at different places 
of the thesis. 

In chapter 3, we review the equilibrium solutions of symmetric four body 
problems given by Roy and Steves (1998) and derive some further solutions. 
In section 3.1 we review three types of analytical solutions for the equal mass 
four body problem which include the Square, the Equilateral triangle and the 
Collinear equilibrium configurations (Roy and Steves, 1998). In section 3.2 we 
give the solutions for two pairs of equal masses where the ratio between the 
two pairs is reduced from 1 to in order to obtain equilibrium configurations 
which involve the five Lagrange points of the Copenhagen problem. We discuss 
four kinds of equilibrium configurations both symmetric and non-symmetric 
which include the Trapezoidal, the Diamond, the Triangular and the Collinear 
equilibrium configurations. The derivation of the Trapezoidal, the Diamond 
and the Collinear equilibrium configurations are a review from Roy and Steves 
(1998), while the derivation of the two Triangular equilibrium configurations 
is original work. 

It would be useful to broaden this understanding to include the stability of 
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general solutions to restricted four body problems. In chapter 4 we look at the 
stability of a restricted case of the four body problem called the Caledonian 
Symmetric Four Body Problem (CSFBP). This symmetrically restricted four 
body problem was first developed by Steves and Roy (1998). 

Steves and Roy (2001) later derived an analytical stability criterion valid 
for all time for the CSFBP. They show that the hierarchical stability of the 
CSFBP solely depends on a parameter they call the Szebehely constant, C , 
which is a function of the total energy and angular momentum of the system. 
This stability criterion has been verified numerically by Szell, Steves and Erdi 
(2003a, 2003b). Szell, Erdi, Sandor and Steves (2004) analyze the connection 
between the chaotic behavior of the phase space and the global stability given 
by the Szebehely constant using the Relative Lyapunov indicator (RLI) and 
the smallest alignment indices (SALI) methods. They found that as the Sze- 
behely constant is increased, making the system hierarchically stable from a 
global point of view, the corresponding phase space becomes increasingly more 
regular. 

In chapter 4 we review their research on the CSFBP and then investigate 
the CSFBP phase space using nearly symmetric, perturbed, initial conditions 
and the general four body equations in order to study if the CSFBP system 
will remain symmetric under slightly perturbed initial conditions. We use 
long time integrations of a million time-steps to study the behavior of the 
slightly perturbed CSFBP phase space. An integrator is specifically developed 
for this purpose using the Microsoft Visual C++ Software. The results of 
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integrations are processed using Matlab 6.5. During the integrations, we record 
the following observations: 

1. For all orbits which have a close encounter, we stop the integration and 
record the type of the collision or the close encounter. 

2. For all orbits which fail the symmetry breaking criterion, we stop the 
integration and record the symmetry as broken . 

3. For all orbits where there are no collisions or breaking of the symmetry, 
the integration continues to the end of 1 million time-steps. 

In Chapter 5 we introduce a stationary mass to the centre of mass of 
the CSFBP and derive an analytical stability criterion for this new five body 
symmetrical system similar to the one of the CSFBP. This stability criterion 
enables us to analyze the effect of the additional central body on the stability 
of the whole system. We call this new problem the Caledonian Symmetric Five 
body Problem (CS5BP). The CS5BP has direct applications in dynamical sys- 
tems where a very large mass exists at the centre of mass with four smaller 
masses moving in dynamical symmetry about it. The four small masses are 
affected by the central mass but are small enough that they do not dynami- 
cally affect the central body. Hence the central body remains stationary and 
dynamical symmetry is maintained. This could occur, for example, in exo- 
planetary systems of a star with four planets or a planetary system with four 
satellites. 
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For completion we analyze the full range of mass ratios of central body to 
the other four bodies: from the large central body system described above to 
a five body system of equal masses to a small central body with four large sur- 
rounding bodies. In the case of five equal masses or a small central body with 
four surrounding large masses, the central body is unlikely to remain station- 
ary as is required by the CS5BP. The model may, however, still be applicable 
to real systems in which the outer bodies are well spaced and stationed far 
away from the central body so that they have minimal effect on the central 
body. 

In chapter 6 we continue our analysis of the CS5BP. We investigate nu- 
merically the hierarchical stability of the CS5BP. The main objective of this 
exercise is to validate the hierarchical stability criterion derived in chapter 6 
and to study the relationship between the number of hierarchy changes and 
the value of C . 

Szell, Steves and Roy (2002) provide a numerical investigation of the hier- 
archical stability of the CSFBP which covered half of the phase space. We pro- 
vide a brief review of these results in section 6.1 before deriving the equations 
of motion for the CS5BP in section 6.2. The comprehensive numerical explo- 
ration was completed using an integrator specially developed for the CS5BP in 
Fortran and C++. We used Matlab 6.5 software to process the results. 3000 
orbits were integrated to 1 million time steps of integrations for each of 63 
values of Co. The total number of orbits integrated were 3000 x 63 which took 
about 70 days of CPU time. During the integrations the following information 
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was recorded for each value of C to be able to comment on the hierarchical 
stability of the system under discussion: 

1. The total number of hierarchy changes and 

2. The types of hierarchy changes 

In section 6.5, we give a complete analysis of the hierarchical stability of the 
CS5BP for a whole range of mass ratios. This includes the completion of the 
analysis of Szell et al. (2002) and Szell, Steves and Erdi (2004a) as the CSFBP 
is a special case of the CS5BP. 

In chapter 7, we generalize the Caledonian Symmetric 5 Body Problem 
(CS5BP) to the Caledonian Symmetric N Body Problem (CSNBP) to derive 
an analytical stability criterion for the Symmetric N Body problem where 
N > 4. 

Finally in chapter 8 we give a brief summary of the results obtained in this 
thesis and in chapter 9 we point out some possible areas for future exploration. 

In this thesis the original research involves: the determination of the trian- 
gular equilibrium solutions of the four body problem (sections 3.3.3 and 3.3.4); 
the stability analysis of the near symmetric coplanar CSFBP (sections 4.4 and 
4.5); the derivation of the analytical stability criterion valid for all time for the 
symmetric configuration of the general five-body problem, the CS5BP, (Chap- 
ter 5); the numerical investigation of the hierarchical stability of the CS5BP 
(sections 6.2 to 6.7) and the derivation of the stability criterion for the CSNBP 
(Chapter 7). 
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Now in chapter 2 we present a brief review of the current state of the 
literature on analytical studies of the four and five body problems. 
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Chapter 2 

The Four and Five Body 
Problem : A literature review 



In our galaxy it is estimated that roughly two-thirds of all stars exist in binary 
systems. Furthermore it is estimated that about one-fifth of these systems 
actually exist in triple systems, while a further one-fifth of these triple systems 
are believed to exist in quadruple or larger multiple systems. Therefore, out of 
the 10 11 stars in the Galaxy, of the order of 2 x 10 9 of them exist in quadruple 
stellar systems (Steves and Roy, 2001). It is therefore of interest to study the 
dynamical behavior of such systems. 

The classical equation of motion for the n-body problem assumes the form 



where the units are chosen so that the gravitational constant is equal to one, 





(2.1) 
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i-j is a vector in three space, 



u= J2 TT=H (2 ' 2) 



is the self-potential, r*; is the location vector of the ith body and mi is the mass 
of the it\i body. 

Because of the greater complexity of higher number of bodies, the main 
focus of the literature for four and five body problems has been on the analyt- 
ical study of the Equilibrium Configurations. We therefore largely concentrate 
our review on this research. We also review the few existing papers which 
address the general four body problem: In particular the papers of Loks and 
Sergysels (1985) and Sergysels and Loks (1987). Publications on the equilib- 
rium solutions of the four body problems and the stability of the Caledonian 
Symmetric Four Body Problem ( Steves and Roy (1998, 2000, 2001), Roy and 
Steves (1998) and Szell (2003)) form the main source of knowledge relevant to 
the research of this thesis. 

In this chapter we review the analytical research on both four and five 
body problems. In section 2.1 a literature review is given for the analytical 
solutions and stability for four body problems and in section 2.2 for the five 
body problems. 
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2.1 The Four Body Problem : A literature re- 
view 

We divide this section further into four subsections. The first subsection deals 
with particular solutions, the Central Configurations, of the Four Body Prob- 
lem. The remaining three subsections deal with the stability analysis of various 
degrees of restricted four body problem: the restricted four body problem, the 
symmetric four body problem and the general four body problem. 

2.1.1 Central Configurations of the Four Body Problem 

To understand the dynamics presented by a total collision of the masses or the 
equilibrium state for a rotating system, we are led to the concept of a central 
configuration. For a system to be the central configuration, the acceleration of 
the ith mass must be proportional to its position (relative to the centre of mass 
of the system), thus f\ = Arj for alH = 1, 2, n. A Central Configuration can 
also be expressed as a critical point for the function U 2 I, where I is the moment 
of inertia. Though one has two different formulations of a central configuration, 
a number of basic questions concerning them remain unanswered, such as "for 
a fixed number of particles with equal mass, do there exist a finite number of 
Central Configurations? 

In this section a review of papers concerning central configurations and 
hence equilibrium solutions which is a special case of central configurations is 
provided. An equilibrium configuration of four-bodies is a geometric configura- 
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tion of four bodies in which the gravitational forces are balanced in such a way 
that the four bodies rotate their centre of mass and thus the geometric configu- 
ration is maintained for all time. Equilibrium solutions of the symmetric four 
body problem are studied in chapter 3 and chapter 4 in greater detail. These 
chapters include a review of Steves and Roy (2001) along with original work. 

The straight line solutions of the n-body problem were first published by 
Moulton (1910). He arranged n masses on a straight line so that they always 
remained collinear and then solved the problem of the values of the masses 
at n arbitrary collinear points so that they remained collinear under proper 
initial projections. This paper is not particularly concerned about the specific 
case of the four body problem but it can be deduced that there are 12 collinear 
solutions in the equal mass four body problem. These solutions are sometimes 
referred to as Moulton Solutions. 

The two papers on the Classification of Relative Equilibria by Palmore(1975, 
1982) presented several theorems on the classification of equilibrium points in 
the planar n-body problem. He then applied his results to the cases of three 
and four body equal mass systems. These highly technical papers show that 
the theorems lead to the existence of 120 possible configurations of relative 
equilibrium for four equal masses, of which 24 exist in the isosceles configura- 
tion. This number takes into account all of the possible permutations of the 
masses in the same configurations including the 12 Moulton collinear solutions. 

Zhiming and Yisui (1988) in their article under the title, 'The Central 
Configurations Of the General 4-Body Problem', obtain the equations of Cen- 
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tral Configurations using a method different from that introduced by Palmore 
(1975). Zhiming and Yisui (1988) also investigate the finiteness of Central 
Configurations for the general four-body problem. They show that for the 
collinear four-body problem there are twelve central Configurations for each 
set of masses. This agrees with the existence of 12 Moulton (1910) solutions. 
They also prove the following very important theorems: 1. There is not any 
point-line type Central Configuration, where three of the four mass-points lie 
on a line while the fourth does not, in the general four-body problem. 2. For 
any four-body Central Configurations if there is a mass point which keeps the 
same distance from two of the other three mass points, the fourth one also keeps 
the same distance from the two. Such a Central Configuration is symmetric. 

Using algebraic and geometric methods, Arenstrof (1982) investigates the 
number of equivalence classes of central configurations in the planar four-body 
problem with three arbitrary sized masses and a fourth small mass m 4 . For 
example the following system of algebraic equations is algebraically reduced 
to three equations in three unknowns and then the problem is reduced to the 
case 7714 = 0. 

N 

fi = h = - = fN = 0, fk = 22-f(qk-Qj), r jk = \qj -q k \>0. 

j=i r « 

His results show that each three-body collinear Central Configuration generates 
exactly two non-collinear Central Configurations (besides four collinear ones) 
of four bodies with small m 4 > ; and that the three body equilateral triangle 
Central Configuration generates exactly eight, nine or ten planar four-body 
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Central Configurations with m 4 = 0. 

Glass (1997) studies the Central Configurations of the classical N-body 
problem and the asymptotic properties of a system of repelling particles. An 
asymmetric configuration obtained in the eight-particle system is described 
and a bifurcation in the four-particle system is investigated. An asymmetric 
configuration with eight particles, which is the smallest number of particles 
with no axis of symmetry, is found. Bifurcation occurs in the general system of 
particles as the parameter (3 is varied, where U'(r) = U'(l)/r 1+l3 . A particular 
example with four particles is discussed and the following four equilibrium 
configurations for this system with (5 — 1 is given. 

1. A configuration with four particles lying along a line through the origin. 

2. A configuration with particles at the vertices of a square. 

3. A configuration with particles at the vertices of an equilateral triangle 
and one particle at the origin. 

4. A configuration with particles at the vertices of an isosceles (but not equi- 
lateral) triangle and one particle within the triangle and on its symmetry 
axis. 

Simo (1978) presents the classification and solutions of the central configu- 
rations of the four body problem using topological proofs. Simo first considered 
the restricted problem of three finite masses plus one infinitesimal mass. He 
showed diagrammatically the location of equipotential curves and the nine 
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equilibrium points in the case of three equal masses and in the case of masses 
mi = 0.2, m 2 = 0.3 and m 3 = 0.5. He then moves on to the problem of four 
arbitrary masses. He found (a) a division of the mass space in regions for 
which there is a different number of relative equilibria, (b) masses for which 
degeneration occurs, and (c) for given masses, the determination of the rel- 
ative equilibria and the relation with the corresponding relative equilibria of 
the equal mass case. 

Several papers (Majorana (1981), Kozak and Oniszk (1998), Gomatan, 
Steves and Roy (1999)) derive the equilibrium solutions and analyze their 
stability for different types of four body problems. Majorana (1981) studies the 
linear stability of the equilibrium points in the restricted four body problem, 
where three bodies, of masses //, /x and 1— 2/i rotate in an equilateral triangular 
configuration (the Lagrange solution), whilst the fourth body of negligible mass 
moves in the same plane. The equations of motion of the particle under the 
influence of the other three bodies were derived which led to the determination 
of eight equilibrium points in the problem. The author found numerically that 
the linear stability of some of these equilibrium points depended on the masses 
of the three bodies. Five of the points were found to be unstable for all values 
of /i, whilst two of the other points were found to be stable for small values of 

Kozak and Oniszk (1998) studied the motion of a negligible mass in the 
gravitational field generated by a collinear configuration of three bodies (of 
masses m ^ mi = m 2 = m). To study the linear stability of the equilibrium 
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points, the authors write the equations of motion in their Cauchy form follow- 
ing the definition of the equilibrium points. Because it is impossible to find a 
compact analytic form for the solutions, they find the approximate solutions 
by: 

1. expanding terms as a power series in m/m and m /m. 

2. using Newton's Iterative method for given values of m/m . 

They proved that the straight line (collinear) configurations are linearly 
unstable for any value m as in the Lagrange case of the three-body prob- 
lem. Intervals are given for the stable and unstable regions of the triangular 
configurations. 

This review on central configurations gives the context for the research of 
chapters 3, where the work by Steves and Roy (2001) on equilibrium solu- 
tions of symmetric four body problem is reviewed and extended to include 
derivations of the triangular equilibrium solutions. Roy and Steves (1998) dis- 
cuss some special analytical solutions of four body problems. They show that 
these solutions reduce to the Lagrange solutions of the Copenhagen problem 
when two of the masses are equally reduced. The content of this paper will be 
discussed in chapter 3 in detail. 

2.1.2 Symmetric Four Body Problem 

It is well known that it is not possible to find analytical solutions for the gen- 
eral four body problem. Therefore restriction methods utilizing assumptions 
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of neglecting the masses of some of the bodies or assumptions which require 
specific conditions of symmetry have been used to reduce the dimensions of 
the phase space to manageable levels while still producing results which are 
meaningful to real physical systems. For example a four body model requiring 
symmetrical restrictions was used by Mikola, Saarinen and Valtonen (1984) as 
a means of understanding multi-star formation in which symmetries produced 
in the initial formation of the star system were maintained under the evolution 
of the system. In this section papers which used specific conditions of symme- 
try to simplify the four body problem analysis will be reviewed. These include 
a series of papers on the Caledonian Symmetric Four Body Problem (CSFBP) 
which are very relevant to the research presented in this thesis, in chapters 5 
to 8. 

Roy and Steves (2000) first introduced a symmetric problem which they 
called the Caledonian Symmetric Double Binary Problem (CSDBP). It is a 
special case of the Caledonian four-body problem introduced earlier by the 
same authors (Roy and Steves, 1998). To form the CSDBP, they utilized 
all possible symmetries. The CSDBP is a three dimensional problem, with 
two pairs of masses, each pair binary having unequal masses but the same 
two masses as the other pair. They have shown that the simplicity of the 
model enables zero-velocity surfaces to be found from the energy integral and 
expressed in a three dimensional space in terms of three distances ri,r 2 and 
ri2, where r\ and r 2 are the distances of the two bodies which form the pair 
from the centre of mass of the four body system. r i2 is the separation between 
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the two bodies. This problem was later renamed as the Caledonian Symmetric 
Four Body Problem (CSFBP). Further details of the CSFBP are given in 
chapter 5. 

Steves and Roy (2001) derived an analytical stability criterion for the 
CSFBP valid for all time. They show that the hierarchical stability of the 
CSFBP depends solely on a parameter they call the Szebehely Constant, Co, 
which is a function of the total energy and angular momentum of the system. A 
four body system is said to be hierarchically stable if it maintains its initial hi- 
erarchy state for all time. This analytical stability criterion was numerically 
verified by Szell, B. Steves and Roy (2002) and Szell, Steves and erdi (2004a) 
for the coplanar CSFBP. They have performed a wide range of numerical inte- 
grations for a variety of values of the Szebehely constant Co. They also show 
that for C greater than a critical value, the system is hierarchically stable for 
all time and will undergo no change in its hierarchical arrangement. They also 
show that the number of hierarchy changes decreases with the increasing value 
of C . 

Szell, Erdi, Sandor and Steves (2004) investigate the chaotic and stable 
behavior of the CSFBP using the relative Lyapunov indicator (RLI) and the 
smallest alignment indices (SALI) which are fast chaos detection methods. 
They analyze the phase space of the CSFBP in detail for different mass ratios 
and for each mass ratio a range of values C . The color coded points on their 
graphs give the nature of the orbits (chaotic, regular, ending in collisions) 
defined by different initial conditions of the CSFBP. They show that the regular 
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and chaotic behavior of the phase space is closely connected with the Szebehely 
constant. The larger the Szebehely constant, the more regular the phase space 
becomes. 

Szell, Steves and Erdi (2004b) give a numerical escape criterion for the 
CSFBP. By approximating the symmetrical four body problem on the verge 
of break up as a two body problem, they derive four energy like parame- 
ters Ei, E 2 , E 3 and E A which are related to each of the four types of possible 
escapes. Numerous orbits were investigated with the help of the numerical 
escape criterion. They found that for stellar systems of nearly equal mass the 
most likely escape configuration is the double binary escape, i.e. the system 
falls apart into two binary-star configurations. For small values of the mass 
ratio, /x, the system can be a model for two stars two planets systems. The 
integrations show that the most likely escape in this case is the escape of the 
two planets. 

Sweatman (2002) discusses a symmetrical four body system where the bod- 
ies are distributed symmetrically about the centre of mass. He also restrict 
the bodies to move upon a line undergoing elastic collisions. It is shown by 
numerical simulations that there is a similarity between this system and the 
one dimensional newtonian three-body problem. The non-linear stability of 
the Schubart-like orbit is studied using long numerical integration. It is shown 
that this orbit is unstable in three dimensions. 

Jiang-Hui, Xin-hao and Lin (2000) study a particular case of the 4-body 
problem in the plane with a rhomboidal configuration. The rhomboidal four 
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body configuration means the four particles with masses mi, m 2 ,m 3 , 777,4 re- 
spectively are located in the plane at the vertices of a rhombus. The particles 
are given symmetric initial conditions in position and velocities with respect 
to the axes in the plane and they thus always keep a symmetric rhomboidal 
configuration under the law of newtonian attraction. Jiang-Hui et al. (2000) 
show that the larger the mass ratio (a = m 3 /mi), the more extended is the 
region of stable motion in the phase diagram. 

Vidal (1999) studies a symmetric equal mass four body problem, The Tetra- 
hedral 4-Body Problem with Rotation, which is a particular case of the general 
four body problem in space. It is obtained when two bodies form a binary pair 
mi and m 2 that are always symmetric to one another with respect to the ver- 
tical axis, and the other two bodies form a pair m 3 and 777,4 that are symmetric 
to one another with respect to the same axis, with all the masses equal to each 
other. The total angular momentum is then forced to be zero by making the 
angular momentum of the binary (mi — m 2 ) equal and opposite to the angular 
momentum of the other binary (m 3 — m.4). He proves that all singularities of 
this model are due to collision. He also proves that the singularities due to 
simultaneous double collisions are regularizable. The set of equilibrium points 
on the total collision manifold is studied, as well as the possible connections 
among them. It is shown that the set of initial conditions on a given energy 
surface going to quadruple collision is a union of twenty submanifolds. 
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2.1.3 Restricted Four Body Problem 

Four body problems found in nature can often be simplified using the reduction 
of variables given by having some masses which are negligible compared to 
others. Thus they have no gravitational effect on the finite masses, yet all 
finite masses have an effect on them. A restricted four-body system of this type 
usually contains three bodies with non-zero masses (primaries) and a fourth 
massless particle. Here the time evolution of the massless particle's orbit can 
be studied. In the case of the restricted three-body problem, (problem with two 
primaries and a third massless particle) the orbits of the two primaries were 
well known but in the case of the restricted four-body problem the motions of 
the three primaries cannot be given analytically. Thus some assumptions on 
their motions must be made. 

The "very restricted four-body problem" was introduced by Huang (I960) 
being a circular restricted four-body system. In the circular restricted four- 
body system, it is assumed that the primaries have circular orbits i. e. orbit on 
a circle. He defined the masses of the bodies in the system be m 1; m 2 , m 3 and 
m 4 such that m-y ^> m 2 + m 3 and m 4 = the massless particle. He assumed 
that r 2 3 <C r 12 ,r 13 for all time, where denotes the distance between the ith 
and jth bodies. If 7712 and m 3 revolve on a circular orbit and their barycenter 
revolves on a circular orbit about the mi body and all bodies are in the same 
plane then the system is defined to be a "very restricted four-body system". 

Similarly to the restricted three-body problem, zero velocity curves can 
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be defined in this case (Huang, 1960), i.e. surfaces can be determined within 
which motion can take place. Outside these surfaces, motion is impossible; 
however in the "very restricted four body problem", the surfaces change their 
shape under the evolution of the system and so osculating zero velocity curves 
must be introduced. Applying this method to the Sun-Earth-Moon-satellite 
very restricted four-body system (m x = m Su n, rn 2 = m Earth , m 3 = m M oon, 
7714 = rrisateiiite = 0) it can be found that any artificial satellite orbiting the 
Moon must have a closed orbit in order to be stable (Huang, 1960). 

Matas (1968) studies a special restricted four body problem. He considers 
four material points, three of them have finite non-zero mass and the fourth 
represents an infinitesimal body. He also requires a constant configuration for 
the three finite masses i.e. the motion of the first three points is given by the 
particular Lagrange's solution of the problem. He obtains an integral which 
is a generalization of the Jacobi's integral. As a consequence of this integral, 
the equation of the surface of zero relative velocity (a generalization of Hill's 
surfaces) is derived. 

Mohn and Kevorkian (1967) consider the motion of a particle in the gravi- 
tational field of three mass points (Sun-Earth-Moon) that obey the restricted 
three body equations. In particular, they concentrate on the two cases where 
the three primaries move according to Hills theory for periodic lunar orbits. 
In these cases the leading terms in the solar perturbation are included in a 
systematic and dynamically consistent manner so that the results can then 
be derived in a form where, for very close configurations of earth, moon, and 
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particle, the solar influence introduces a perturbation upon a motion which 
would otherwise obey the restricted three-body problem for earth, moon, and 
particle. 

Hadjidemetriou (1980) studies the motion of a small planet moving in the 
gravitational field of the Sun- Jupiter-Saturn system. He first gives the initial 
conditions of a periodic orbit in the general three body problem of the Sun- 
Jupiter-Saturn system, using a suitable rotating frame of reference. He then 
derives the equations of motion for a massless body moving in the Sun- Jupiter- 
Saturn system in the rotating frame. The approximate initial conditions for a 
periodic orbit of the massless body are then determined. The stabilities of the 
periodic orbits are examined using numerical techniques. 

Michalodimitrakis (1981) introduces the restricted four body problem by 
generalizing the restricted three body problem. He studies a special case of the 
restricted four body problem, the circular restricted four body problem, which 
can be considered as a generalization of the circular restricted three body 
problem. Using the method of Szebehely (1967) he finds some equilibrium 
points including four on the x-axis (collinear points) and two on the y-axis 
where all the equilibrium points are symmetrically located with respect to the 
origin. He gives several monoparametric families of periodic orbits. 

Giacaglia (1967) gives the regularization of the restricted four body prob- 
lem, where the same transformations as in the restricted three body problem 
(Szebehely and Giacaglia, 1964) were used to remove the singularities. The 
author then extends the transformation of the equations of motion to the prob- 
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lem of a particle moving under the influence of n equal mass bodies which move 
unperturbed on the vertices of an n-polygon. 

2.1.4 General Four Body Problem 

The general four body problem is the most general case without any restric- 
tions. The Newtonian equations of motion may be written 



They are second order ordinary differential equations for the unknown fi(t) 
functions. It is a 24th degree system since there are four bodies with three 
location coordinates and three velocity coordinates. With the help of the 
energy integral, momenta integrals and the elimination of the node it can be 
reduced to a 14th degree system. 

The equations can not be solved analytically; however some analytical 
statement can be made by examining the energy and momenta manifolds of 
the system. Using the concept of zero velocity curves, some hyper-surfaces 
of zero velocity can be constructed which determine the regions of possible 
motion (Loks and Sergysels, 1985, Sergysels and Loks, 1987). Since these 
hypersurfaces are many dimensional, the visualization of the results are very 
difficult. In the general three body problem, Zare (1976), Zare (1977) and 
Marchal and Saari (1975) showed that the energy H and angular momentum 
c integrals can combine to form a stability criterion c 2 H which, for a given 
critical value c 2 H crit , produces a phase space for the three body system which 
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contain topologically separate subregions. Thus the hierarchically arrange- 
ment of the system which exists in one subregion of the real motion cannot 
physically evolve into the hierarchical arrangement of the system existing in 
another subregion of the real motion. 

In the two papers by Loks and Sergysels (1985) and Sergysels and Loks 
(1987), the c 2 H stability criterion in the three body problem was extended 
to the general four body problem by finding the zero velocity hypersurfaces 
which define the regions in which motion can take place. In their first paper, 
the authors began by describing the kinetic energy and the angular momentum 
of four planar bodies, using generalized Jacobi variables. By transforming the 
equations they then derived expressions for the kinetic and potential energies of 
the system. From these equations an expression for zero velocity hypersurfaces 
was found. In their second paper, the authors applied Sundman's inequality 
for n bodies to the four body problem in the generalized Jacobi coordinates 
in order to derive zero velocity surfaces in three dimensional space. The zero 
velocity surfaces then define the limits of three dimensional regions in which 
the motion is constrained. The cross section of these regions were then plotted 
for an equal mass four body problem. These papers are highly relevant to the 
work done by Roy, Steves and the author. 

Martinez and Simo (1999) consider simultaneous binary collisions in the 
general planar four body problem. They prove that these collision are reg- 
ularizable in the sense of continuity with respect to initial conditions using 
a blow-up of the singularity. They discuss subproblems of the general case 
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which include the particular example of the trapezoidal solution of the equal 
mass four body problem and double isosceles four body problem. Numerical 
evidence is given that the differentiability of the regularization should be, in 
general, less than 8/3. 

Asteroids can be found around the 7 6 secular resonance with Saturn. The 
longitudes of perihelion of these asteroids rotate with the same angular velocity 
value as Saturn's longitude of perihelion. Williams (1979) listed the names of 
the supposed 76 asteroids. Performing more than a 1 Myr numerical integra- 
tion on the elliptical restricted four-body model i.e. the restricted four body 
problem where the three primaries' orbits are ellipses, so that the asteroid is 
allowed to move in an elliptical orbit. Froeschle and Scholl (1987) found that 
there were only two asteroids from Williams's list located deeply in the 76 
resonant region. In order to check the validity of their results, they repeated 
the integration including the effects of all planets (beside Pluto) in the Solar 
System. They found some differences which suggests that in some cases nei- 
ther the restricted three-body model nor the restricted four-body model can 
describe the complete behavior of the asteroid belt. Thus the investigations of 
the five or more body systems is very important to consider. 
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2.2 The Five Body Problem: A literature Re- 
view 

In the past Celestial Mechanics has been mainly devoted to study the three 
and four body problems. Very little has been published with regards to the 
five body problem. Because of the greater complexity of higher numbers of 
bodies, the main focus of the literature has been on the analytical study of cen- 
tral configurations and relative equilibria of the five body problem. A relative 
equilibrium is a special configuration of masses of the n-body problem which ro- 
tates about its centre of mass if given the correct initial momenta. In rotating 
coordinates these special solutions become fixed points, hence the name relative 
equilibria. After an extensive search for literature on the five body problem 
the author is led to believe that our work in chapter 5 and 6 is the first to 
address the issue of the analytical stability of five body problem. Therefore all 
the literature review is on equilibrium solutions and relative equilibria configu- 
rations. The following is a review of the most relevant, perhaps all the papers 
on the five body problem. 

The easiest and most accessible relative equilibria are those configurations 
which contain large amounts of symmetry. Roberts (1998) discusses relative 
equilibria for a special case of the 5-body problem. He considers a configuration 
which consists of four bodies at the vertices of a rhombus, with opposite vertices 
having the same mass, and a central body, see figure (12. ip . He shows that in 
this five body problem for the masses (m 1; m 2 , m 3 , m 4 , m 5 ) = (1, 1, 1, 1, —1/4), 
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Figure 2.1: Set up for the 1+rhombus relative equilibria 

there exist a one parameter family of degenerate relative equilibria where the 
four equal masses are positioned at the vertices of a rhombus with the re- 
maining body located at the centre. As the parameter k varies, one pair of 
opposite vertices move away from each other while the other pair move closer, 
maintaining a fixed length between consecutive vertices. He also shows that 
the number of relative equilibria equivalence classes is not finite. 

Mioc and Blaga (1999) discusses the same problem introduced above in the 
post Newtonian field of ManevQ They prove the existence of monoparametric 
families of relative equilibria for the masses (mo, l,m, l,m), where mo is the 
central mass, and prove that the Manev square five body problem, where k — 1 



1 In the 1920's, Manev proposed a gravitational model based on the A/r + B/r 2 potential, 
with A = [i, a gravitational parameter of the field generating body, B = 3/x 2 /(2c 2 ), with c = 
speed of light. Manev considered that this model could be used as a substitute to general 
Relativity (Mioc and Stavinschi, 2000) 
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with masses (m , l,m, l,m) admits relative equilibria regardless of the value 
of the mass of the central body. A continuum of such equilibria (as in the 
Newtonian field) does not exist in the Manev rhomboidal five-body problem. 

Albouy and Llibre (2002) deals with the central configurations of the 1+4 
body problem i.e. they study the central configurations without collisions that 
are the limit of the central configurations of the five body problem in space, 
when the mass of one of the bodies goes to infinity. He considers four equal 
masses on a sphere whose centre is the 'big' mass. They find four symmet- 
ric central configurations and prove that they all have at least one plane of 
symmetry. He also conjectures that there are exactly five classes of central con- 
figurations, as one can always exchange some bodies and apply some isometry 
or change of scale. 

Ollongren (1988) studies a restricted five body problem having three bodies 
of equal mass m placed on the vertices of the equilateral triangle; they revolve 
in the plane of the triangle around their gravitational centre in circular orbits 
under the influence of their mutual gravitational attraction; at the centre a 
mass of (3m is present where f3 > 0. A fifth body of negligible mass compared 
to m moves in the plane under the gravitational attraction of the other bodies. 
All bodies are considered to be point masses. He discuss the existence and 
location of the Lagrangian equilibrium points. 

Ollongren shows that there are 9 Lagrangian equilibrium points. There 
is a critical value of f3 — 43.181 for which the outer lagrangian points on 
the negative x-axis become stable. For (3 less than this value, all Lagrangian 
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points in the plane are unstable and for (3 larger than this value, the 3 outer 
Lagrangian points are stable and the fifth body will carry out a libration motion 
around them. They also conclude that the central mass has a stabilizing effect 
on the motion of the fifth small body provided the mass of the central body is 
large enough. 

Kalvouridis (1999) considers an n-body ring problem with n — 1 primaries 
of equal mass (m) arranged in equal arcs on an ideal ring and a central body of 
a different mass ((3m) located at the centre of mass of the system. He adds a 
further negligible mass to the system and derives the zero velocity curves and 
zones of stationary solutions for the negligible mass. He does this for different 
values of n which include the five body case and the Copenhagen case. He 
shows that the stationary solutions are arranged on the crossing points of 
concentric circles with radial lines forming equal angles between them. Their 
number depends on the specific values of the parameter v — n — 1 and (3 

2.3 Summary 

Because of the greater complexity of higher number of bodies, the main focus 
of the literature has been on the analytical study of the Equilibrium Configu- 
rations of the four and five body problems. Specially for the five body problem 
nothing is published yet on the analytical stability of the five body problem 
and the authors work in chapter 6 and 7 may be the first of its kind. 

A review of the papers which address the four body problems relevant 
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to our work is given in section 2.1. In particular the papers of Loks and 
Sergysels (1985) and Sergysels and Loks (1987) are amongst the most relevant. 
Publications on the equilibrium solutions and the stability of the Caledonian 
Symmetric Four Body Problem ( Steves and Roy (1998, 2000, 2001), Roy and 
Steves (1998) and Szell (2003)) form the main source of knowledge relevant 
to the research of this thesis. In section 2.2 a comprehensive literature review 
is given for the five body problem. The work of Roberts (1998) is of more 
interest to us as he discuses a symmetric five body problem which is similar to 
our problem. In future it will be interesting to search for the relative equilibria 
of the CS5BP using the method of Roberts (1998). It is also possible to find 
equilibrium solutions for the five body problem using the same method as Roy 
and Steves (1998). 
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Chapter 3 

Equilibrium Configurations of 
the Four Body Problem 

In this chapter we discuss Equilibrium Configurations of four-body problems. 
An equilibrium configuration of four-bodies is a geometric configuration of four 
bodies in which the gravitational forces are balanced in such a way that the 
four bodies rotate together about their centre of mass and thus the geometric 
configuration is maintained for all time. 

In this chapter we present a detailed summary of the method of derivation 
of the equilibrium configurations given by Roy and Steves (1998). We give all 
the analytical solutions given by them including the equal mass cases and the 
cases with two pairs of equal masses. We also give two analytical solutions 
of our own called the triangular solutions of the four body problem. Section 
3.1 contains a review of the Lagrangian solutions of the restricted three body 
problem, in particular the solutions to the Copenhagen problem. The Copen- 
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hagen solutions are useful for the analysis of equilibrium solutions in the four 
body problem as they are the solutions which occur when two of the masses 
are reduced to zero. In section 3.2 we review three types of analytical so- 
lutions of the equal mass four body problems which include the Square, the 
Equilateral Triangle and the Collinear equilibrium configurations (Roy and 
Steves, 1998). In section 3.3 we give solutions for two pairs of equal masses, 
where the ratio between the two pairs is reduced from 1 to in order to obtain 
the five Lagrange points of the Copenhagen problem. We discuss four kinds of 
equilibrium configurations both symmetric and non-symmetric which include 
the Trapezoidal, the Diamond, the Triangular and the Collinear equilibrium 
configurations. The Trapezoidal, the Diamond and the Collinear equilibrium 
configurations are a review of Roy and Steves (1998), while the Triangular 
equilibrium configurations are original work derived as an extension of Roy 
and Steves (1998). In section 3.4 we summarize our results. 



3.1 The Lagrangian Solutions of the Restricted 
three body problem 

It is well known that in the general three-body problem there exists five spe- 
cial equilibrium configurations, where either the three masses mi, m 2 , m 3 are 
collinear or they occupy the vertices of an equilateral triangle. 

The Restricted three body problem has five corresponding special equilib- 
rium solutions where the particle of infinitesimal mass resides at one of five 
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points, usually denoted by the letters L\ y L 2 , L 3 , L 4 and L 5 . The first three 
points are collinear with the two primary masses while the last two points form 
equilateral triangles with the two primary masses (Figure I3TTT) . 

In the Copenhagen problem (Figure [3TT|) . with m x = m 2 = m, if a frame of 
coordinates rotating with constant angular velocity n is centred on the centre 
of mass C of the two finite mass bodies and Pi and P 2 are placed at (-0.5,0) 
and (0.5,0) respectively, the lagrange solutions in the rotating frame become 
L 3 = -L 2 = (1.19841, 0), Li = (0, 0), L 4 = -L B = (0, 0.866) (See Roy (2004)). 

3.2 The Equal mass cases- A Review 

In this section we discuss equilibrium configurations for the four-body problem 
where the four bodies are of equal mass. These were discussed by Roy and 
Steves (1998). We consider the following cases: A) Four equal masses arranged 
in a square; B) Four equal masses arranged in an equilateral triangle; C)Four 
equal masses arranged in a line. In all the three cases, the origin is taken to 
be the centre of mass of the system and is assumed to be stationary. 

3.2.1 The Equilibrium Configuration of the four body 

problem with four equal masses making a square 

First we find the equations of motion for four particles of mass m; (i = 
1,2, ■• - 4) whose radius vectors from an unaccelerated point C are , while 
the distances between the particles are given by (see any standard text for 
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more detailed discussion, for example Roy (2004)) where 



(3.1) 



Newton laws of motion are then written as 



TO,; r,; 



^ \ TO,TO.j 

G z = l,2,3,4. 



(3.2) 



We choose our units so that G becomes unity. Also let ^ = Pij, then ( 13. 2ft 
becomes 

4 

" (3.3) 



to. 



r'i= G mivrijPij i = 1,2,3,4. 



Equation (13.31) is the form of the equations of motion used for all the derivations 
of equilibrium configurations found in chapter 3. Simplifying equation (13. 3p 
and putting m.j = M, we get the following four equations of motion for the 
equal mass four body problem. 



(3.4) 



r'i= M(p 12 + pn +Pu), 
f 2 = M(p 21 + p 23 + P24), 
r 3 = M(p 3 i +P32 + P34), 

r 4 = M(p 4 i +p 4 2 +P43)- 

Now specific to deriving the square equilibrium configuration, let the four 
particles of mass M lie at the vertices of a square of side of length a, as shown 
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Figure 3.2: Square Equilibrium Configuration of the Four Body Problem 



in Figure (13.21) It is clear from the geometry that 



r 4 = -r 2 , 



ri = -r 3 , 

r 34 — — r 12) 
1*23 = — r 14. 



Also 



ri2 = r 23 = r 34 = r 41 = a. 



(3.5) 



(3.6) 



From Figure (13.21) we can easily find that 



ri 3 = ^24 



b = V2a. 



(3.7) 



Using equations (13 .4p . (13.51) . (13. 6p and (13 .7p . we get the following equations 
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(3. 



of motion for the four-body problem under discussion. 

r 2 =-f(2 + ^)r 2 , 

r 3 =-f(2 + ^)r 3 , 
r 4 =-f(2 + ^)r 4 . 

Since the co-efficient of all the differential equations are the same and negative 
a simple harmonic motion solution with constant angular velocity for all masses 
is possible. We can write all the equations of motion in a compact form as 



follows, 



r- = - n n, i = 1,2,3,4. 



(3.9) 



where 



n 



M . 1 . 



(3.10) 



Equation (I3.9P is a second order linear differential equation with the fol- 
lowing solution 



r i = r jo cos nt H sin nt, 

n 



(3.11) 



where r i0 and rjo are the radius and velocity vectors respectively at t — 
and n is the angular velocity at which the square configuration rotates about 
its centre of mass. 
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3.2.2 Equilibrium Configuration of the four body prob- 
lem with four equal masses making an equilateral 
triangle 

We now consider the equilateral triangle equilibrium configuration of the four- 
body equal mass problem (Roy and Steves, 1998). Let three particles of mass 
M lie at the vertices of an equilateral triangle, with the fourth particle of the 
same mass at the centroid of the triangle. (See Figure ( 13.31) ) As this is an 
equilateral triangle therefore 

r\2 = r 23 = r 31 = a. (3.12) 

Also, as the fourth mass lies at the centroid of the triangle which is also the 
centre of mass of the system, 

r i = r 2 = r 3 = ~~i= and r 4 = 0. (3.13) 
V3 

Using equations (13.41) , (13.121) and (I3.13p . we get the following second order 
differential equations as our equations of motion. 

r 1 =-M( 3 + 3 3 / 2 )r 1 
r 2 =-f (3 + 3 3 / 2 )r 2 
r 3 =-f (3 + 3 3 / 2 )r 3 
r 4 = 

As we have the same co-efficients for all differential equations which are also 
negative, again simple harmonic motion solution is possible. The following is 
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(3.14) 



the compact form of equations (13.261) . 



r,;= - 




(z = 1,2,3) 



(3.15) 



where 




M 



T ( 3 + 33/2 ) - 



(3.16) 



Equation (I3.15P is a second order linear differential equation and has the fol- 
lowing solution 



where r i0 and r i0 are the radius and velocity vectors respectively at t — and 
n is the angular velocity at which the equilateral triangle configuration rotates 
about its centre of mass. 

3.2.3 Equilibrium Configuration of the four body prob- 
lem with four equal masses lying along a straight 
line 

We now consider the collinear equilibrium configuration of the four-body equal 
mass problem. In this case, the solution must be symmetrical about the centre 
of mass C. Then 



Using equations (13 .4p and (I3.18P we get the following equations of motion 



r i0 cosnt H sinnt, 

n 



(3.17) 



r 4 — — r l, r 3 — — 



r 2 and let r 2 = ari. 
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Figure 3.3: Equilibrium Configuration of Four Body Problem making an equi- 
lateral triangle 



r _ M (l , 2(l+a 2 ) . 
1_ \* (1-a 2 ) 2 ) Fl > 



r 2 - (4^2 - (TT^f )) r 2, 



r 3 - -ff (4^2 - (IZ^f )) r 3> 



1-4- (| + frSp) r 4, 



(3.19) 



where R\ = -R4 and i?2 — Rz- 

For a rigid rotating geometry we must have R\ = R 2 = R3 = Ra- Now 



equating Ri to R2 we get 
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6 + 2a 2 



= 0. 



(3.20) 




After further simplification we get 




Therefore 



(a 7 + 6a 5 - a 4 + 25a 3 + 2a 2 - l) = 0. 



(3.22) 



a = 0.3162 is the only real solution for the above equation. This value of a 
gives Ri = 2.966 (i = 1,2,3,4). Hence 



3.3 Solution for two pairs of equal masses 

In this section we discuss equilibrium configurations of two pairs of equal 
masses for the four-body problem. The first two subsections will present a 
review of Roy and Steves (1998) work i.e. the Trapezoidal equilibrium con- 
figuration and Diamond equilibrium configurations, while the last two subsec- 
tions introducing the triangular equilibrium configurations are new analysis 
following the procedure of Roy and Steves (2001). 



r i = r i0 cos nt H sin nt 

n 



(3.23) 



is the required solution, where the angular velocity n is given by: 




(3.24) 
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Figure 3.4: Collinear Equilibrium Configuration of the Four Body Problem 

3.3.1 The Trapezoidal equilibrium configuration for the 

four-body problem- A Review 

We now consider the Trapezoidal equilibrium configuration of the four-body 
problem (Roy and Steves, 1998). It has two pairs of equal masses i.e. M and 
m. In Figure ( 13. 5ft . the geometry is taken to be symmetrical about the line AB 
joining the centre of mass A of the pair P2 and P3, each of mass m, and the 
centre of mass B of the pair Pi and P4, each of mass M. 

Let fi = m/M < 1. Let C be the centre of mass of the system then //r^ + 
r B = 0. Let 

r 23 = -ar u ; \r A - r B \ = f3r u and r = - r B . (3.25) 
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Figure 3.5: Trapezoidal Equilibrium Configuration of Four Body Problem 
We can easily see from the Figure (I3.5I) that 



ri 2 = r 34 = (/5 2 + i(l-a) 2 ) 1/2 r 14 , 
r 24 = r 13 = {[3 2 + \(l + a) 2 ) 1/2 r 14 , 



r B 



(3.26) 



Using the relation given in equations f!3.25j) in conjunction with the geom- 
etry ( 13.26ft . we get the following vectors which describe the locations of the 
four bodies, 



r 2 - TT^ r +f r 4i> 



r 3 



T4 



(3.27) 



44 



The magnitudes of the vectors given in equations (13.271) are the following 



n 

^3 = r 2 , 



1/2 



(3.28) 



r 4 = ri 

The equations of motion of the system are the following 

r'i= M [p (pi2 + pi 3 ) +Pu] , 

f 2= M [p 2 l + PP23+P24] , 

r 3 = M [p 3 i + PP32+P34] > 
r 4 = M [p 4 i + p (P42+P43)] • 
We use equations (13.261) . (13.271) and (13.291) to obtain differential equations 
for r and ri4. 



(3.29) 



M , 

r i4=--r( 2 + P 

'14 



1 — a a + 1 

a 3 b 3 



r 14, 



(3.30) 



where 



and 



(3.31; 



(3.32) 



Now for r we proceed as follows, 



2p 

r 4 + r a = -— — r. 

1 + p 



(3.33) 



Thus 



M 



14 



1 

If' 



(3.34) 
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If the coefficients of r and r 41 are equal and negative in equations (I3.3(jp and 
(I3.34p . then the solutions to these equations give a rigid and rotating geometry. 
Therefore equating the coefficients of r and r 4 i in equations (I3.30P and (13. 34ft 
we get the following 

< 1+ ">(^)- 2 -"(V + i ^)= - < 3 - 35 > 

or 

i±^ + ^-2 = 0, (3.36) 

where a and b are functions defined by equations (13.311) and (I3.32p . 

To obtain a unique solution for a given value of the mass ratio /x, we require 
a second relation in a and (3, which is obtained by using 

r 23 = -ar 14 , (3.37) 

or 



r 3 - r 2 + a(n - r 4 ) = 0. (3.38) 
Differentiating twice and using the equations of motion, it is found that 

- (1 - an) (a + 1) + \ (1 + an) {a - I) - 2a ( 1 - )) r 4i = 0, (3.39) 
b 6 a 6 V a 6 J J 

as r 4 i ^ 0, hence we can write 

(Jp (1 - a//) (a + 1) + ^ (1 + afj) (a - 1) - 2a (l - ^)) = 0. (3.40) 

Solving equations (13.361) and (13.401) simultaneously for a given value of \i (0 < 
/i < 1) we can find a and (3. 

The following two special solutions give the boundary of the family of 
solutions. 
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1. n = 1 (Case of four equal masses) 

Solving equations (13.361) and (13.401) simultaneously for [i — 1 gives a = 1 
and (3 = 1. This is the square solution for the four-body problem. 

2. /i = 0. (Case of the Lagrange solution of the Copenhagen Problem) 

Solving equations ( 13. 36ft and ( 13.401) simultaneously for /i = (i.e. the 
two smaller masses become equal to zero) gives a = and (3 = y^Jl 
and hence the problem reduces to the famous Copenhagen Problem and 
the Lagrange point L 4 is reached. 

Figure (I3.6P gives the locations of the four bodies for a range of /x from 1 
to 0. It shows there is a continuous family of solutions for fi 6 (0,1). The 
bigger points in the figure indicate the locations of the two large masses, which 
remain at size M and the location of the other two bodies when their masses 
are also M (fi = 1). The smaller points indicate the location of the other two 
bodies when /i reduces to zero. As the origin of the system is at the centre of 
mass of the system, the two primaries M appear to move up in the graph to 
their positions at (0.5,0), (—0.5,0) as fi approaches 0. The arrow shows the 
change in the equilibrium solution as /i is varied from 1 to 0, i.e. from the 
square solution to the triangular. 
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Figure 3.6: The evolution of all the four masses when /i is varied from 1 to 0, 
for the Trapezoidal equilibrium solution. (The Origin is located at the centre 
of mass of the system). 



48 




P 2 (m) 



P/M) 



P 3 (M) 



P 4 (m) 



Figure 3.7: The Diamond Equilibrium Configuration of the Four Body Problem 

3.3.2 The Diamond Equilibrium configuration of the 

four-body problem- A Review 

We now consider the diamond equilibrium configuration of the four-body prob- 
lem (Roy and Steves, 1998). It has two pairs of equal masses: a pair of larger 
masses M and another pair of smaller masses m. It is clear from Figure (13. 7p 
that the geometry is symmetrical about the line -P2-P4 and P\Pz- Because of the 
high degree of symmetry only one parameter needs to be introduced, namely 
a, where r 2 = ctr\. Thus 



r 3 = — ri and r 4 = — r 2 . 



(3.41) 
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Therefore only the equations of motion for ri and r 2 need to be studied. 



r'i= M [p (pia + pu) + pi 3 ) 



(3.42) 



r 2 = M [p 21 + P23 + Pp2i] 



After some reduction, by utilizing the symmetry conditions, we obtain 



(3.43) 



M 

1 rf 



M 

r 2 = 3 



1 

4 

p 



2p 



[1 + a 



2^3/2 



ri, 



r 2 - 



(3.44) 



(3.45) 



4a 3 (1 + «2)3/2 

Now proceeding along the same lines as for the Trapezoidal solution and search- 
ing for solutions as the masses m are reduced to zero, we find at p = 1, the 
square solution, and at p = 0, which is the equilateral triangle solution of the 
Copenhagen problem. See Figure (13.81) . Between p = 1 and p = there lies a 
continuous family of solutions. 



3.3.3 Triangular Equilibrium Configuration of four-body 

problem: Case I 

Consider the Triangular Equilibrium Configuration given in Figure (13.91) . where 
two large masses M lie at the two vertices, Pi, P3, of a triangle and two smaller 
masses m lie at the vertex P 2 and at P4 on the line of symmetry of the triangle. 
Let A and B be, respectively, the centre of mass of the pair P 2 and P 4 and the 
pair Pi and P3. Let CA = and CB = r#. From the centre of mass relation 

r B + pr A = 0, (3.46) 
50 




Figure 3.8: The evolution of all the four masses when [i is varied from 1 to 
for the diamond equilibrium solution. (The origin is located at the centre of 
mass, which is also always the halfway point between two primaries.) 
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where 



m 

fi — — < 1. 
M 



(3.47) 



Let 



r\2 = ar 13 and r i4 = fir 13 . 



(3.48) 



We also have from the symmetry conditions that 



ri2 = r 2 3] Ui = r 43 ; and n = r 3 . 



(3.49) 



Now its easy to show that 



n = -/^ - 2^3, 
r 2 = C 2 r A , 
r 3 = -fir A + |r 13 , 
r 4 = C A r A , 



where 



and 



C 4 



(3.50) 



-1 + 4a 2 - /i + 2^a 2 + 2///3 2 - (1 + //) y/(-l+4a 2 ) (-1 + 4(3 2 ) 



2 (a 2 - (3 2 ^ 



(3.51) 



'l - 4/3 2 + /I - 2/ia 2 - 2/i/3 2 + (1 + //) ^(-1 + 4a 2 ) (-1 + A(3 2 ) 



2 (a 2 - p 2 ) 



(3.52) 



Now using the symmetry conditions we get the magnitudes of all the vectors 
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Figure 3.9: Triangular Equilibrium Configuration of the Four Body Problem. 



Case- 1 



involved, which are listed below 

r 2 = C 2 r A , 



r 3 = n 



r 4 = C 4 r A , 



ri2 



ri3 (^-l^) 1 / 2 ^' 

^24 = (C 2 - C 4 ) r A . 
Equations ( 13.501) give 

1 



(3.53) 



(3.54) 



Differentiating equation (I3.54p twice and using equations (I3.29P we get the 
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1.1 



Figure 3.10: Graphs of Fi(a,j3) = and F 2 (a, /3) = for fi = 0.9, showing 



two intersections. 



following equation of motion for 



M 



13 



O." 



/? 3 



i\4, 



(3.55) 



Now further using equations (I3.29j) we get the following equation of motion for 



ri3 



■ 13 



M 



13 



2 + p 



cr 



1 



ri3, 



(3.56) 



Other equations of motion derived in a similar manner, which will be needed 
for further analysis, are the following 



r 2 = 



r 4 = 



M 



// 



(C 2 -C 4 )^(/3 2 -l/4)3/2 



A/(C 4 +aO 



/ X 7 



(C 4 +m) 



' (C 2 - C 4 ) 2 (/3 2 -l/4)3/2 /3 3 

For a rigid, rotating system solution we require 

Ca — Cb = 0, 
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*2, 

r 4 . 



(3.57) 



(3.58) 



which gives 



F 1 {a,/3) = ^| + ^-2 = 0. (3.59) 



For a second relation we proceed as follows. We know from the centre of mass 
relation that 

r A = i(r 2 +r 4 ), (3.60) 

and also 

tA = ^r. (3.61) 

U 2 — O4 

Equating (I3.60p and (I3.6ip and then differentiating the resulting equation twice 
with respect to "time" we get 

(C 4 - C 2 + 2) r 2 +(C 4 -C2-2) r 4 = 0. (3.62) 

Now using equations (13.571) in conjunction with the above equation we obtain 
the second relation needed for the rigid body motion, 



2/i(C 4 + fx) 3 (C 4 - C 2 + 2 ) (C 2 + /i) 

K '2 - C 4 ) 2 (/? 2 - 1/4)3/2 + a 3 

fi{C 4 -C2-2) (C 4 + m) 



F 2 {a,p) = — 2 H = h (3.63) 

' ' " ' 



/? 3 

To find the rigid body solution we need to solve equations (13.591) and (3.63), 
simultaneously, for all values of a and (3. These two equations are highly non- 
linear and therefore its algebraic solution is not possible. We used Mathe- 
matica to find numerical solutions for it. Figure (I3.10p . which is the graph of 
Fi(a, (3) = and F 2 (ct,l3) = 0, given by equations (3.59) and (3.63) respec- 
tively, for n = 0.9 shows that there are two solutions for each \i value. These 
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two (a,/3) solutions correspond to two different triangular equilibrium config- 
urations possible for the given fx. Two special solutions give the boundary of 
the family of solutions: 

1. (/, = l (Case of four equal masses). We get the two well known solutions 
namely the isosceles triangle solution (see figure (13. 11a )) and the equilat- 
eral triangle solution (Simo, 1978), (see figure (13.11b )). The equilibrium 
solutions given by a and f3 are listed in Tables (13.11) and (3.2). We call 
the case 1 solutions that start with the isosceles triangle at fx = 1, solu- 
tion 1, and the case 1 solutions that start with the equilateral triangle 
at (i = 1, solution 2. 

2. (i = (Case of Lagrange solutions of the Copenhagen problem). In case 
1 solution 1-isosceles triangle, as fx — > 0, both P 2 and P4 go to L 4 (see 
Table (13. II) ). In case 1 solution 2-equilateral triangle, as fx — > 0, P 2 goes 
to L 4 and P4 goes to L\ (see Table (3.2). 

Between fx = 1 and 0, there is a continuity of solutions, each fx having 
two equilibrium solutions. Please note that there are some numerical errors in 
table 3.1 and 3.2 of the order of 10~ 7 as fx goes to 1. These are given in Tables 
(13.11) and (3.2) and two families of the equilibrium configurations are shown in 
Figure (3.11a) and (3.11b). Note that Figures (3.11) have as their origin the 
point halfway between the two primaries of masses M. Thus unlike in Figure 
3.6, the centre of mass is a point that moves as fx is reduced from 1 to 0. 
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Table 3.1: Equilibrium solutions in the triangular case 1, solution 1, Isosceles 
triangle. 



u 




1 1 


n = 


74 


r 2 


1 
1 


l.UO 


U. Oc/U 


U.OOUO lOOl 


U.UIOIOZUO^ 




n q 


1 1Q 
1. io 




n ^7i n^Q^ri9 


fl n7Q^zLn9DQ 
u.u 1 yo^Uiuy 


U.U 1 ^Z^tUUJU 


0.8 


1.17 


0.683 


0.545071889 


0.126825011 


0.719321922 


0.7 


1.196 


0.709 


0.518599486 


0.175497276 


0.759292937 


0.6 


1.213 


0.732 


0.492688262 


0.227166212 


0.797696755 


0.5 


1.22256 


0.754 


0.46877787 


0.284370156 


0.835638092 


0.4 


1.224 


0.776 


0.448824987 


0.349063876 


0.872837665 


0.3 


1.217 


0.799 


0.436231551 


0.423284096 


0.909610336 


0.2 


1.199 


0.825 


0.436287216 


0.510720961 


0.944271801 


0.1 


1.162 


0.861 


0.455927961 


0.621403477 


0.969385705 


0.01 


1.074 


0.932 


0.494662843 


0.777927333 


0.941914337 


0.001 


1.034 


0.967 


0.499478426 


0.826835512 


0.904206852 


0.0001 


1.015 


0.985 


0.499948997 


0.84857412 


0.883216869 


0.00001 


1.007 


0.993 


0.499994953 


0.857924321 


0.874090305 
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1.2 - 




-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 



(a) 




-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 



(b) 

Figure 3.11: The evolution of all four masses when /x is varied from 1 to zero in 
the triangular equilibrium case-1 (a). Solution 1-Isosceles triangle (b) Solution 
2- Equilateral triangle (The Origin is located halfway between the two primaries 
and thus the centre of mass moves as \i is varied.) 
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Table 3.2: Equilibrium solutions in the triangular case 1, solution 2-Equilateral 
triangle 



fj, 


«2 


P2 


7"i 


r 4 


?"2 


1 


1 


0.57735 


0.577350202 


4.03 x 10 u ' 


0.577350404 


0.9 


0.910515913 


0.536533511 


0.496437204 


0.031715456 


0.53463326 


0.8 


0.874052 


0.521807 


/~\ A -1 t~\ /~\ /~\ ^7 -I -I -I 

0.418007111 


0.043212884 


0.524428572 


0.7 


0.851383 


0.513011 


0.328494841 


0.050703648 


0.523586719 


0.6 


0.837823 


0.507424 


0.222436467 


0.055784156 


0.53000428 


0.5 


0.832221 


0.503911 


0.122366589 


0.058662486 


0.543953079 


0.4 


0.834933 


0.501822 


0.236552793 


0.05890323 


0.567038183 


0.3 


0.847573 


0.500701 


0.589645065 


0.055537618 


0.602359626 


0.2 


0.84101 


0.5002 


1.128549593 


0.043388245 


0.618706196 


0.1 


0.7873491 


0.500194 


0.414519964 


0.01434931 


0.579930294 



3.3.4 Triangular Equilibrium configuration of four-body 

problem: Case II 

We now consider the Triangular Equilibrium Configuration given in Figure 
(I3.12p . The geometry is taken to be symmetrical about the line P^B. In this 
case it is the masses P2 and P4 that retain the original mass M with the 
masses Pi and P3 decreasing in mass though remaining equal in mass m so 
that fi=^<l. 

Points A and B are the centers of mass of the pair P2 and P4 and the pair 
Pi and Pi- Then letting ya = CA and r# = CB, we have 

r A + fiv B = 0. (3.64) 
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P,(m) 



Figure 3.12: Triangular Equilibrium Configuration of Four Body Problem. 
Case- 1 1 



Let 



r = P2-P4 = i"24, r_B = CB = ar and P\B = /3r, 



(3.65) 



(3.66) 



then from the centre of mass relations and the geometry of the figure, we get 
the following list of vectors representing all the four bodies 

ri = ar + ~r 31 , 
r 2 = - (| + fjux) r, 
r 3 = ar - |r 31 , 
r 4 = (| — fiaj r. 

Using Figure (13.121) and the symmetry conditions the magnitudes of the 
vectors involved are 

n =r 3 = (^Ja 2 + /3 2 ) r, 
r 2= {\ + r, 
r 4 = (| - /ia) r, 
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(3.67) 



Also 



Tl2 



ri4 



^(I + a(l +/ x)) 2 + ^j r = r 23 , 
^(4 + a(l + /i )) 2 + ^ r. 



(3.68) 



Using the same techniques as before we get the following differential equations 
as equations of motion, 



M 



2ar 3 
Xr, 



l + 2a(l + /i) 



1 - 2a (1 +fi) 



((i + a (1 + ^)) 2 + ^) 3/2 + a (i + At )) 2 + pi} 



3/2 



(3.69) 



r 3 i 



M 



Yr 31 . 



+ 



_ ((i + a (1 + /i)) 2 + /3 2 ) 3/2 + a (1 + Z^)) 2 + /? 2 ) 



+ 



3/2 4/3 3 



r 3 i, 



(3.70) 



For a rigid body motion we must have X — Y = 0, which gives us 



jia + 1 /2 



Ha-1/2 



((§ + a (1 + /i)) 2 + /5 2 ) 3/2 ((-§ + a (1 + /i)) 2 + ^ 



3/2 4/3 3 
(3.71) " 



To completely find the rigid body solution, we need a second relation for 



which we make use of the following 



, r 2 + r 4 

= — liar and ya = — - — , 



(3.72) 



this gives us 



(fia + 1/2) r 4 - (fia - 1/2) r 2 = 0. 



(3.73) 



Differentiating equation (I3.73P twice with respect to "time" and then after 
some simplifications we get the following 



G 2 (a,ff) 



1 ^^(1 + ^ + 1/2^ 



(3.74) 
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1 \ ( g{l + n)-\/2 \ 
2 +fJX *) { b ) - « = 0, 



where 



and 



A=[[- + a(l + n)) +(3 



3/2 



B 



((4 + «(l + /i)) 2 + /? 2 j 





(3.75) 



(3.76) 



(b) 

Figure 3.13: Graphs of the functions Gi(a,P) = and G 2 (a,P) = for a) 
fi = 0.998 b) \i = 0.9971 

To find the equilibrium solutions we need to solve equations (13.711) and 
(3.74), simultaneously, for values of a and /3, given a range of /i from 1 to zero. 
These two equations are highly non-linear and cannot be solved analytically. 



62 




Figure 3.14: Evolution of all four masses when /x varies from 1 to zero in the 
second triangular case (The Origin is located at the Centre of mass) 

Therefore we find numerical solutions using Mathematica. There are two so- 
lutions for each /i value, see figure (13. 13a ) for the example when \i = 0.998. 
These solutions are very close to each other as they have a difference of order 
1CT 2 at most. Figure (13.131) is the graphs of G\(a, (3) = and G^a, 0) = for 
H = 0.998 and \i = 0.9971. Figure (13.13b ) shows that there is no solution for 
/i < 0.9972. The solutions between fi = 1 and /i = 0.9972 are shown in Figure 
(I3.14p . At /i = 1 we get the equilateral triangle solution, see Figure (13.141) . 
as we would expect. In fact case-ii is a continuation of the family of solutions 
found in case-i, where /i in case-i goes to 1/0.9972. 
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Table 3.3: Equilibrium solutions in the triangular case 2 





a 


P 


r\ 


r-i 


r 4 


1 


0.5 
0.528 


0.866 
0.857 


0.999978 
1.006 


1 

1.028 




0.028 


0.999 


0.503 
0.526 


0.865 
0.857 


1.0006 
1.0055 


1.002497 
1.0254 


0.0024 
0.025 


0.998 


0.50757 
0.522 


0.863 
0.859 


1.0012 
1.0052 


1.0066 
1.021 


0.0066 
0.021 


0.9975 


0.519 
0.511 


0.86 
0.8627 


1.0045 
1.0027 


1.0177 
1.0097 


0.0177 
0.0098 


0.9974 


0.5122 
0.518 


0.8624 
0.8605 


1.0031 
1.0044 


1.0109 
1.0167 


0.01087 
0.01665 


0.9973 


0.5158 
0.5144 


0.8612 
0.8617 


1.0039 
1.0036 


1.0144 
1.0130 


0.0144 
0.0130 


0.99725 


0.51272 
0.5174 


0.8622 
0.8607 


1.003 
1.00424 


1.0113 
1.0160 


0.0113 
0.01598 


0.9972 


0.5004 
0.5035 


0.8662 
0.865 


1.00035 
1.0009 


0.999 
1.0021 


0.001 
0.0021 
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Figure 3.15: Collinear Equilibrium Configuration Case-I 

3.3.5 Collinear equilibrium configurations for two pairs 

of equal masses- A Review 

In this section we discuss all possible arrangements of the two pairs of equal 
masses along a straight line (Roy and Steves, 1998). The first two of the 
four arrangements are symmetric while the last two arrangements are non- 
symmetric. 

Case-I 

In this symmetrical arrangement of two pairs of different masses, the larger 
masses (M) lie in the middle of the line and the smaller masses (m) lie at the 
corners apiece, as shown in the Figure (13.151) . 

This arrangement is symmetric about the center of mass C. 
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Let r 2 = ari, ji = m/M < 1 and p^- = r^/r 3 -. By symmetry r 4 = — r x and 



r 3 = -r 2 . 

Now using the general equations of motion in conjunction with the above 
considerations we get, 

ri= M(pi2 + P13 + PPm), 

r 2 = M(p2i +P23 + ^24) ■ 
Using all the symmetry conditions we obtain the following final form of the 

equations of motion 



(3.77) 



M 

r l— JJLlI - !, 



(3.78) 



where 



Also 



Ri 



+ 



+ 



l -ay {1 + ay 4 



r "2— jR-2 r 2, 



(3.79) 



(3.80) 



where 



1 u / 1 
+ ' 



4a 3 ' a V(l + «) 2 (1 -«) S 
For a rigid rotating geometry we must have 



(3.81) 



R 1 -R 2 = 0. 



(3.82) 



Therefore 



1 1 

+ 



1 fi ( 1 



(1-a) 2 (1 + a) 2 4a 3 a\(l + a ) 2 (1 - a) 2 J ' -J 
After further simplification we get 

rr , . 2 (1 + a 2 ) + 4u If 1 \ 

HA a) = — , — T -f- =0. 

U ; (l-« 2 ) 4 V «V 



+ ^ = 0. (3.83) 



(3.84) 
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0.36 



0.34 
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Figure 3.16: Variation of the parameter a for all values of /x in Case-I of the 
collinear equilibrium configurations of two pairs of equal masses 
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Solving the above equation for p = 1 (the equal mass case) and p = (i.e. 
m = 0) we get a = 0.44 and a = 0.316 respectively. All other solutions for 
different mass ratios can be found easily. See figure (I3.16p . which is graph of 
Hi(a) and shows there is a continuous family of solutions for p e (0, 1). 



In this second symmetrical arrangement of the two pairs of different masses, 
the two smaller masses are near the centre of mass while the two larger masses 
are farther away. 

Again proceeding along the same lines as before, let r 2 = ar±. By symmetry 
r 4 = — i"i and r 3 = — r 2 . The mass ratio is /i = m/M < 1 and pij = Tij/rfj. 
The equations of motion are 



Using the symmetry conditions, we get the following simpler form of the 
equations of motion 



Case-II 



r 1= M(/ipi2 + /ipi3 + Pu) 



(3.85) 



r 2 = M(p 21 + pp 2 3 + p24>- 



r i— Q"-Ri r i ; 



(3.86) 



where 




(3.87) 



and 



r 2— o"-R2l"2, 
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Figure 3.17: Collinear Equilibrium Configuration Case-II 

where 

^((TT^-(I^F + ^)- (X89) 

For a rigid body motion we must have 

R 1 -R 2 = o. (3.90) 

Hence 

K (T^ + (T^) + ^ K (IT^f ~ (T^f + 4&) =0 ' (3 ' 91) 

After further simplification we get 

*■<<•) = i (i - £) + tt^< 2 + "f 1 + Q ' 2 » = °- < 3 - 92 > 
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Figure 3.18: Variation of the parameter a for all values of /x in Case-II of the 
collinear equilibrium configurations of two pairs of equal masses 
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Figure 3.19: Collinear Equilibrium Configuration Case-Ill 

Solving equation (13.921) we get a = 0.275 for /i = 1, which is the equal 
mass case also obtained for case 1, \i — 1. There is no solution for fi = but 
as fi approaches zero a converges to zero which is the L\ Lagrange solution. 
See figure (13.181) . which is the graph of H 2 (a) and shows there is a continuous 
family of solutions for /i G (0, 1). 

Case-Ill 

From Figure (13.191) we can see that this is a non-symmetric arrangement of 
the four bodies. All arrangements previous to this, have had some kind of 
symmetry which has considerably simplified the problem. 
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Let 

r 2 = ari; r 3 = -/?n; r 4 = -7 
By the centre of mass relation 



^rriiTi = 0. 



i=l 

Thus 

mr 1 + Mr 2 + mr 3 + Mr 4 = 
After substituting /i = m/M we get, 

^ (ri + r 3 ) + r 2 + r 4 = 0, 

from equations (13.931) and (13.961) we get 

(fi (1 - P) + a - 7) rj = 0. 

As ri 7^ therefore 

// (1 - /3) + a - 7 = 0, 

gives 

7 = // (1 — 0) + a. 

The equations of motion are 

r'i= M [p 12 + PP13+P14] > 

r 2 = M [//p 2 l + PP23+P24] ) 

r 3 = M [fip 31 + P32+P34] > 
r 4 = M [/xp 4 i + p 42 +pp 43 ] , 

where p - = 
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Now using equation (I3.93P and the centre of mass relation, it can easily be 
shown that 

M , 
r 1= -— R x x u (3.101) 



where 



Also 



where 



and 



where 



and 



where 



n 



1 a 1 , 

i?i = ? + — + K. 3.102 

(1-af (1 + /3) 2 (l+ 7 ) 2 



M 

r 2 = 3 R 2 r 2 , (3.103) 



R * = £ (7-^2 + 77—71 ) + -7-^—7^ ( 3 - 104 ) 



(a + j3) (1 — a) J a (a + 7) 



M . . 

r 3 = -— R 3 r 3 , (3.105) 



M 

r 4= 3-i?4r 4 , (3.107) 



7 ( A *((l + 7 ) a + ( 7 - i 9)V ' (7 + ar 



^4 = - ( H ( 7— —2 + 7 772 + 7—^2 ) • ( 3 - 108 ) 



For a rigid rotating solution we must have R\ = R 2 = R3 = -R4 > i.e. 
Ri = R > and constant. We therefore require pi = Ri — Ri = 0, i = 2, 3, 4 for 
the set of values of a, (3 and 7. Now from equation (13.991) each relation Pi = 
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Figure 3.20: Collinear Equilibrium Configuration Case-IV 

is a curve in the plane. It is easy to show that if any two of the p^s are 
zero the third will also be zero (see Roy and Steves (1998)). Thus to find a 
solution we need to solve P2 = and ps = in the a, (3 plane to obtain the 
values of a and (3 and hence 7. When the masses of Pi and P 3 are reduced to 
zero they migrate to the L 3 and L\ Lagrange point respectively. The third P 4 
of mass M migrates to the point 1/2 with Pi remaining at —1/2. 

Case-IV 

This is the second and last non-symmetric case of four collinear configurations 
under discussion. This will be treated the same way as case-Ill. 
Let 

r 2 = axi;r 3 = -/3ri;r 4 = -7^. (3.109) 
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From the centre of mass relation 



ri + r 2 + p(r 3 + r 4 ) = 0, 



(3.110) 



where p = m/M. Now from equations (13. 109[) and (I3.110P we get 



a = p (/3 + 7) — 1. 



(3.111) 



The equations of motion are 

r'i= M [pi2 + p (P13+P14)] » 
r 2 = M [p 2 i + p (P23+P24)] , 
r 3 = M [p 3 i + P32+PP34] , 
r 4 = M [p 4 i + P42+PP43] • 
It can easily be shown, as in the previous case, that 



(3.112) 



(3.113) 



where 



Also 



1 / 1 1_ 

1 ~ (1 - a) 2 + ^ {(1+P) 2 + OTjf 



2 / • 



(3.114) 



M 

r 2— jR2 r 2, 



(3.115) 



where 



and 



R*=- 



1 / p 



a. 



{a + PY (I (a + 7)' 



M 

r 3— 3 -fL3 r 3) 



2 / ' 



(3.116) 



(3.117) 
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where 

R & = \( + — * ^— 9 I , (3.118) 

and 

M . . 

r 4 = 3^4, (3.119) 

where 

IU = -( + + — 9 ] • (3.120) 

7V(l + 7) 2 (7-/5) 2 (7 + «) V 

For a rigid rotating solution we must have R\ = R2 = R3 = R4 > i.e. 
Ri = R > and constant. We therefore require pi = Ri — R± = 0, i = 2, 3, 4 
for the set of values of a, /3 and 7. From equation (13.1111 ) each relation pi = 
is a curve in the /3, 7 plane. We will solve the following two equations for j3 
and 7 to obtain the solution (a,/3)which gives the rigid rotating geometry. 

R 3 - r x = and R A - R 1 = 0. (3.121) 

The two smaller bodies migrate to the L2 Lagrange point when /1 — > and 
at = 1 gives the collinear equal mass solution. 

3.4 Summary and Conclusions 

In this chapter we discussed special analytical solutions for the coplanar four 
body problem of equal and non-equal masses (Roy and Steves, 1998). The 
non-equal mass cases have two pairs of equal masses where the ratio between 
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the two pairs is reduced from 1 to 0. Three special arrangements are discussed 
for the four equal masses which include: 

1. Four equal masses making a square. 

2. Four equal masses arranged at the vertices of an equilateral triangle with 
the fourth mass at the centroid of the triangle which is also the centre 
of mass of the system. 

3. Four equal masses lying along a straight line symmetric about the centre 
of mass. 

In section 3.2 we allow two of the masses to reduce symmetrically to m 
while the other pair of bodies remain at the original mass M. We define 
pi = m/M < 1. In the cases studied by Steves and Roy (1998) they found a 
continuous family of equilibrium solutions occurred as fi was reduced from 1 
to 0. In all, for /i — 1, the equal mass case solutions were obtained which is 
a good check that the equations giving the families of solutions are correct as 
/i — 1 is a special case. Then in the limit, as /x — > 0, they always obtain one of 
the Lagrange five equilibrium points Li,L 2 ,L 3 , L 4 and L 5 of the Copenhagen 
problem (Roy and Steves, 1998) as locations for the two masses being reduced 
to zero. We complete the analysis of Roy and Steves (1998) to include two 
more examples of equilibrium solutions of the four body problem i.e. 

1. The Triangular equilibrium configuration of the four body problem with 
the two bigger masses making the base of the triangle. 
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2. The Triangular equilibrium configuration of the four body problem with 
the two smaller masses making the base of the triangle. 

In the Triangular Case-I for fj, — 1 we get the two well known equal mass 
solutions i.e. the equilateral triangle solution and the isosceles triangle solution 
and as /x is reduced to zero the isosceles triangle evolves so that P 2 and P 4 
approaches L 4 . The equilateral triangle solution evolves so that P 2 approaches 
L 4 and P 4 evolves to L\. In the Triangular Case-II for /i — 1, the equilateral 
triangle solution is obtained. There are no solutions for /i < 0.9972. This is 
the only case where no continuous family of solutions exist from /i — 1 to 0. 

We conjecture that the equilibrium solutions discussed in this chapter are 
all the equilibrium solutions which exist for the symmetric case where two of 
the masses reduce symmetrically to m while the other pair of bodies remain 
at the original masses M. 
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Chapter 4 



The Stability Analysis of the 
Nearly Symmetric Caledonian 
Symmetric Four Body Problem 
(CSFBP) 

In this chapter we study the stability of a more general case of the four 
body problem called the Caledonian Symmetric Four Body Problem (CSFBP) 
(Steves and Roy 2001). 

To tackle the complicated nature of the general four body problem, different 
restriction methods are used, the neglecting of the mass of some bodies being 
the most common (see chapter 1 for some examples). Another method of 
restriction, not as common but very effective, is the introduction of some 
symmetry conditions (Steves and Roy, 2001). This type of restriction method 
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reduces the dimensions of the phase space very effectively, whilst still producing 
a model which can be very close to real systems. 

The Caledonian Symmetric Four Body Problem (CSFBP) is a four body 
system with a symmetrically reduced phase space. This symmetrically re- 
stricted four body problem was developed by Steves and Roy (1998) and later 
they derived an analytical stability criterion for it (Steves and Roy , 2001). 
After 2001 different aspects of this problem were studied by Szell (2003), Szell 
et al. (2002) and Szell, Erdi, Sandor and Steves (2004). 

In this chapter we investigate the stability of the symmetric nature of the 
CSFBP by using nearly symmetric, slightly perturbed, initial conditions and 
the general four body equations to be able to see if the CSFBP system re- 
mains nearly symmetric. Covering a comprehensive range of initial conditions, 
we integrate each orbit for a million time-steps, to determine the evolution of 
the system, recording when and if its symmetry is broken. An integrator is 
specifically developed for this purpose using the Microsoft Visual C++ Soft- 
ware. The results of integrations are processed using Matlab 6.5. During the 
integration we record the following observation. 1. We stop the integration 
when there is a close encounter and record the type of the collision or close 
encounter which we color code on the graphs shown. 2. We stop the inte- 
gration when it fails the symmetry breaking criterion and color code it to be 
shown on the graph. 3. If there are no collision or symmetry breaking then 
the integration continues until the 1 million time-steps which we also note and 
color code it for the graphs. 
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In section 4.1 we introduce the CSFBP which is a review of Steves and Roy 
(2001). In section 4.2 we derive the equations of motion for the general four 
body problem to be used for integrating the CSFBP. In section 4.3 we discuss 
the initial conditions of the CSFBP. Section 4.4 discusses the general four body 
integrator, specifically developed for this analysis by the author to integrate 
the general four body problem. A comprehensive set of orbits covering the 
initial phase space of the CSFBP are integrated in order to determine their 
state of symmetry after a long time. In section 4.5, the procedure of analysis 
is explained. In section 4.6 we discuss the results of the different integrations 
performed for the equal mass case of the CSFBP. These results are compared 
with that of Szell, Erdi, Sandor and Steves (2004) in section 4.7. In section 

4.8, the results of the different integrations performed for the n — 0.1 case of 
the symmetric and nearly symmetric CSFBP's, are analysed. These results 
are compared with those of Szell, Erdi, Sandor and Steves (2004) in section 

4.9. Finally in section 4.10 conclusions to the chapter are given. 



4.1 The Caledonian Symmetric Four Body Prob- 
lem (CSFBP)- A Review 

Steves and Roy (1998, 2000, 2001) have recently developed a symmetrically 
restricted four body problem called the Caledonian Symmetric Four Body 
Problem (CSFBP), for which they derive an analytical stability criterion valid 
for all time. We will give its brief introduction in this section. 
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Let us consider four bodies Pi,P 2 ,P3,P<i of masses m 1 ,m 2 ,m 3 ,m4 respectively 
existing in three dimensional Euclidean space. The radius and velocity vectors 
of the bodies with respect to the centre of mass of the four body system are 
given by and r-j respectively, i = 1,2,3,4. Let the centre of mass of the 
system be denoted by O. The CSFBP has the following conditions: 

1. All four bodies are finite point masses with: 

mi = r«3 = M, 777,2 = W4 = 171 (4-1) 

2. Pi and P3 are moving symmetrically to each other with respect to the 
centre of mass of the system. Likewise P 2 and P4 are moving symmetri- 
cally to each other. Thus 

ri = -r 3 , r 2 = -r 4 
Vi = h = -h, V 2 = r 2 = -r 4 , (4.2) 

This dynamical symmetry is maintained for all time t. 

3. At time t = the bodies are collinear with their velocity vectors perpen- 
dicular to their line of position. This ensures past-future symmetry and 
is described by: 

n x r 2 = 0, i"! • fx = 0, r 2 ■ r 2 = (4.3) 
Figure (14.1b ) gives the initial configuration of the CSFBP. 
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(b) 

Figure 4.1: a. The initial configuration of the CSFBP b. CSFBP for 
t > 
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4.2 The Equations of motion 

The classical equations of motion for the general n-body problem are given by 

-y-r v , . = 1,2,3... (4.4) 



rrii rf- 

— ' r v 



where ri = {x^yi) and = r j — ri. For a general four-body problem we will 
get the following equations of motion from (14.41) above 



m 2 r 12 , m 3 ri 3 m^r u 
r 1= — 3— + — 5— + — 3— 4.5 

'12 '13 '14 

mir 21 m 3 r 23 , m 4 r 24 , . 

r 2 = — 3— + — 3— + — 3— (4.6) 

'21 '23 '24 

mir 31 m 2 r 32 m 4 r 34 
r 3 = — 3— + — — + — 3— (4.7) 

31 '32 ' 34 

mir 41 m 2 r 42 m 3 r 43 
r 4 = -A^ + -A^ + -4^. 4.8 

ftiO /pO /y^O 

'41 '42 '43 

Let n = jj i.e the smaller mass divided by the larger mass. We take M — 1 
and therefore /1 = m. Equations (14. 5 p to (14. 8p become 



/ i"i2 . r 14 \ r 13 

ri= // — + — + — (4.9) 

V r 12 r 14/ r 13 

r 21 . r 23 . A tr 24 /a tn\ 

r 2 =^- + ^- + — 3- (4.10) 
'21 '23 '24 

r 31 . ( r 32 . r 34 \ //i 1 1 \ 

r 3=ir + /^L3- + i3- ( 4 - n ) 

'31 V ' 32 ' 34 / 

r 4 =^ + ^ + ^. (4.12) 

r 41 r 42 r 43 

The above equations can be further simplified by using all possible symmetries 
of the CSFBP (Szell, 2003) but we will not use the symmetric equations of 
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motion as we want to use the most general form of the equations of motion to 
be able to see what happens when we marginally break the symmetry. 
The potential function may be written as 

U = Gtt^ M^<<, (4-13) 

i=l 3=1 *J 

where 



Then the energy E of the system may be written as 

E = T-U (4.14) 

where T is the kinetic energy given by 

I 4 

i=i 

4.3 The Initial Conditions 

The numerical integration of (14.91) to (14.121) require initial values for ri, r 2 , r 3 
and r4. In order to satisfy the initial conditions of the CSFBP we immediately 
have 

r x = (xi, 0) = -r 3 and r 2 = (x 2 , 0) = -r 4 (4.16) 

We provide the values of X\ and x<i by hand by varying x\ and x<i in the 
intervals (0, 1.5] and (0, 1.5] with a step-size of 0.008. The initial velocities 
Vi of P\ and V 2 of P 2 are calculated using the following relations which are 
derived using the energy and angular momentum equations. 
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B - VB 2 - AAC 
~2A 



Viv = — > (4-17) 



where 

2 



A = l + -[ X -1) (4.18) 



5 = -2* (4.19) 
(/xx 2 ) 

C=^-U-Eo (4.20) 



— -M^-i (4.21) 

Where V\ y and V^j, are the y components of Vi and V 2 respectively, c = 
J is the angular momentum of the system, E'o is negative of the energy E 
and U is the potential given in section 4.2. The x components of Vi and V 2 
are set to zero for t = 0. Co is the Szebehely constant (Steves and Roy, 2001). 
More detailed analysis of how Cq can be used to determine the stability of the 
four body problem can be found in chapter 5. But for now, we simply set Co 
to be key values as used by Szell, Erdi, Sandor and Steves (2004). 

Our aim in this chapter is to examine the motion and stability of a variety 
of the CSFBP systems each with a different Szebehely constant C . We also 
want to compare systems with the same Szebehely constant. Thus we select 
the initial conditions so that they result in the same Cq. For a given Co, the 
variables V\ y and can be calculated from equations (14. 17ft and (14.211) . 

We investigate the following sets of mass ratios of the CSFBP with several 
values of C for each set of mass ratios. 
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1. fl= 1 

(a) C = 40, 46, 60 

2. n = 0.1 

(a) C = 0.6,0.5,0.9 

4.4 The integrator 

The equations of motion for the CSFBP, equations (14. 9 j) to (14.121) . are highly 
non-linear second order coupled differential equations. It is not possible to find 
their analytical solution. Therefore we have to use some numerical technique. 
To do so we need to develop an integrator which is a software in which the input 
data are the initial parameters of the differential equations and the output 
data are the solution of differential equations belonging to the input data. To 
develop the integrator we first need to find an accurate and fast numerical 
method and then develop an environment where the input and output data is 
efficiently handled. 

The numerical method, we chose, for this integration project is a 15th order 
method with a adaptive step size control called the Radau method of Everhart 
(Everhart, 1985). This method makes use of Gauss- Radau spacing. It varies 
the time-step used according to the rate at which the variables are changing 
and thus it makes an efficient integrator that deals well with close encounters. 
We used Microsoft Visual C++ 6.0 to develop the integrator for the general 
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four body problem and to construct the environment of the integrator. This 
integrator can easily be generalized for n body problems with n being an even 
integer. 

4.5 Procedure of Analysis 

Szell (2003) discussed the chaotic and regular structure of the phase space of 
the Caledonian Symmetric Four Body Problem. They gave some very interest- 
ing results about the relationship between the global hierarchical stability of 
the CSFBP and the chaotic and regular nature of its phase space. We further 
their analysis by perturbing the symmetric initial conditions and using the 
general four body equations of motion. We integrated the system with nearly 
symmetric initial conditions as well as symmetric initial conditions, while using 
the general four body equations. We studied two indicative CSFBP systems 
stipulated by /i — 1 (a quadruple stellar system) and ji = 0.1 (a binary star 
system with two very massive planets or brown dwarf stars). For each given 
mass ratio, //, we kept the energy of the system constant, effectively providing 
a relative scale of size for the system and then chose a range of increasing 
Szebehely constants (Steves and Roy, 2001). 

For any given /i, C and E each CSFBP orbit is uniquely determined by its 
initial values r\ and r 2 . Therefore the nature of each orbit can be depicted by 
integrating it for some reasonable time. For each orbit, we follow its evolution 
searching for collision or the breaking of the symmetry of the problem. Here, 
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'collision' is denned to be a close encounter between two of the bodies such that 
the conservation of energy fails i.e. there is a difference between the energy 
at t — and the energy at the current time of the order greater than 1CT 15 . 
The symmetry is defined to be broken when {x\ + x 3 ) 2 + (y 1 + y 3 ) 2 > 1CT 4 or 
(x 2 + x 4 ) 2 + (y 2 + y 4 ) 2 > 1CT 4 . We constructed r\r 2 graphs giving different 
colors to each type of orbit. See table 5.1 for a listing of the different color 
coded final events possible. We have colored the orbit black if symmetry is 
broken, red if there is a collision between Pi and P 2 (12 collision), yellow if 
there is a collision between P 1 and P 3 (13 collision), magenta if there is a 
collision between Pi and P 4 (14 collision), blue if there is a collision between 
P 2 and P 3 (23 collision), green if there is a collision between P 2 and P 4 (24 
collision), cyan if there is a collision between P3 and P4 (34 collision) and grey 
if nothing happens during the whole process of integration. We will refer to 
such orbits as stable. The collisions here are not physical collisions. It is the 
situation when any two of the four bodies come very close to each other and the 
energy conservation fails. In such situation we stop the integration and label 
the event as a collision orbit. Please note that we cannot conclude that such 
orbits are stable as we do not know what will happen to them if we integrate it 
for longer. Also it is not possible to conclude that orbits which do not remain 
symmetric are unstable as it is possible to have stable orbits which are not 
symmetric. 

For fj, — 1 for each (Co, E ) we produce nine rir 2 graphs which include: 

1. Three r\r 2 graphs with symmetric initial conditions i.e. 
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(a) 10 4 integration time 

(b) 10 5 integration time 

(c) 10 6 integration time 

2. Three rir 2 graphs with a perturbation of 10~ 6 in the x component of Pi 
to the symmetric initial condition (a,b,c same as above). 

3. Three rir 2 graphs with a perturbation of 10~ 5 in the x component of Pi 
to the symmetric initial condition (a,b,c same as above). 

For /j, — 0.1, for each (C , E ) we produce six rir 2 graphs: 

1. Three r\r 2 graphs with symmetric initial conditions (a,b,c same as above). 

2. Three r\r<i graphs with a perturbation of 10~ 5 in the x component of Pi 
to the symmetric initial condition (a,b,c same as above). 
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Table 4.1: Criterion and color codes for the different categories of orbits 



Criterion 


Nature of Orbit 


Color code 


Sundman inequality not true 


Forbidden region to real motion 


white 


(xi + x 3 ) 2 + ( Vl + y 3 ) 2 > or 


Non-symmetric 


Black 


(x 2 + x 4 ) 2 + (y 2 + y 4 ) 2 > 






ri2 ~ 


12 type of collision 


red 


r 13 w 


13 type of collision 


yellow 


ri4 w 


14 type of collision 


magenta 


r 23 w 


23 type of collision 


blue 


r 24 w 


24 type of collision 


green 


r 34 ~ 


34 type of collision 


cyan 




stable 


grey 
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4.6 Results: The nature of orbits in the 

space for the equal mass case of the CSFBP, 

fi = 1 

The value of E was fixed to be -7 and Co was chosen to be 40, 46 and 60. The 
initial values of r\ and r 2 were varied from to 1.5 with a step size of 0.008. 

Recall that we use the general four body integrator and it's not necessary 
to have symmetric collision. For example a 12 collision does not necessarily 
mean a 34 collision and vice versa. To have guaranteed symmetric collisions 
one must use the equations of motion of the CSFBP (Szell, 2003). For color 
coding of the different types of orbits see table 14.11 

In this equal mass case, an interchange of the r\ and r 2 produces the same 
physical orbit. Thus the categorization of the orbits is symmetric with respect 
to the rir 2 line with the exception that the 24 and 13 collisions are reversed. 
The initial ordering of the equal mass case for r\ ^> r 2 is 1243 which becomes 
2134 for r 2 ^> r±. Therefore the resulting categorization will be the same 
except that some of the collisions will be reordered. Thus comparing the 
collision types from r\ ^> r 2 space to r 2 ^> r\ space, 24 collisions becomes 13 
collision and vice versa. 12, 14, 34 and 23 collisions remain the same. 

We can separate the graph into three regions: The double binary region 
around r\ ~ r 2 (DB) and two single binary regions around r 2 ^> r\ (SB1) and 
n >r 2 (SB2). 
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Figure ( 14. 3ft shows the results of integrations for Co = 40 with symmetric 
initial conditions. Figure (14. 4D and (14. 5 h shows the results with perturbed 
initial conditions with a perturbation of 10~ 6 and 10~ 5 respectively in the x 
component of Pi. 

What is immediately clear from the comparisons of figure (14.3b ). (14.3b ) and 
(14.3b ). is that the longer we integrate this CSFBP system the more unstable 
orbits we unearth. The most chaotic region for Cq = 40 are the single binary 
regions, as most of the orbits end up in collisions after a very small integration 
time, see figure (14.3b ). Most of the collisions in the SB1 region are 24 type 
collisions. There are also some 12 and 14 type collisions. The orbits in this 
region always start in the 24 type hierarchy state which is the reason for most 
of the orbits being 24 type collisions. In the SB2 region we have the same 
number of 13 type collisions as we had of the 24 type collisions in SB1. The 
double binary region is surrounded by 12 and 34 type collisions. See Figure 

Figure (14.3b ) shows the results for the same problem but for a longer inte- 
gration time i.e. 10 5 . The few orbits in the SB1 and SB2 region which survived 
the shorter version of integration have ended up in collisions except for a small 
grey island at the bottom i.e. close to the origin. This is a clear indication 
of the chaotic nature of these regions as opposed to the double binary region 
which is mostly grey. Again most of the collisions in the SB1 and SB2 regions 
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are of the type 13 and 24. There are now also a large number of 14 type 
collisions both in the SB1 and SB2 region. There is a group of orbits near the 
island of the grey region in the SB1 and SB2 regions which fail our symmetry 
criterion. The DB region now has many collisions near the origin but still its 
mostly grey. The types of collisions we have in the DB region are either of the 
type 12 or 34. See figure (14. 3b ). 

When integrated for 1 million time steps of integration time, the grey re- 
gions of figure ( 14.3b ) in the DB region become black, except for a small island 
of grey region in the middle (figure (14.3b )). The grey areas at the beginning 
of the SB1 and SB2 regions of figure (14.3b ) survive this very long integration 
time and remain stable. 

Figures (14.41) and (14.51) show the analysis of the same orbits discussed above 



but with slightly perturbed initial conditions. It is immediately clear from 
the comparison of figures (I4.3p . (14.41) and (14. 5p that there is no significant 
difference between the results obtained using symmetric initial conditions and 
with perturbed initial conditions. 

C = 46 

Figure (14. 6p shows the results of integrations for C = 46 with symmetric initial 
conditions. Figure (14.71) and (14.81) shows the results for the same but with 
perturbed initial conditions with a perturbation of 10 -6 and 10 -5 respectively 
in the x component of Pi. 

Like the previous case the shorter version of integration shows a stable 
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double binary region surrounded by collisions of 12 and 34 types (figure HTtjb ) . 
The single binary regions (SB1 and SB2) has fewer collisions than we had in 
the previous case but a large number of orbits fail the symmetry breaking 
criterion. All the collision in the SB1 region are the 24 type. The SB2 region 
has the same characteristics because of the symmetry, but with the collisions 
being 13 type collisions. See figure (14.6b ). 

In figure (14.6b ) we show the same results but for a longer integration time 
i.e. 10 5 . In the DB region, like the previous case there are now a lot of collisions 
near the origin. All these collisions are either 12 or 34 type collision with 12 
being the most frequent. Most of the collisions in the SB1 region are 24 type. 
There are a few 12 and 14 type collisions. Most of the orbits fail the symmetry 
breaking criterion except for some near the origin of the graph. Similarly in 
the SB2 region most of the collisions are 13 type with the same number of 
12 and 14 type collisions and non-symmetric orbits as in the SB1 region. See 
figure (14.6b ). 

When integrated for 1 million time steps of integration almost all of the 
grey area in the DB region become black except for a small island of grey region 
between the collision and the non-symmetric orbits. Most of the collisions are 
12 and 34 types but there are some 13 type collisions as well. The grey areas at 
the beginning of SB1 and SB2 regions survive this long integration and remain 
stable. Most of the collisions in the SB1 area are 24 type collisions. There are 
some 14, 34 and 12 type collisions. See figure (14.6b ). 

Figure (14. 7p and H4.8[) shows the analysis of the same orbits discussed above 
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but with slightly perturbed initial conditions. It is immediately clear from the 
comparison of figures (I4.6p . (14.71) and (14.81) that there is no significant differ- 
ence between the results obtained using symmetric initial conditions and with 
perturbed initial conditions. The small numbers of single binary collisions in 
the double binary area and of double binary collisions in the single binary area 
are indications of the fact that we are heading towards hierarchical stability 
with increasing values of the Szebehely constant. 

The main characteristics of the orbits for both Cq = 40 and Co = 46 
remain the same. In the single binary regions we have comparatively more 
orbits failing the symmetry criterion in the case of Cq = 46 than we had for 
Cq = 40. For longest integration time we have slightly bigger island of stable 
orbits in the double binary region than the previous case of smaller Cq. 

C = 60 

Figure (14. 9 p shows the results of integrations for Co = 60 with symmetric initial 
conditions. Figure (I4.10p and (14.111) shows the results for the same but with 
perturbed initial conditions with a perturbation of 10~ 6 and 10~ 5 respectively 
in the x component of Pi. 

The single binary and double binary regions are completely disconnected as 
Co = 60 is larger than the critical value (Steves and Roy, 2001). In the shorter 
version of integration (figure 14.9a ) the phase space appears to be stable as 
most of the regions are grey. Both DB, SB1 and SB2 regions are surrounded 
by collisions. In the DB region collisions are either of type 12 or 34. The 
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collisions in the SB1 region are the 24 type collisions as no other type of 
collisions are possible because the phase space is disconnected. Similarly the 
only type of collisions in the SB2 area are 13 type. 

As in all other cases when integrated for a little longer (10 5 integration 
time) the grey regions in the SB1 and SB2 regions turn black except for a 
small island near the origin (figure fl~9b ). The SB1 region is surrounded by 
24 type collisions and the SB2 region is surrounded by the 13 type collisions. 
In the DB region there are large numbers of collisions around the T\ = r2 line 
near the origin. All these collisions are either 12 or 34 type collisions. When 
integrated further i.e. for 1 million time steps of integration there is no change 
in the graph which shows that the CSFBP system is more stable with larger 
Co values. 

Figures (14.101) and (14.111) show the analysis of the same orbits discussed 
above but with slightly perturbed initial conditions. It is immediately clear 
from the comparison of figures (14. 9p . (I4.10p and (14.1 ip that there is no sig- 
nificant difference between the results obtained using symmetric initial condi- 
tions and with perturbed initial conditions, which shows that perturbing the 
symmetric initial condition does not effect the final evolution by much. There- 
fore this equal mass CSFBP system is stable to small perturbations in the 
x— coordinate of the initial conditions. The absence of single binary collisions 
in the double binary area and of double binary collisions in the single binary 
area shows that this system is hierarchically stable. 

The main characteristics of the orbits for all the values of the Szebehely 
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constant i.e. C = 40, 46 and 60 remain the same. In the current case, C = 60, 
we have a much larger island of stable and symmetric orbits in the double 
binary area. 

The main conclusion from the figures is that as the value of C increases, 
the phase space of the CSFBP becomes more stable and the small perturbation 
of initial conditions does not effect the final evolution of the equal mass CSFBP 
by much. As we have stated earlier, non-symmetric orbits are not necessarily 
unstable. Therefore it will be interesting to see, in the future, the behavior 
of these orbits without the symmetry restrictions using our general four body 
integrator. 

4.7 Comparison with Szell et.al (2004) analy- 
sis in fi = 1 case 

We used the general four body equations to study the long-term evolution of 
the CSFBP system for different C values with symmetric and nearly sym- 
metric initial conditions. Szell, Erdi, Sandor and Steves (2004) investigated 
the CSFBP with the symmetric initial conditions and symmetric equations 
of motion, searching for the regular and chaotic regions using the fast chaos 
detection methods. Therefore it will be interesting to see if there are any simi- 
larities between the two investigations. Before giving the comparisons we first 
review their analysis for the /i — 1 case. Please note that Szell, Erdi, Sandor 
and Steves (2004) label the bodies in order 1 2 3 4 which means that m 1 = m 4 
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and m 2 = m 3 while we label the bodies in the order 1 2 4 3 which means that 
mi = m 3 and m 2 = m 4 . Because of the difference in labelling we have the 
following difference in labelling for the hierarchy states: 13 hierarchy state in 
Szell, Erdi, Sandor and Steves (2004) case is 14 in our case and their 23 hierar- 
chy state is 24 in our case. From now on we shall call the notation arising from 
1 2 3 4 ordering of Szell, Erdi, Sandor and Steves (2004), the CSFBP notation. 
All the CSFBP original research ((Roy and Steves, 1998), (Szell et al, 2002). 
etc) uses this CSFBP notation. Whenever we compare our work with this 
original CSFBP research, we shall indicate clearly which the CSFBP notation 
is being use. If no indication is given, it should be assumed that it is in our 
new notation. 

4.7.1 The stable and Chaotic behavior of the CSFBP- 
A Review 

Szell, Erdi, Sandor and Steves (2004) used fast chaos detection methods namely 
the Relative Lyapunov Indicator (RLI) and the Smaller Alignment Indices 
(SALI) to investigate the connection between the chaotic behavior of the phase 
space and the global stability given by the Szebehely Constant. 

In order to determine the relationship between the global hierarchical sta- 
bility of the CSFBP and the chaotic and regular nature of its phase space 
Szell, Erdi, Sandor and Steves (2004) studied two indicative CSFBP systems 
stipulated by jj, = 1 (a quadruple stellar system) and fj, — 0.1 (a binary star 
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system with two very massive planets or brown dwarf stars) by choosing a 
fine grid containing of the order of 150, 000 pairs of initial values of r± and 
T2 covering all regions of real motion in the r\T2 plane. The RLI and SALI 
numbers for each CSFBP orbit, using 10, 000 iterations for each case were nu- 
merically calculated. They chose a range of increasing Szebehely constants Co 
ie. Co < C cr iti (unstable); C cr i t i < Co < Ccrit2 (23 hierarchy stable-CSFBP 
notation) and Co > C cr u2 (stable). 

The chaotic and regular nature of the orbits is identified in Figures ( 14. 2 p 
and (I4.12p for both methods of fast chaos detection. 

We now give a brief review of their results for \i = 1. A review of their 
results /J, — 0.1 case will be given in section (14.9. ip . 

4.7.2 Review of Szell et.al (2004) results for the equal 

mass case, fi = 1 

They chose Co = 10,40 and 60 in this case. Note all hierarchies described in 
this section are in the CSFBP notation used by Szell et.al (2004). 

C = 10 (Figure 14.21 left). In this case C is well below the critical value of 

Ccriti — C cr it2 = 46.841. 

In this case the regular part of the double binary region (DB) is sur- 
rounded by a thick chaotic cloud. There are a large number of different 
type of collisions denoted by different colors. The single binary (SB) 
regions are well separated from the DB region by 12 collisions (collisions 

100 



between P% and P 2 ) (red) and contain only collision orbits indicating 
chaotic behavior. It can be concluded that the picture is very chaotic ex- 
cept for the central region of the DB area. The collision orbits dominate 
the phase space. 

C = 40 (Figure H~2| middle). In this case C is just below the critical value 

of C cr it\ = C cr it2- 

In this case the SB regions are comparatively more chaotic than the DB 
region. Most of the collisions are near the DB area. For large r\ values, 
the orbits are more regular, being light grey in both the RLI and SALI 
graphs. In the SB1 region (where r\ ^> f2), real physical orbits are likely 
to be found for greater T\ values and similarly in the SB2 (where r\ <C r 2 ) 
region, for greater r 2 values. 

C = 60 (Figure H~2l right). In this case C is much greater than the critical 
value C cr iti = C cr it2- 

The SB and DB regions are now disconnected. The phase space is regular 
except close to the boundaries of real motion where chaotic motion and 
collisions exist. The DB, SB1, and SB2 regions are surrounded by 12, 
23 (collisions between P 2 and P3) and 14 (collisions between Pi and P4) 
type of collisions respectively. Only one type of collision is seen in each 
disconnected region, corresponding to the one type of hierarchy state 
that is possible in that region. 
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Figure 4.2: /i — 1, E — 7. Graphs of the ri, r 2 initial conditions for C = 10, 
40, 60 left to right. The top line of three graphs show the RLI categorization 
of the different (ri,r2) orbits, while the bottom line of three graphs show the 
SALI categorization for the same Co values. The colors indicate collisions: red 
- "12" type, green - "13" type, yellow - "14" type, blue - "23" type. The light 
grey regions are regular, the dark grey regions are undetermined, the black 
regions are chaotic, and the white regions are forbidden for real motion. 
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4.7.3 A comparison of the Szell et.al (2004)results with 

our own CSFBP results 

We now compare the results of Szell et.al (2004) based on symmetric initial 
conditions and symmetric equations of motion with our results given in section 
5.6 based on symmetric initial conditions and the general four body equations 
of motion. In both cases the main features of the CSFBP remain the same. 

1. The stability of the CSFBP system increases as the value of the Szebehely 
constant increases. 

2. The regions of real motion are always surrounded by collision orbits. 

3. At the junction of the single binary and double binary regions almost all 
of the orbits are collision orbits 

From Szell et.al (2004) at Cq = 40, (Figure T4.21 middle), there is a dark 
cloud of black regions along the r\ = r2 line near the origin , which indicate 
strong chaos. From our work, in figure (14 .3b and 14.3b ) these orbits end up in 
collisions. Also in both cases the SB1 and SB2 regions contain most of the 
collision orbits. 

Similarly from Szell et.al (2004) at Co = 60, (Figure T4.2[ right), there is 
a dark cloud of black regions indicating chaos along the r\ = T2 line near the 
origin. In our work these regions (14. 9p after a long term integration ends up in 
collisions which shows instability in that region. In both Figures (14. 2[ right) 
and (14.9b ) the SB1 and SB2 regions are shown to be regular and symmetry is 
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maintained at least at first; however when integrated for longer times most of 
the SB1 and SB2 regions becomes black (non-symmetric orbits), figures (14. 9b ) 
and ( 14. 9b ). This does not mean that they are no longer regular orbits. It will 
be interesting to further investigate these regions in the future. Overall the 
double binary region in both the analysis are comparatively more stable than 
the single binary regions. We now look at the evolution of symmetric and 
nearly symmetric CSFBP for fi = 0.1. 



104 



(a) 




0.5 1 1.5 

(b) 

1 -5 mm riim!!iMIMiiliHif!i!lilii!!liiil!ii'! 
: H\'.V 

■i: ,« tllJifiilP 

Si, 
SI* 



■'" IH.i..-ii.I."'. ■-■.!!.■ « M. it :l r'.^U -i i .U R.H -.HI - ;:H 

0.5 1 1.5 

(c) 

Figure 4.3: /i = 1, Co = 40, E Q = —7 with no perturbation. Integration time a. 
10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: red -12 
type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, cyan 
-34 type, black -symmetry breaking and grey stable. 

105 



1 y 

y 



(b) 

1.5 mmm 
B 
v., 

R. : 

• : ::; : 

ill 



■'" IH.i..-ii.I."'. ■-■.!!.■ « M. it :l r'.^U -i i .U R.H -.HI - ;:H 

0.5 1 1.5 

(c) 

Figure 4.4: n = 1,C = 40, E = —7 with perturbation of 10~ 6 . Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.5: /i = 1,C = 40, E = —7 with perturbation of 1CT 5 . Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.6: fj, = 1,Cq = 46, E Q = —7 with no perturbation. Integration time a. 
10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: red -12 
type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, cyan 
-34 type, black -symmetry breaking and grey stable. 
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Figure 4.7: /i = 1,C = 46, E = —7 with perturbation of 10~ 6 . Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.8: /i = 1,C = 46, E = —7 with perturbation of 1CT 5 . Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.9: fi = 1, Co = 60, i^o = —7 with no perturbation. Integration time a. 
10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: red -12 
type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, cyan 
-34 type, black -symmetry breaking and grey stable. 
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Figure 4.10: /i = 1,C = 60, E = —7 with perturbation of 10~ 6 . Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.11: /i = 1,C = 60, -E = ~ 7 with perturbation of 10~ 5 . Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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4.8 Results: The nature of the orbits in the 
initial r\T2 space for fi = 0.1 

We now return to our investigation of the stability of the symmetric nature of 
the CSFBP when nearly symmetric, slightly perturbed initial conditions and 
the general four body equations are used. Here we study the case of /i = 0.1. 
We aim to understand which regions of the phase space maintain dynamical 
symmetry for long-time. The value of Eq was fixed to be -1.2 and Co was 
chosen to be 0.6, 0.7 and 0.9. The initial values of r\ and r 2 were varied from 
to 1.5 with a step size of 0.008. 

Recall that we use the general four body integrator and its not necessary 
to have symmetric collision. For example a 12 collision does not necessarily 
mean a 34 collision and vice versa. To have guaranteed symmetric collision 
one must use the symmetric equations of motion of the CSFBP (Szell, 2003). 
For color coding of the different types of orbits possible see table 14.11 

Unlike the previous case an interchange of r\ and r 2 does not produce the 
same physical orbit. Thus the categorization of the orbits is not symmetric 
with respect to the r\ = r 2 line in this case. The graphs can be seen in Figure 
(I4.13P to (14.181) . Two sets of graphs are given for each C value, one with 
symmetric initial condition and one with perturbed initial conditions. 

The graphs in this case can be separated into four regions: two double 
binary regions around r% « r 2 , r 2 < r\ (DB1) and r 2 > r\ (DB2) and two 
single binary regions, r 2 <C r\ (SB1) and r\ <C r 2 (SB2). The SB2 region 
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is very small as compared to all other regions. Also the SB2 region is very 
chaotic as most of the orbits end up in collisions. 

C = 0.6 

Figure (I4.13P shows the results of integration for Cq = 0.6 with symmetric 
initial conditions. Figure (14.141) shows the same results with perturbed initial 
conditions with a perturbation of 10~ 5 in the rc-component of Pi. 

What is immediately clear from the comparison of Figures (14. 13a ). (14. 13b ) 
and (14. 13b ) is that the longer we integrate this CSFBP system the more un- 
stable orbits we unearth. The most chaotic region for Co = 0.6 is the SB2 
region where every orbit ends up in collision or fails the symmetry breaking 
criterion. The region between n = and n = 0.2 is also very chaotic as most 
of the orbits are collision orbits. The collision types in this area are 12 and 14 
type collisions. There are some collisions along the ri-axis which are mostly of 
type 24. There are also some 12, 13, 14 and 34 type collisions. The SB1 and 
DB1 regions are separated by a thick cloud of collisions. The majority of these 
collisions are 12 and 24 type collisions. There are also some 14 and 34 type 
collisions. All of the collisions in SB2 region are of type 13. The double binary 
region is surrounded by 12 and 34 type collisions. The collisions around the 
ri = r 2 line are also 12 and 34. 

When we integrate longer, Figure (4.13b) the number of collisions increases 
in all regions except the SB2 region which didn't have any stable orbits. The 
double binary region is now surrounded by a thicker cloud of collisions which 
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are all of 12 or 34 type collisions, See Figure (14.13b ) and (14.13b ). When 
integrated for 1 million time steps of integration we get some orbits near the 
meeting point of the SB2 and DB2 region which fail the symmetry breaking 
criterion, see Figure (14.13b ). The number of these orbits is very small as 
compared to the equal mass case. The most stable region in this case is the 
SB1 region. In the DB regions the number of non-collision orbits are much 
higher than the collision orbits. We can therefore conclude that this case of 
CSFBP is more stable than the equal mass case. 

Figure (I4.14p shows the analysis of the same orbits discussed above but 
with perturbed initial conditions. Most of the orbits give the same results as 
before i.e. with symmetric initial conditions except the type of collision change 
along the r\ = r 2 line in Figure (14.14b ). In the perturbed case we get some 
24 type collisions along r\ = r 2 line which were only 12 and 34 type for the 
symmetric case. 

Co = 0.7 

Figure ( 14.151) shows the results of integration for Co = 0.7 with symmetric 
initial conditions. Figure ( 14.161) shows the same results with perturbed initial 
conditions with a perturbation of 10~ 5 in the a;-component of Pi. 

The behavior of this case is similar to the previous case. Overall we have 
fewer number of collisions than we had in the previous case but the difference 
is not very big. The only noticeable difference is in the area which separates 
the SB1 and DB1 regions. The number of collisions drops by third from the 
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previous case. See Figures (14. 15j) and (I4.16P 



C = 0.9 

Figure (14. 17ft shows the results of integration for Cq = 0.9 with symmetric 
initial conditions. Figure (14.181) shows the same results with perturbed initial 
conditions with a perturbation of 10~ 5 in the x-component of Pi. 

The single binary and double binary regions are completely disconnected 
as Co = 0.9 is greater than the critical value (Steves and Roy, 2001). There 
are no collision orbits in the SB1 region for all versions of integrations. In the 
longest version of the integration there appears an island of orbits between 
7*i = 1 and r*i = 1.5 which fail the symmetry breaking criterion. In summary, 
despite a number of symmetry breakings, the SB1 region is the most stable 
in this case too. See Figure (14.171) . The double binary region is surrounded 
by collisions of 12 and 34 type. In the shorter version of integration we have 
a very few number of collisions but when integrated for longer the number 
of collisions increases on the boundaries of the double binary region and also 
along the 7*1 = r 2 line in the DB2 region. Surprisingly there is no symmetry 
breaking which may be because of the hierarchical stability we have in this 
case. The SB2 region as in all other cases is the most unstable and none of the 
orbits are stable. The collisions in the SB2 region are all 13 type collisions, as 
no other types of collisions are possible. See Figure (14.171) . 

Figure (14.181) shows the analysis of the same orbits but with perturbed ini- 
tial conditions. There is no significant difference between the results obtained 
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using symmetric initial conditions and with perturbed initial conditions. The 
only difference is the increase in the number of non-symmetric orbits. For 
example Figure ( 14. 17b ) has no such orbits and the SB1 region is completely 
grey but we have a big island of such orbits in Figure (14.18b ) which is almost 
of the same size as we have in Figure (14.17b ). Also the island of these orbits 
in the SB1 region in Figure (14.18b ) is thicker than in Figure ( 14.17b ). 

The absence of single binary collisions in the double binary area and of 
double binary collisions in single binary areas show that this system is hier- 
archically stable. The main characteristics of the orbits for all values of the 
Szebehely constant i.e. Co = 0.6, 0.7 and 0.9 remains the same. In the current 
case, Co = 0.9, we have much less number of collision than we had in any other 
case of the CSFBP discussed in this chapter. 

4.9 Comparison with Szell et.al (2004) analy- 
sis in the \i = 0.1 case 

Before comparing the results we give a brief review of the results of Szell, Erdi, 
Sandor and Steves (2004) in the — 0.1 case. 
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4.9.1 Review of Szell et.al (2004) results for the non- 
equal mass case, fi = 0.1 

In this case Co was chosen to be 0.4, 0.8 and 0.9. For fi = 0.1 C cr ui = 0.7792 
and C crit 2 = 0.8886. Note again all hierarchies described in this section are 
the CSFBP notation used by Szell et.al (2004). 

The graphs can be separated into four regions: two double binary regions 
around r 2 ~ ri, r 2 < r\ (DB1) and r 2 > T\ (DB2), and two single binary 
regions, r 2 <C r\ (SB1) and n <C r 2 (SB2). In the following each region is 
analyzed for the given values of C . 

C = 0.4 (Figure 14.121 left). In this case C is well below the critical value 

C C ritl ■ 

The SB2 region is dominated by the " 14" (yellow) type of collision orbits. 
Some chaotic (black) orbits can be seen at around r 2 ~ 0.9. 

The DB2 region is very chaotic as suggested by the presence of large 
number of collisions. There a large number of " 12" (red) and "13" (green) 
type of binary-binary collisions which form a regular pattern. Some " 14" 
(yellow) collisions can also be found between them. The DB1 region 
has very similar behavior to the fi = 1 case. It contains a massive 
regular region with a chaotic and collision region at the boundaries of 
real motion. 

The SB1 region has a very large number of collision orbits, where " 13" 
(green) and "12" (red) collision orbits are dominant. Clear big non- 
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collision regular islands (light grey) can be found between them. It can 
be concluded that the phase space is quite chaotic with some regular 
islands. 

Co = 0.8 (Figure 14.121 middle). In this case C cr ui < Co < C cr i t 2- 

The "23" hierarchy region is now totally disconnected from the rest. 
There is a dramatic change in the SB1 region. It is now predominantly 
regular. Chaotic regions can be seen only at the boundaries. The DB1 
region is still the most chaotic. Some important changes can be found in 
the DB2 region. The non-collision area has now grown to the full length 
of the DB2 region and is more regular as is indicated by the white spots 
in the RLI graph. It can be concluded that the phase space shows more 
regular behavior. 

Co = 0.9 (Figure 14.121 right). In this case C is greater than C crit 2 and all 
hierarchy regions are disconnected from each other which means complete 
hierarchical stability. The collision orbits have almost disappeared and 
most of the orbits are regular. 

Overall, there is a connection between the global hierarchical stability cri- 
terion and the chaotic behavior of the phase space, namely as Co increases and 
the system becomes more hierarchically stable, the phase space becomes more 
regular. 
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4.9.2 A comparison of the Szell et.al (2004) results with 

our CSFBP results 

We now compare the results of Szell et.al (2004) based on symmetric initial 
conditions and symmetric equations of motion with our results given in section 
5.8 based on symmetric initial conditions and the general four body equations 
of motion. The following are the important and notable similarities: 

1. The stability of the phase space of the CSFBP depends on the C value. 
The number of collision orbits decreases as the value of Cq increases. 

2. The SB1 region is the most stable as it has very few collision orbits. 

3. For Co < C cr it there is always a big island of collision orbits at the 
junction of the single binary and double binary regions. 

From Szell et.al (2004)at C = 0.9, (Figure 5.12, right) and our case for the 
integration time of 10 4 (Figure 5.17) we have almost identical results in terms 
of the collision orbits. In the SB2 region Szell, Erdi, Sandor and Steves (2004) 
have some chaotic orbits which are also non-symmetric orbits in our case. The 
SB1 regions are identical except for the existence of a few collision orbits along 
r 2 ~ in the Szell, Erdi, Sandor and Steves (2004) case, which are symmetric 
in our case. 
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Figure 4.12: jj, = 0.1, Eq = 1.2. Graphs of the n, r2 initial conditions for 
Co = 0.4, 0.8, 0.9 from left to right. The top line of graphs show the RLI 
categorizations and the bottom the SALI. The color coding is as before in 
Figure H21 



122 



1.5 



(a) 



(b) 



Figure 4.13: /i = 0.1, Co = 0.6, Eq = —1.2 with no perturbation. Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.14: \i = 0.1, Co = 0.6, Eq = —1.2 with perturbation of 10~ 5 . Integra- 
tion time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of 
orbits: red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green 
-24 type, cyan -34 type, black -symmetry breaking and grey stable. 
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(c) 

Figure 4.15: /i = 0.1, Co = 0.7, Eq = —1.2 with no perturbation. Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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(b) 



(c) 

Figure 4.16: \i = 0.1, Co = 0.7, Eq = —1.2 with perturbation of 10~ 5 . Integra- 
tion time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of 
orbits: red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green 
-24 type, cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.17: fi = 0.1, Co = 0.9, Eq = —1.2 with no perturbation. Integration 
time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of orbits: 
red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green -24 type, 
cyan -34 type, black -symmetry breaking and grey stable. 
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Figure 4.18: /i = 0.1, Co = 0.9, E = —1.2 with perturbation of 10~ 5 . Integra- 
tion time a. 10 4 time steps b. 10 5 c. 10 6 . The colors indicate categories of 
orbits: red -12 type, yellow -13 type, magenta -14 type, blue -23 type, green 
-24 type, cyan -34 type, black -symmetry breaking and grey stable. 
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4.10 Summary and Conclusions 

In this chapter we investigated the stability of the symmetric nature of the 
Caledonian Symmetric Four Body Problem (CSFBP) by using nearly symmet- 
ric, slightly perturbed, initial conditions and the general four body equations 
to see if the CSFBP system remain symmetric. 

We analyzed the phase space in detail in the \i = 1 and /i = 0.1 cases. For 
the /i values we drew graphs of the phase space of the CSFBP with different 
Co values. Each point on the graphs, Figures (14.31) to (14. 18ft . describes a dif- 
ferent set of initial conditions of the CSFBP. In each graph we denote different 
kinds of collisions with different colors. The orbits which break the symmetry 
criterion are colored black and the orbits which maintain the symmetry are 
colored grey. The main behavior of the graphs were discussed in detail. 

We found that the stability of the CSFBP orbits is dependent on the value 
of the Szebehely constant. The larger the value of the Szebehely constant the 
more stable is the CSFBP system. For the equal mass case of the CSFBP, the 
initial conditions in the r\ — r2 space were very chaotic as most of the orbits 
ended in collision orbits. The single binary regions were the most chaotic as 
almost all of the orbits ended in collision. In the /x = 0.1 case, the SB2 region 
is the most chaotic and was comparatively very small. The SB1 region is the 
most stable as overall there were few collisions. To compare both the cases, 
/j, = 1 and /i = 0.1, the /i = 0.1 case is the most stable. In the /i = 1 case 
a large number of orbits fail the symmetry breaking criterion. These orbits 



129 



are not necessarily unstable therefore it will be interesting to investigate, in 
future, the equal mass case of the CSFBP without the symmetry restrictions. 
The effect of the symmetry breaking is considerable in only one case when for 
/i = 0.1 a large island of non-symmetric orbits appear in Figure (14. 18b ) which 
occurs after 10 5 integration time steps in the SB1 region. In all other cases of 
both fi = 1 and fi = 0.1, there are not many breakings of the symmetry and 
therefore it can be concluded that the CSFBP is stable to small perturbations. 

In both fi = 1 and /i = 0.1 case, the results were compared with Szell, Erdi, 
Sandor and Steves (2004). It is found that the main features of the CSFBP 
remain the same in both the analysis. 

1. The stability of the CSFBP system increases as the value of the Szebehely 
constant increases. 

2. The regions of real motion are always surrounded by collision orbits. 

3. At the junction of the single binary and double binary regions almost all 
of the orbits are collision orbits 

4. The SB1 region is the most stable in /i = 0.1 it has very few 
collision orbits, while in \x = 1 case the double binary regions are the 
most stable orbits. 

5. in all the comparisons made, the orbits which maintained the dynamical 
symmetry over long times were recognized as regular orbits in the Szell 
et.al (2004) investigations 
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In the next chapter we introduce a stationary mass to the centre of mass 
of the CSFBP system to derive an analytical stability criterion for the new- 
system and also to see the effect of the central stationary mass on the stability 
of the CSFBP. Later in chapter 6 we will numerically investigate its hierarchical 
stability 
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Chapter 5 

Caledonian Symmetric Five 
Body Problem (CS5BP) 

Steves and Roy (1998, 2000, 2001) have recently developed a symmetrically 
restricted four body problem called the Caledonian Symmetric Four Body 
Problem (CSFBP), for which they derive an analytical stability criterion valid 
for all time. They show that the hierarchical stability of the CSFBP depends 
solely on a parameter they call the Szebehely Constant, Co, which is a function 
of the total energy and angular momentum of the system. This stability cri- 
terion has been verified numerically by Szell, Steves and erdi (2003a, 2003b), 
while the relationship between the chaotic behavior of the phase space of the 
CSFBP and its global stability as given by the Szebehely constant is analyzed 
by Szell, erdi, Sandor and Steves (2003). 

Our aim in this chapter is to introduce a stationary mass to the centre of 
mass of the CSFBP, to derive analytical stability criterion for this five body 
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system and to use it to discover the effect on the stability of the whole system 
by adding a central body. To do so we define a five body system in a similar 
fashion to the CSFBP which we call the Caledonian Symmetric Five body 
Problem (CS5BP). 

Our motivation for studying this restricted or symmetric five body prob- 
lem comes from the three and four body problems, where restriction methods 
utilizing assumptions of neglecting the masses of some of the bodies or assump- 
tions which require specific conditions of symmetry have been very successful 
in reducing the dimensions of the phase space to manageable levels while still 
producing results which are meaningful to real physical systems. For example a 
four body model requiring symmetrical restrictions was used by Mikola, Saari- 
nen and Valtonen (1984) as a means of understanding multi-star formation in 
which symmetries produced in the initial formation of the star system were 
maintained under the evolution of the system. It is hoped that a restricted 
five body model of this type which can easily be studied, may shed light on 
the general five and four body problems in the same way that the restricted 
three body problem has helped to deepen our understanding of the general 
three body problem. 

The Caledonian Symmetric Five Body Problem has direct applications in 
dynamical systems where a very large mass exists at the centre of mass with 
four smaller masses moving in dynamical symmetry about it. The four small 
masses are affected by the central mass but are small enough that they do not 
dynamically affect the central body. Hence the central body remains stationary 
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and dynamical symmetry is maintained. This could occur, for example, in 
exoplanetary systems of a star with four planets or a planetary system with 
four satellites. 

For completion we analyze the full range of mass ratios of central body 
to the other four bodies: from the large central body system described above 
to a five body system of equal masses to a small central body with larger 
bodies surrounding it. In the case of five equal masses or a small central 
body surrounded by larger masses, the central body is unlikely to remain 
stationary as is required by the CS5BP. The model may still be applicable 
to real systems in which the outer bodies are well spaced and stationed far 
away from the central body so that they have minimal effect on the central 
body. The CS5BP model with /i — i.e. the special case of four bodies called 
the CSFBP, however, has no central body required to remain stationary. It is 
therefore a realistic model of symmetric four body systems: either four equal 
masses stellar systems or two binary stars and two planets systems. 

In section 5.1 we define the CS5BP in such a way that the CSFBP be- 
comes a special case of the CS5BP. The equations of motion and Sundman's 
inequality, the key to the derivation of a stability criterion, are given in section 
5.2. We then derive in sections 5.3 to 5.5, the analytical stability criterion for 
the CS5BP. Here we show that the hierarchical stability depends solely on the 
Szebehely constant which is a function of the total energy and angular mo- 
mentum. In section 5.5 we explain the role played by the Szebehely constant 
C in determining the topological stability of the phase space. The topological 
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Figure 5.1: The initial configuration of the CS5BP 



stability of the phase space for a wide range of mass ratios is discussed in 
section 5.6. The conclusions are given in the final section 5.7. 

This work constitutes a team effort of Steves, Shoaib and Roy. My contri- 
bution to this work was to adopt the four body model of the CSFBP by adding 
a central mass to it, derive the subsequent stability criterion, and explore the 
stability of the CS5BP for wide range of mass ratios. I also wrote the first 
draft of the paper in the accompanying material to this thesis. I improved this 
paper according to the suggestions of Steves and Roy and it has subsequently 
formed this chapter. 



135 



5.1 Definition of the Caledonian Symmetric 
Five Body Problem (CS5BP) 

The formulation of the CS5BP follows very closely from the formulation of the 
CSFBP. The CSFBP is defined by Szell et.al (2004) in the following manner. 

"The main feature of the model is its use of two types of symmetries. 1. past- 
future symmetry and 2. dynamical symmetry. Past future symmetry 
exists in an n-body system when the dynamical evolution of the system 
after t = is a mirror image of the dynamical evolution of the system 
before t — 0. It occurs whenever the system passes through a mirror 
configuration, i.e. a configuration in which the velocity vectors of all the 
bodies are perpendicular to all the position vectors from the centre of 
mass of the system (Roy and Ovenden, 1955). Dynamical symmetry 
exists when the dynamical evolution of two bodies on one side of the 
centre of mass of the system is paralleled by the dynamical evolution of 
the two bodies on the other side of the centre of mass of the system. The 
resulting configuration is always a parallelogram, but of varying length, 
width and orientation." 

The above is a description of the CSFBP. However, it can be easily seen that 
all symmetries are maintained if a fifth body of any mass is placed at the centre 
of mass of the system and required to be stationary there for all time. This 
constitutes the CS5BP. 
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Let us consider five bodies P ,Pi,P 2 ,P 3 ,P4 of masses m ,mi, m 2 ,m 3 ,m 4: 
respectively existing in three dimensional Euclidean space. The radius and 
velocity vectors of the bodies with respect to the centre of mass of the five 
body system are given by and respectively, i = 0,1,2,3,4. Let the 
centre of mass of the system be denoted by O. The CS5BP has the following 
conditions: 

1. All five bodies are finite point masses with: 

m 1 = m 3 , m 2 = m 4 (5.1) 

2. P is stationary at O, the centre of mass of the system. P 1 and P 3 are 
moving symmetrically to each other with respect to the centre of mass of 
the system. Likewise P 2 and P 4 are moving symmetrically to each other. 
Thus 

ri = -r 3 , r 2 = -r 4 , r = 0, 
Vi = f! = -r 3 , V 2 = r 2 = -r 4 , V = r = 0, (5.2) 

This dynamical symmetry is maintained for all time t. 

3. At time t — the bodies are collinear with their velocity vectors perpen- 
dicular to their line of position. This ensures past-future symmetry and 
is described by: 

n x r 2 = 0, ri • f 1 = 0, r 2 • r 2 = (5.3) 
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Figure (15.1 1) gives the initial configuration of the CS5BP. 

It is useful to define the masses as ratios to the total mass. Let the total 
mass M of the system be 

M = 2 (mi + m 2 ) + m (5.4) 

Let /ij be the mass ratios defined as ^ = ^ for i — 0, 1, 2, 3, 4. Equation (15. 4p 
then becomes 

2 (m + M2) + Mo = 1 (5.5) 

and 

< Mo < 1, 0<Mi<0.5, 0<^ 2 <0.5 (5.6) 

We simplify the problem yet further by studying solely the coplanar CS5BP, where 
the radius and velocity vectors are coplanar. Figure (15.21) gives the dynamical 
configuration of the coplanar CS5BP at some time t. In the next sections we 
will derive an analytical stability criterion for the CS5BP in such a way that 
the CSFBP becomes a special case of the CS5BP. This criterion will be applied 
to the CS5BP for a wide range of mass ratios. 

5.2 The Equations of Motion and Sundman's 
Inequality for the CS5BP 

Taking the centre of mass of the system to be at rest and located at the origin, 
the equations of motion of the general five body system may be written as 

m i r i = V i U, 1 = 0,1,2,3,4 (5.7) 
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Figure 5.2: The configuration of the coplanar CS5BP for t > 
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where Vj = i^- + j^- + k^-. i,j,k are unit vectors, along the rectangular 
axes O x ,O y , O z respectively, x^y^ Zi are the rectangular coordinates of body 
Pi and O is the centre of mass of the system. 

Given (15. 2p . the differential equations for the CS5BP reduce to 

mi ^ = ViU, 2 = 0,1,2 (5.8) 

For the CS5BP the force function U can be written as 



U = G 



1 /m? ml\ ( 1 1 

- — + — + 2mim 2 — 



2 V r i r 2 J ' \ri2 a/2 (rf + r§) - r\ 2 



, mi m 2 
+2m — + — 
T\ r 2 



(5.9) 



where ri2 = ri — r 2 . Since mo is stationary at the centre of mass, it does not 
contribute to the kinetic energy T, the angular momentum c or the moment 
of inertia I. Equations for these quantities are given below 

T = m\f\ + rri2fl, (5.10) 

c = 2 (miri x r'i + m 2 r2 x r 2 ) , (5-H) 
/ = 2 (mirf + m 2 rj) . (5.12) 
Sundman's inequality (Roy, 2004) may be written as 

u + E >Yf (5 ' 13) 

For a given energy and angular momentum, Sundman's inequality gives bound- 
ary surfaces in space which define regions of real motion. For the CS5BP, 
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these boundary surfaces can be written solely in terms of r!r 2 r 12 space as 
follows. Let Eq = —E. Then Sundman's inequality for the CS5BP becomes 

1 1 



G 



1/2 2 

1 / vn \ m 2 



2 \ r l r 2 



ru a/2 {r{ + rf) -r 



12 



-2mo ['- + -) 
n r 2 J 



> 



+ E 



a / 2 2 \ 1 "^0 ' 

4(mir^ + m 2 r 2 J 
Introducing the mass ratios /ij, ( 15. 14ft becomes 



(5.14) 



GM 



2 \ r l r 2/ v ri 2 v 2 ( r i+ r 2) - 



12 



ri r 2 



> 



4M (/iirf + /x 2 r| 



(5.15) 



Let us now introduce dimensionless variables p\, p 2 and pi 2 , and a new dimen- 
sionless constant Co, called the Szebehely constant (Steves and Roy, 2001) as 
follows 



En 



Pi -- 

Pl2 



GM 2 

Eq 



n, 



P2 



En 



GM 2 



ti 2 C 



c 2 E n 



GM 2 u G 2 M 5 ' 

where E ^ 0. Then Sundman's inequality takes the following form 

1 fu\ f4\ I 1 1 



(5.16) 



2 VPi Pa 



P12 a/2 {pi + p 2 ) - p\ 2 



.Pi P2/ 4{nip( + fjQPi) 



+ 1. 



(5.17) 



Additionally we have the following kinematic constraints on the problem 



\n - r 2 \ < r 12 <r 1 + r 2 



(5.18) 



which in the dimensionless variables becomes 



I Pi - P2| < Pl2 < Pi + P2 



(5.19) 
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Recall that given //} and /x , p-2 can be always determined using (15 .5p . Thus for 
a CS5BP system with given mass ratios of fii and /x and Szebehely constant 
(i.e. angular momentum and energy combination) Co, the equalities of ( 15. 17ft 
and (I5.19P define a surface in dimensionless coordinate space P1P2P12 which 
confine the possible motions. 

If any region of the possible space p\p2Pn is totally disconnected from any 
other, then the hierarchical arrangement of the system in that region can- 
not physically evolve into the hierarchical arrangements possible in the other 
regions of real motion. Thus a CS5BP system existing in that hierarchical 
arrangement would be hierarchically stable for all time. The topology or dis- 
connectedness of the boundary surface given by ( 15.171) and ( 15. 19ft can therefore 
provide a stability criterion for the system. 



5.3 Regions of motion in the CS5BP 

In this section, we construct explicit formulae for the boundary surface of real 
motion, enabling us to draw it in section 4, and in section 5 to identify the 
critical points for which the topology and therefore the stability of the system 
changes. 

It is useful to parameterize the surface in terms of variables (i = 1,2) 
and x%2, 

Hi = — where p n = max(p l5 p 2 ) an d £12 = ^— ■ (5.20) 

Pn Pn 

This allows us to study two halves of the phase space separately, depending 

142 



on the relative magnitudes of p\ and p 2 . 



1. if pi > p 2 , then p n = pi, y x = 1; y 2 = g. and z u = 



2. if p 2 > Pi, then p n = p 2 ; 2/2 = 1; 2/i = {J and x 12 = 



In the new variables Sundman's inequality takes the following form 



1 

Pn 



2 Vvi !fc/ I 112 V2(Sr+aI) 



\2/l 1/2 



> 



Co 



+ 1. 



In the new variables the kinematic constraint of ( 15.19P becomes 



(5.21; 



\Vi — 2/2 I < %12 < 2/1 + 2/2 



(5.22) 



Taking the equality sign in (15.211) . we obtain the following quadratic equation 
which defines the boundary surface between real and imaginary motion, 



p 2 n - Ap n + B = 0, 



(5.23) 



where 



A(yi,y 2 ,x 12 ) 



1 ^ 9 ( 1 

- — + — + 2pip 2 — + 

2 V 2/i 2/2/ 



X12 v/2 (y\ + y\) - ^ 



and 



V2/i 2/2 



4 (pxj/f + ^ 22/ 2) 

The solution of the above quadratic equation is 



1 C(y 1 ,y 2 ,x 12 ) 

2 V Pivl + P22/2 



1± Wl 



Co 



C(y h y 2 x 12 ) 



(5.24) 



(5.25) 



(5.26) 
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where 



C(yi, 2/2,^12) = A 2 fay* + fj, 2 yl) 



(5.27) 



Thus for case 1. p x > p 2 



Pi (2/2,^12) = \\\ Ci<V ^ Xl2 \ [ 1± 4/1- ^ C ° > 
2 V Pl + PbVz V V ^(2/2,^12) 



(5.28) 



where 



C(y2, X12) = + // 2 yl) 



-, 2 



, 2 ^(^ + ^ (l 4i)-J +2 ^' + S), 



(5.29) 



with the constraints < 2/2 < 1 and from (I5.19P 



1 - 2/2 < X12 < 1 + 2/2 



(5.30) 



For a given pi, the values of p 2 and pi 2 can be reconstructed by 



P2 = 2/2P1 P12 = ZiaPi- 



(5.31) 



For case (2) p 2 > p\ 



p2(yi,x 12 ) = - 



1 /C(2/i,£i 2 ) 



2 V A*i2/i + 92 



1± Wl 



Co 



C(yi,aJi2) 



(5.32) 



C(2/i,a;i2) = (/ii2/i + P2) 



-, 2 



. 2 ^(i + 7^ifHs) +2 ''»te + ^). 

(5.33) 



with the constraints < y\ < 1 and from (I5.19P 



1 - 2/1 < x 12 < 1 + 2/1 



(5.34) 
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For a given p 2 , the value of p\ and pi 2 can be reconstructed by 

Pi = 2/1)02, P12 = a;i2P2- (5.35) 

The cases 1 and 2 , given above, provide an explicit set of formulae for de- 
termining points (pi, p2, P12) on the boundary surface. The parameters (yi, x± 2 ) 
or (y 2 , X12) are simply varied from to 1. Here y-y and y 2 are the gradients of 
straight lines through the origin O in the p\Op 2 plane. x\ 2 is the gradient of 
a straight line through the origin O in either the p\Op\ 2 or the p 2 0p\ 2 plane. 

Steves and Roy (2000,2001) showed that for the CSFBP, real motion takes 
place in four tube-like regions of the P1P2P12 space that are connected to each 
other near the origin for Cq = 0. Each tube represents a particular hierarchical 
arrangement. This is also the case in the CS5BP. See figure (15.31) . 

The four type of hierarchies possible in the CS5BP are as follows: 

1. '12' type hierarchy, figure (15 .3b ). A double binary hierarchy where P\ 
and P 2 orbit their centre of mass C\ 2 . P3 and P4 orbit their centre of 
mass C34. The centre of masses C\ 2 and C 34 orbit each other about the 
centre of mass of the five-body system C. 

2. '14' type hierarchy, figure ( 15.3b ). A double binary hierarchy where Pi 
and P4 orbit their centre of mass C14. P 2 and P3 orbit their centre of 
mass C23. The centre of masses C14 and C23 orbit each other about the 
centre of mass of the five-body system C. 

3. '13' type hierarchy, figure (15.3b ). A single binary hierarchy where Pi and 
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Figure 5.3: The four possible hierarchies in the CS5BP 
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P3 orbit their centre of mass C in a small central binary. The P 2 and P4 
orbit around the central binary. 

4. '24' type hierarchy, figure (15. 3d ). A single binary hierarchy where P 2 and 
P4 orbit their centre of mass C in a small central binary. The P\ and P3 
orbit around the central binary. 

5.4 Projections in the p\-p2 plane of real mo- 
tion in the p\pip\i space 

In the CSFBP, Steves and Roy (2001) found that as C is increased, forbidden 
regions near the origin grow to the point where they meet the boundary walls 
of the tubes, resulting in disconnected regions. They found that it was possible 
to study the connectivity of the regions of motion by projecting the intersection 
of the boundary surface with the extreme values of p\ 2 on to the p\Op 2 plane. 
This is also the case for the CS5BP. 

5.4.1 Maximum extension of the real motion projected 
in pip2 space 

Equation (I5.19P and its equivalent in the new variables (6.22) give the extreme 
values of pi 2 and x± 2 respectively. Intersection of the kinematic constraints 
with the boundary surface (6.17) or (6.26) give curves projected on the p\p 2 
plane which show the maximum widths of the four tubes as three arms. The 
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two tubes located near pi~p2 lay one on top of the other in the projection, 
giving one arm near pi~p2- 

The projection curves of the maximum widths of the tubes can be found 
by substituting the equality of (6.22) into the boundary surface (6.26). Both 
limits Xi2 + =1/1 + 1/2 an d Xi 2 _ = \yi — 2/2 1 give the same equations, indicating 
that the maximum widths at the Xu + upper location and the lower location 
x\2_ are identical. 

The equations giving the maximum projections are: 



Recall that there are only two independent mass ratios required since 2(pi + 

fJa) +Po = 1- 

For case 1, p\ > P2, (15. 36ft becomes 




(5.36) 



where 




(5.37) 




(5.38) 



where 




(5.39) 
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The corresponding variable p 2 is given by 



P2 = V2p\- 



For case 2, p 2 > Pi (6.36) becomes 



^M = y-%Ml±Jl-7^T] I- I") 



where 



Ce(s/i) = (Pi2/i 2 + P2) 



1 /Vi , ,2 



2 



Pi , \ , < 1 



+2/i h p 2 + 4pi/i 2 



n 2 



.1/1 ' r "V ' " r ir "i -z/1 2 

The corresponding variable pi is given by 



(5-41) 



Pi = 2/1P2 (5.42) 

Figure (15.51) shows two typical examples of such projections for a given pi,po 
and Co- 

Real motion is possible only in the white regions. For Co7^0,a forbidden 
region (grey) exists at the origin, which grows as Co is increased to meet the 
forbidden region surrounding the exterior of the arms. The p 2 << Pi arm 
represents the '24' type of hierarchy, the p 2 ~ pi arm represents both the '12' 
and '14' type of hierarchies and the pi << p 2 arm represents the '13' hierarchy. 
The connected projection in figure (5a) for C = 0.039 indicates that change 
from one hierarchical arrangement to another is possible. Figure (5b) gives the 
projection for Co = C cr u = 0.055, a critical value, at which the regions of real 
motion become disconnected. For C > C crit all allowed regions are totally 
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disconnected and evolution from one hierarchy to another is impossible. The 
system is therefore hierarchically stable. 

5.4.2 The minima of the boundary surface of real mo- 
tion projected in p\p2 space 

The minima of the boundary surface give information on the three dimensional 
shape of the surface. Projection of the curves indicating where the minima are 
located in pip 2 space are useful in identifying when, as C is increased, holes 
first appear within the boundary surface. Motion is still possible from one tube 
to another by moving around the central holes. See Steves and Roy (2001) for 
further details. 

For a given yi,y2, the minima of p n with respect to xyi occur at x\ 2 = 
\Jy\ + y\- The projection of the minima onto the p\p 2 plane are given by 



where 




(5.43) 



(5.44) 



For case 1 p\ > p 2 (6.43) becomes 

l/^Mjl ± A- «• ), (5.45) 
2 y Hi + p 2 yi \ y C m (y 2 ) J 
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where 



C m (y 2 ) = (hi + H2yf) 



2 1^1 + V2 



+2/i + 



For case 2 p 2 > Pi (6.43) becomes 



where 



Cm(yi) = (a*ii/i + P2) 



1 />? 



2 V2/1 



+ P2 



\yi J VW+i 



(5.46) 



(5.47) 



(5.48) 



5.5 The Szebehely ladder and Szebehely Con- 
stant 

Through the projections in the pip 2 plane given by the maximum extensions 
and of the minima of the boundaries of real motion in the p\p 2 p\2 space, we 
can study the topology of the boundary surfaces and thus gain knowledge on 
the hierarchical stability of the system. The topology changes as Co increases. 
The critical values of C at which the space becomes disconnected therefore 
provide stability criterion. 

The value of p n (yi,y2), for the maximum extensions and the minima pro- 
jections, explicitly depends on the value of C(yi,y 2 ) i.e. it has two roots, a 
single double real root or imaginary roots, if C is greater than, equal to or less 
than C respectively. 
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The Equations C e (y 2 ), C' e (yi) ,C m (y 2 ) , C' m (yi) therefore give information 
on the point at which the topology of the projections changes. The critical 
changes occur at C(y\, y 2 ) = C , the single real root solution. For example C e 
can be evaluated for the range of y 2 from to 1. Recall that y 2 is the gradient 
of a straight line through the origin O in the p±Op 2 plane. 

The minimum value of C e (y 2 ), C e (min), is the first value of Co, as it is 
increased, where there is only one solution (pi, p 2 ) to the maximum projection 
curve. For Co > C e (min), there are no solutions (pi,p 2 ) at y 2 = min and the 
projection becomes disconnected indicating a stable system. 

The minima of the four C-functions indicating the point of change in the 
topology can be thought of as the rungs of a ladder, which Steves and Roy 
(2001) call the Szebehely ladder. The rungs of the ladder Ri = C m (min); 
R 2 = C^(min); R 3 = C e (min) and _R 4 = Cg(min) are dependent only on the 
masses of the system and are invariant to the initial conditions or the c and E 
of the system. 

Both y\ and y 2 lie in the range to 1. Hence we can combine relations 
C e (y 2 ), Cg(yi), C m (y 2 ), C' m (yi) in the same figure plotting C against y. See 
figure (5.4), where the minima of the curves of the four C— functions form the 
four rungs of the Szebehely ladder. The stability of the system depends on 
the location of the Szebehely constant C for the system with respect to these 
rungs. Thus for 

1. Co > Ri, there is a hole of forbidden motion near the central region 
connecting the four tubes within the boundary surface for pi > p 2 . 
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2. Cq > i?2, there is a hole of forbidden motion near the central region 
connecting the four tubes within the boundary surface for p 2 > p±. 

3. Cq > i?3, the arms in the projection of the maximum extensions for 
Pi > P2 are disconnected and the '24' hierarchy is stable. 

4. Cq > i?4, the arms in the projection of the maximum extensions for 
P2 > Pi are disconnected and the '13' hierarchy is stable. 

When Cq > max(i?3, R4), all arms are disconnected and the system is hierar- 
chically stable. 

1. If H2 > Hi, then i? 4 > i? 3 > i? 2 > R\. 

2. If /ii > /i 2 , then R 3 > i? 4 > i?i > R 2 . 

3. If Hi = fi2, then R 2 = Ri and -R3 = R±. 

Thus the critical value of Cq at which the whole system becomes hierarchically 
stable for all time is 



Hi = H2 is the special case of equal masses where Cq > (C cr i t = i?3=i?4) 
gives total hierarchical stability at one critical point. Otherwise, hierarchical 
stability occurs in two stages Co > (C cr i tl = R3) and Co > (C cr i t2 = -R4). If 
/io = 0, then we have the special case of the CSFBP discussed by Steves and 
Roy (2001). The critical value of Co at which the system becomes hierarchically 
stable for all time is given by (15.491) . R3 and R± are purely functions of /xq 




(5.49) 
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and fj,±. For Hi = fi 2 , they are functions only of fi . Figure (6.8) plots these 
critical values as a function of /i . For C > C crit (fi ), hierarchical stability 
is guaranteed. Figure (6.8) shows that C cr u(fio) has a maximum of 0.065667 
at fi = 0.183. Thus if C > 0.065667, all CS5BP's with ^ = // 2 will be 
hierarchically stable. 

We now show several examples for a range of mass ratios of how rungs of 
the Szebehely ladder can be computed solely from /ii and /i . Then using the 
value of the Szebehely constant Co for the system, which depends on the initial 
conditions, the hierarchical stability of the system can be determined. 



5.6 The Stability of the CS5BP systems with 
a range of different mass ratios 

5.6.1 The Equal Mass CS5BP 

The equal mass CS5BP has /ii = fi 2 = Ho = 1/5. In this case there exists only 
two rungs of the Szebehely ladder; since C m = C' m and C e = C' e , viz. 



where < y < 1. 

The minimum values of C m (y) and C e (y) form the two rungs of the Szebe- 
hely Ladder. The minimum value of C m is 0.039222 = R\ and the minimum 




(5.50) 



(5.51) 
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Figure 5.5: — /i 2 — Ho'- The projection of the boundary surface onto the 
pi - p 2 plane a. C = Ri = 0.0392219 b. C = R 4 = 0.0655514 

value of C e is 0.065551 = _R 4 . R\ and i? 4 occur at y = 1 and y = 0.472 
respectively. 

Figure (15 .5b ) shows the projection of the maximum extensions for Cq — R\. 
The phase space remains connected but a small forbidden region exists near the 
origin. This forbidden region grows as Cq is increased until at Cq = R4 at the 
highest rung of the ladder, the phase space becomes disconnected. See figure 
(15. 5b ). Note that the gradient y of the line OA (passing through the point 
of single solution) is 0.472. The five body equal mass CS5BP is hierarchically 
stable for values of C greater than i? 4 = 0.06555. 

5.6.2 Four equal masses with a varying central mass /io 

In this case there exists only one independent mass ratio since given /io, 
/ii = l/4(l-/x ) from (6.5) and fi 2 = Pi- Since fii = fi 2 , only two rungs of 
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Figure 5.6: The projections of the boundary surfaces onto the p\ — p 2 plane 
for Co = 0, H\ = fi2 and a range of /iq from to 0.96 
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the Szebehely ladder exist i.e. C m = C' m and C e = C' e 




(5.52) 



C e = p 1 (l + y 2 ) + (1 + -)+^^) • (5.53) 

Figures (15. 7p and (15.91) show the projections of the maximum extensions onto 
the p\Pi plane for two typical cases: 

1. a small central mass; p\ = p<i = 22 1 qq 5 , Po = 0.01 and 

2. a large central mass; /ii = p 2 = 0.01, /xo = 0.96. In each figure, two 
values of Co, Co = Rx, and Co = -R4, have been selected. 

For CS5BP's with a small central mass, the largest region of real motion is 
found to be the central arms where p\ ~ p^. This indicates that such systems 
are most likely to be moving in double binary hierarchies of type '12' and '14', 
figure flSTTj) . 

CS5BP's with large central bodies, however, have their largest region of 
real motion in the arms p\ ~ and p2 ~ 0, indicating that the single binary 
hierarchies '13' and '24' will be the dominant systems. 

It is interesting to study the effect of placing a small mass at the centre of 
mass and increasing its mass from to 1. Figure (15. 6p shows the projections 
for C = and for a range of p . 

For /i = to 0.2, i.e a small central mass, the double binary hierarchies 
dominate, with single binary hierarchies more prevalent as po increases. At 
p =0.2, i.e. the five body equal mass case, the areas of real motion are of 
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Figure 5.7: \i\ — /j, 2 — \qI , Ho = 0.01: The projection of the boundary surface 
onto the Pl - p 2 plane at a. C = R x = 0.0295707 b. C = Ra = 0.048036 



0.07 




^0 

Figure 5.8: The critical values of C , C crit , at which the CS5BP becomes 
hierarchically stable as a function of /x , where /ii = /i 2 
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Figure 5.9: — [12 — 0.01, /i = 0.96: The projection of the boundary surface 
onto the Pl - p 2 plane at a. C = Ri = 0.0000301 b. C = R 4 = 0.0000323 

relatively equal sizes for the double binary and single binary hierarchies, sug- 
gesting neither is dominant. When comparing the area of real motion available 
for /iq — (the four body equal mass case), with that of /1q = 0.2 (the five 
body equal mass case), we see that the addition of a fifth body of equal mass 
at the centre increases the area of real motion in both single binary and double 
binary hierarchies. It thus most likely increases the chances of moving from 
one hierarchy to another. Thus it is likely to be more hierarchically unstable; 
as would be expected. 

For hq = 0.2 to 1, i.e a larger central mass, single binary hierarchies dom- 
inate, with double binary hierarchies becoming virtually non-existent for /i 
close to 1. no close to 1 represents a star with four planets or a planet with four 
satellites. In such situations, it is highly unlikely that the four small bodies 
will form two binary pairs orbiting the central body. 
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Figure 5.10: Hi = 0.195, /i 2 = 0.3 and Ho = 0.01. The projection of the 
boundary surface onto the p\ — p 2 plane at a. Co = R\ = 0.0281 b. 
C = R A = 0.0270 

5.6.3 Non-equal masses i.e. \i\ ^ fi2 ^ f^o 

With in 7^ fi2 7^ Ho we now have two independent mass ratios Ho-, Hi, since 
Hi = 1/2(1 _ Ho) — Hi, by (6.5). We also have four rungs of the Szebehely 
ladder i.e. C m ^ C' m and C e ^ C' e 

Figure (15.111) and (I5.13P give two typical examples of projections for 1) 
< Hi < Hi\ and 2) /i 2 < Ho < Hi- I n eacn figure, two values of C , 
Co = i?3 and Co = R4 have been selected to show the two stages of increasing 
hierarchical stability. For example, in figure (15. 11 j) Hi < Hi, therefore when 
R3 < Co < Ra, the arm p 2 ~ becomes disconnected first and any system in 
a '24' hierarchy will be stable. If the system is in any other hierarchy it is still 
free to change to all hierarchies, but the '24' hierarchy. See Figure (15.11a ). 
Once Co > R4, all arms become disconnected and the system is hierarchically 
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Figure 5.11: [i\ = 0.195, /i 2 = 0.3 and /i = 0.01. The projection of the 
boundary surface onto the pi — p 2 plane at a. C = R 3 = 0.0439 b. 
C = R A = 0.0470 
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Figure 5.13: p,\ = 0.3, jj 2 = 0.1 and jj = 0.2. The projection of the boundary 
surface onto the p\ — p 2 plane at a. Co = R3 = 0.051 b. Co = -R4 = 0.0553 

stable for all hierarchical arrangements. See figure (15.11b ). 

Note that in figure (15.131) . where p\ > p 2 the arm p\ w becomes discon- 
nected first. Thus for R 3 < Co < R4, any system in a '13' hierarchy will be 
stable. 

The critical value of Co at which the whole system becomes stable is given 
by (I5.49P and is only a function of /i and p\. 

Figure (5.14) plots these critical values as a function of po and p\. For 
Co > C cr it(pjo, fii), hierarchical stability is guaranteed. Figure (5.15) shows a 
cross-section of Figure (5.14) for the value /j,q = 0.2, describing the two curves, 
the critical values -R3 and R4 as functions of p\. Figures (5.14) and (5.15) 
show that C cr it(po, pi) has a maximum of 0.0659 at (/Jo,£*i) = (0.2,0.223) and 
(^o,/xi) = (0.2,0.185). Thus if C > 0.0659, all CS5BP's regardless of their 
mass ratios, will be hierarchically stable. 
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Figure 5.14: The critical values of C ,C crit , at which the CS5BP becomes 
hierarchically stable as a function of /io, l^i 
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Figure 5.15: The critical values of C , C crit , at which the CS5BP becomes 
hierarchically stable as a function of where /x = 0.2 
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.7 Difference of Notation with Steves and Roy 
(2000, 2001) explained 

1. Steves and Roy (1998, 2000) began numbering their CSFBP in numerical 
order (1234), thus the symmetric pairs were 1) Pi and P4 and 2) P2 and 
P3. 

In this current work it was realized that when the four body problem 
was generalized to higher number of bodies, mathematically it would 
be advantageous to number the bodies so that if the first body in the 
pair was the }th, the second body in the pair would be numbered n + j, 
where 2n is the total number of bodies in the system. For the four and 
five body problems this meant the symmetric pairs became: 1) P\ and 
P3 and 2) P 2 and P4. Thus when comparing Steves , Roy and Szell's 
original CSFBP results with the results of this thesis, the following are 
equivalent labels given in table (15.11) : 

2. For the CSFBP, Steves and Roy (2001), need to define only two different 
masses m and M. They, therefore need only one mass ratio to describe 
the system. Thus in the four body symmetrical system, fi is defined as 
fi = m/M. when more than four bodies are included in the problem it is 
easier to use a more general system of mass ratios. Thus fii is chosen to 
eb the ratio of the ith body to the total mass of the system. Therefore we 
have a scaling difference between the Steves and Roy original notation 
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Table 5.1: Difference in labelling of hierarchies with Steves and Roy (1998, 
2000) 



Original Notation 


Current notation 


12 hierarchy (double binary (DB)) 


12 hierarchy (double binary (DB)) 


13 hierarchy (double binary (DB)) 


14 hierarchy (double binary (DB)) 


14 hierarchy (single binary (SB)) 


13 hierarchy (single binary (SB)) 


23 hierarchy (single binary (SB)) 


24 hierarchy (single binary (SB)) 


and our notation of 





= (2^T2) " 2 = (2^W (5 ' 54) 

where \x = m/M as defined by Steves and Roy (2001). Hence for /i — 1 
i.e. the equal mass four body problem we have \i\ = 0.25, fj, 2 = 0.25 and 
A*o = 0. 

3. For the same reasons as above, we have a scaling difference for the Sze- 
behely Constant Co for which we give the following conversion formula 

C A = (2fi + 2) C s , (5.55) 

where Ca is the Szebehely Constant given by Steves and Roy and Cs is 
the Szebehely Constant given by ourselves. 
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5.8 Conclusions 

Steves and Roy (1998, 2000, 2001) have recently developed a symmetrically 
restricted four body problem called the Caledonian Symmetric Four Body 
Problem (CSFBP), for which they derive an analytical stability criterion valid 
for all time. We introduced a stationary mass to the centre of mass of the 
CSFBP and derived analytical stability criterion for the resulting five body 
system (CS5BP). The stability criterion was then used to determine the effect 
of adding a central body on the stability of the whole system. 

The critical value of C at which the whole system becomes hierarchically 
stable for all time is 



/ii = /i2 is the special case of equal masses where Co > (C cr i t = Rz=Ra) 
gives total hierarchical stability at one critical point. Otherwise, hierarchical 
stability occurs in two stages C cr m = Rs < Co > C cr u 2 = Ra and Co > 
(Ccrit 2 — R&)- R3 an d R4 are purely functions of /j,q and jii. For /ii = /i 2 , 
they are functions only of ji^. Figure (6.8) plots these critical values as a 
function of /x - For C > C crit (fi ), hierarchical stability is guaranteed. Figure 
(6.8) shows that C cr j t (/i ) has a maximum of 0.065667 at /i = 0.183. Thus if 
Co > 0.065667, all CS5BP's with /ii = /i 2 will be hierarchically stable. 

For the cases of non-equal masses and C cr ui < Co < C cr u2, one hierarchy 
type depending on the relative size of the masses becomes disconnected and is 
therefore stable. For the cases of either /i > ji 2 or /i > Hi, the 13 hierarchy 




(5.56) 
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state is hierarchically stable while in all other cases the 24 hierarchy state is 
hierarchically stable. 

For /jl — to 0.2, i.e a small central mass surrounded by larger masses, the 
double binary hierarchies dominate, with single binary hierarchies becoming 
more prevalent as /i increases. At /i = 0.2, i.e. the five body equal mass case, 
the areas of real motion are of relatively equal sizes for the double binary and 
single binary hierarchies, suggesting neither is dominant. When comparing the 
area of real motion available for /i — (the four body equal mass case), with 
that of fj, = 0.2 (the five body equal mass case), we see that the addition of 
a fifth body of equal mass at the centre increases the area of real motion in 
both single binary and double binary hierarchies. 

In the CS5BP, Figures (5.14) and (5.15) show that C cr i t (no, ^i) has a max- 
imum of 0.0659 at (yu , A*i) = (0.2,0.223) and (n ,Hi) = (0.2,0.185). Thus if 
Co > 0.0659, all CS5BP's regardless of their mass ratios, will be hierarchically 
stable. 

We also show in section 6.4 that when /x > /ii or /i > /i2, the 13 hierarchy 
state is hierarchically stable for R% < Co < Ra while in all other cases the 24 
hierarchy state is hierarchically stable for R 3 < C < Ra- 

The analytical stability criterion derived here tells us about the complete 
hierarchical stability at the critical value of C . It does not tell much about 
what happens in between C = and C crit . To find out this answer in chapter 
6 we will numerically integrate many CS5BP systems by using the analytical 
stability criterion as a guide. We will then analyze the stability of the differ- 
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ent hierarchy states by studying the frequency of changes from each specified 
hierarchy to another specified hierarchy. 
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Chapter 6 

Numerical investigation of 
Hierarchical Stability of the 
Caledonian Symmetric Five 
Body Problem (CS5BP) 

In the last chapter we derived an analytical criterion for the topological sta- 
bility of the CS5BP. In this chapter we investigate the hierarchical stability of 
the CS5BP numerically, and compare the results with that of the analytical 
stability criterion. We verify that for C > R 4 , the CS5BP system is hierarchi- 
cally stable and for R 3 < Co < R4 the CS5BP system is partially stable. It is 
also shown that with increasing value of C , the system becomes more stable. 

In section 6.1 we provide a brief review of the numerical investigation of the 
hierarchical stability of the Caledonian Symmetric Four-Body Problem (Szell, 
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2003). In section 6.2 we determine the equations of motion of the coplanar case 
of the CS5BP and describe our criterion for detecting the hierarchy changes. 
In section 6.3 we briefly describe the integrator software designed by ourselves 
to solve the differential equations given in section 6.2. The equations for the 
initial conditions are derived in section 6.4. In section 6.5 we give a detailed 
analysis of the hierarchical stability of the CS5BP in which we also discuss the 
CSFBP as a special case of the CS5BP. Here we are interested in understanding 
the relationship between the Szebehely constant Co and hierarchical stability. 
Section 6.6 contains the conclusions. 

6.1 Review of the Hierarchical Stability of the 
Caledonian Symmetric Four Body Prob- 
lem (CSFBP) 

Let a be the angle between P\ and P2 in figure (I6.ip . With the help of a we can 
determine the hierarchical evolution of the system (Szell, 2003). The CSFBP 
has four different hierarchy states which are defined by Szell as follows. Please 
note that the definition of hierarchy states for the CS5BP is slightly different 
from those of the CSFBP, see section 5.7 for details. 

1. 12 type hierarchy. A double binary (DB) hierarchy where Pi and P2 
orbit their centre of mass C\2- P3 and P4 orbit their centre of mass C34. 
The centre of masses Cyi and C34 orbit each other about the centre of 
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mass of the four-body system C. 

2. 13 type hierarchy (14 in the case of the CS5BP). A double binary (DB) 
hierarchy where P\ and P3 orbit their centre of mass C13. P 2 and P4 
orbit their centre of mass C 24 . The centre of masses C13 and C 24 orbit 
each other about the centre of mass of the four-body system C. 

3. 14 type hierarchy (13 in the case of the CS5BP). A single binary (SB) 
hierarchy where P\ and P 4 orbit their centre of mass C in a small central 
binary. The P 2 and P 3 orbit around the central binary. 

4. 23 type hierarchy (24 in the case of the CS5BP). A single binary (SB) 
hierarchy where P2 and P3 orbit their centre of mass C in a small central 
binary. The P 1 and P 4 orbit around the central binary. 

Szell (2003) uses the following criterion to determine the hierarchical posi- 
tion of the four body system under discussion. 

Let k be an integer. When a oscillates inside the interval (2kn — |, 2kii + 
|), then we have a 12 type hierarchy. When it oscillates inside the interval 
(2kn + 7r — |, 2kn + n + |), k e Z, interval then we have a 13 type hierarchy. 

In the case when a is a monotonically increasing function and there exists 
t\ and t 2 > h for which a(ti) = (Ak ± 1)| and a(i 2 ) = (4/c ± 3)| we have 
either a " 23" type hierarchy or " 14" type hierarchy. If r\ > r 2 it is a " 23" type 
hierarchy and if r\ < r 2 it is a " 14" type hierarchy 

Now we will briefly discuss the results of Szell's numerical integration into 
the frequency and type of hierarchy changes for different mass ratios. We will 
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Figure 6.1: The Caledonian Symmetric Four Body (CSFBP) model 

particularly use later the results for \i = 1, when comparisons are made with 
the CS5BP numerical investigations. 

6.1.1 Hierarchy changes for different mass ratios of the 
CSFBP 

In this section we give a short summary of the results of Szell (2003) for \x = 1, 
H = 0.1, fj, = 0.01 and /i = 0.001 of the CSFBP. 

H = 1 Case : There are no hierarchy changes for largest critical value of Co 
i.e. i?3 = i?4. There are numerous hierarchy changes for all smaller values of 
Co. The number of hierarchy changes reduces as Co increases. Thus it can be 
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said that the system is more unstable for small C values. The most unlikely 
changes are 12 — > 13 and 13 — > 12, which is a change from one double binary 
hierarchy to another double binary hierarchy. The most likely changes for Co 
less than Ri are 23 — > 14 and 14 — > 23 ( single binary to single binary hierarchy 
changes and vice versa). For bigger C these hierarchy changes decrease and 
the 23 — > 12 and 23 — > 13 (SB to DB and vice versa) hierarchy changes become 
dominant. Near the critical value R 2 the number of hierarchy changes decreases 
as one would expect. Only the 23 — > 12 and 23 — > 13 type of hierarchy changes 
remain considerable. 

/i = 0.1 Case : There are no hierarchy changes to and from the 23 hi- 
erarchy state near i? 4 . With small C values 14 — > 23 and 23 — > 14 type of 
hierarchy changes are dominant. As in the previous case 12 — > 13 and 13 — > 
12 are the most unlikely hierarchy changes. 

H = 0.01 Case : There are no 12 -> 23, 13 -> 23, 23 -> 13 and 23 -> 14 
hierarchy changes for C near i? 4 . With small C , the 14 — > 23 and 23 — > 14 
hierarchy changes are dominant. 

/i = 0.001 Case : There are very few hierarchy changes for this mass ratio 
which indicates that the system is very nearly hierarchically stable for all Cq 
values chosen. 

The overall behavior of the system for different mass ratios indicate that 
as the value of ji decreases the system becomes hierarchically more stable. 

We will now derive the equations of motion and initial conditions in order 
to perform a similar numerical investigation of the CS5BP problem. 
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6.2 The Equations of Motion of the CS5BP 
and Hierarchy Changing Criterion 

Note that from now we will use the CS5BP notation for numbering the five 
bodies. See section 5.7. The classical equations of motion for the general 
n-body problem are given by 

— p-^-ry, z = 0,l,2,3... (6.1) 

where = (xi, yi) and = — r^. For a five-body problem we will get the 

following equations of motion from (16. ip above 

mir 01 , m 2 r 02 m 3 r 03 , m 4 r 04 la n , 

r = — 3— + — 3— + — 3— + — 5— (6.2) 

ryt 'J / yt '.J :yt ' J ryt 'J 

'01 '02 '03 '04 

m Q r w m 2 r 12 m 3 r 13 m 4 r u 
r 1= — 3— + — 3— + — 3— + — 3— (6.3) 

'10 '12 '13 '14 

m r 20 , 77nr 2 i m 3 r 23 m 4 r 24 
r 2 = — 3— + — 3— + — 3— + — 3— 6.4 

'20 '21 '23 '24 

mor 3 o , m!r 31 m 2 r 32 m 4 r 34 
r 3 = — 3— + — 3— + — 3— + — — (6.5) 

ryt 'J v-* '.J . v-< T_) . 'J 

'30 '31 '32 '34 

m r 40 m 4 r 4 i m 2 r 42 m 3 r 43 
r 4 = — ^ 1 7. 1 = 1 5 — (6.6) 

'40 '41 '42 '43 

By utilizing all the symmetries of the coplanar CS5BP, see chapter 5 for details, 
we reduce the number of equations from fifteen to four. As m is stationary 
at the origin for all time, therefore we have ro = tq = ro = 0. Also, as 
r*i = — r 3 and r 2 = — r 4 therefore we need to solve for r 4 and r 2 or r 3 and r 4 
only. The equation of motion for the coplanar CS5BP in the simplified form 
are the following 

If. mi\ fr 1 -r 2 r 4 + r 2 \ 

p i= (™° + -) r i " ™* {— + -JT- ) ( 6 - 7 ) 
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r 2 = --3 ( m + — J r 2 - m 2 — — + — — ) (6.? 

'2 v ^ ' V '12 '23 



As ri = (xi,yi) and r 2 = (£2,2/2), equations ( j6. Tf) and ( 16.8f) can be rewritten 



as 



1 / mi\ 
x a = 5- m + — xi - 

(x 2 + 2/ 2 )^ V AJ 
f x 1 -x 2 X1 + X2 . , 1V 

m 2 2~T 2~T l"'^ 

V ((xi - a; 2 ) 2 + (yi - 2/2) ) 2 ((xi + x 2 ) 2 + {yi + 2/2) ) 2 



{xl + ylY 



m 1 \ 
m + — j y 1 - 



m 2 



2/i - 2/2 



+ 



2/i +2/2 



2 3 2 3 

((xi - x 2 ) 2 + (yi -y 2 ) ) 2 (Oi + x 2 ) 2 + (yi +y 2 ) )3 



5.10) 



•i' 2 



0&2 + 2/2, 



-3 i m + — ) x 2 



m 2 \ 
4 / 



£ 2 - £1 



+ 



£1 + £2 



((xi - x 2 ) 2 + (yi - y 2 ) 2 ) a ((xi + x 2 ) 2 + (y a + y 2 ) 2 ) I 



2/2 



1 



1^2 + 2/1 



2^1 



m + 



m 2 \ 
4 / 



2/2 



mi 



2/2-2/1 



2 3 / > 

((xi - x 2 ) 2 + (yi -y 2 ) )s ((xi + x 2 ) 2 + (j/i + j/ 2 ) ) 2 



+ 



2/i +2/2 



5.12) 



To monitor the hierarchies and its changes from one state to another we 
solve numerically equations (16. 91) to ( 16.121) at any time t to give the locations 
of the four bodies in the phase space with respect to each other. We will 
explain section 6.3 how we solve these equations. Once we find the position 
co-ordinates of all the four bodies, then it is not difficult to determine what 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 6.2: Regions of allowed real motion (white) in the p\ — P2 — P12 space 
projected onto the p\ — p 2 plane 
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hierarchy state the CS5BP system is in. See section 6.7 for the definition of 
hierarchy states for the CS5BP. To monitor the hierarchy changes during the 
integration time we make use of figure ( 16. 2(1 . Please note that figure (16.21) is 
different for each set of initial conditions. Region I and III represent single 
binary areas while region II represent double binary areas. If the system is in 
region I we say that it is in the 24 type of hierarchy state, if it is in region II 
we say that it is in the 12 or 14 hierarchy state depending on the distances 
I-P1-P2I and |PiPzi|. If I-P1-P2I < I-P1-P4I) then the system is in the 12 type of 
hierarchy state. Otherwise it is in the 14 type of hierarchy state. If the system 
is in region III then we say that it is in the 13 type of hierarchy state. But if 
the system is in region IV then we do not recognize any hierarchy state until 
it enters one of the other three regions. When the system transfers from one 
region to another, we record a hierarchy change. For example if the system 
goes out of region I and enters region IV and then region III we record this as 
a 24 — > 13 (2313) hierarchy change and vice versa. 

Note that this criterion for recognizing a hierarchy change is different from 
that of Szell (2003). Our criterion is based on the sizes of r\ and relative 
to each other while Szell's criterion is also based on changes from libration to 
rotation and vice versa of the angle a between Pi and P^. In practical terms, 
both methods produce similar frequencies of hierarchy type changes. 
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6.3 The Integrator 

The equations of motion for the CS5BP, equations (16. 9p to (I6.12p . are highly 
non-linear second order coupled differential equations. It is not possible to find 
their analytical solution. Therefore we have to use some numerical technique 
for integrating this system of the differential equation. The software must take 
as input data, the initial parameters of the differential equations and give as 
output data, the solution ( position and velocity coordinates) of differential 
equations belonging to the input data. To develop the integrator we first 
needed to find an accurate and fast numerical method and then develop an 
environment where the input and output data is efficiently handled. 

The numerical method, we chose, for this integration project is a 15th order 
method with a adaptive step size control called the Radau method of Ever- 
hart (Everhart, 1985). This method makes use of Gauss- Radau spacing. We 
have adopted the FORTRAN code from the example of Everhart's three body 
problem. We used Microsoft Visual C++ 6.0 to construct the environment of 
the integrator. 

6.4 Initial Conditions 

The numerical integrations of (16.91) to (16.121) require initial values for ri and 
r 2 . In order to satisfy the initial conditions of the CS5BP we immediately have 

rx = (xi,0) and r 2 = (x2,0) (6.13) 
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We choose a grid of values of x\ and %2 so that they range within the 
intervals (0, 3] and (0, 2.5], separated by an increment size 0.05. The initial 
velocities Vi of Pi and V 2 of P2 are calculated using the following relations. 



X\cm\ — X2A/m 1 m 2 (— c 2 + 4(x 2 m 1 + xlm 2 )(E + U)) 



(6.14) 



12/ - 



2m\{m\x\ + 1122X2) 



x 2 cm 2 + Xx ^mim 2 (-c 2 + 4(a; 2 mi + xlm 2 ){E + U)) 



(6.15) 



2y - 



2m 2 (mixf + 7712X2) 



Where V\ y and V 2y are the y components of Vi and V 2 respectively, c = 



J 3^ is the angular momentum of the system, Eq is negative of the energy E 
and U is the potential given in chapter 4. The x components of Vi and V 2 
are set to be zero for t = 0. 

Our aim in this chapter is to examine the motion and stability of a variety 
of the CS5BP systems each with a different Szebehely constant Cq. We also 
wish to compare systems with the same Szebehely constant. Thus we select 
the initial conditions so that they result in the same Cq. For a given Co, the 
variables V\ y and V 2y can be calculated from equations f)6.14p and (16.151) . 

We investigate the following sets of mass ratios of the CS5BP with several 
values of Cq for each set of mass ratios. Diagram of the initial configuration 
for each set of mass ratios are given in table (6.1), along with section in which 
each system is discussed. 




1. The Four Body Cases i.e the CSFBP 
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(a) Hi — /i 2 — 0.25 and /i = 0, in CSFBP notation fj, — 0. 
C = {0.01,0.028,0.031,0.038,0.042,0.046} 

(b) fa = 0.05, fi2 = 0.45 and /i = 0, in CSFBP notation \l = 0.1. 
C = {0.012,0.013,0.015,0.018,0.019} 

(c) /i! = 0.00495, n 2 = 0.49505 and = 0, in CSFBP notation fj, = 
0.01. 

Cq = {0.008,0.0082,0.0085,0.0086,0.0087} 

(d) ^ = 0.0004995, fi 2 = 0.4995095 and /i = 0, in CSFBP notation 
H = 0.001. 

Co = {0.0078,0.00786,0.00788,0.00789,0.0079} 

2. Equal mass case of the CS5BP i.e. /ii — /jl 2 — A*o = 0.2 
Co = {0.03,0.06,0.07} 

3. Four equal masses with a varying central mass /x 

(a) — 2 ioq 5 an( i A*o — i^o j Four equal masses with a very small 
central mass 

Co = {0.026,0.04,0.05} 

(b) fj,i — f/,2 — g and /i = |; Four equal masses with a comparatively 
larger central mass but still smaller than the outer bodies 

Cq = {0.036,0.055,0.063} 

(c) /ii = /i 2 = an d A*o = ^jj Four equal masses with a large central 
mass and smaller outer bodies. 
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C = {0.00003,0.000031,0.000033} 

4. Three equal masses and two increasing symmetrically 

(a) Hi = hq = 0.326 and /i 2 = 0.11 

C = {0.02,0.027,0.028,0.029,0.03} 

(b) n 1 = n Q = 0.15 and jj, 2 = 0.275 

C = {0.03,0.036,0.05,0.06,0.063} 

(c) /i! = [j, = 0.01 and ^2 = 0.485 

5. Non-equal masses 

(a) /i! = 0.195, n 2 = 0.3 and /i = 0.01 
Co = {0.02,0.0281,0.04,0.044,0.055} 

(b) /ii = 0.3, [/,% = 0.1 and Ho = 0.2 

C = {0.03,0.0346,0.05,0.055,0.06} 

(c) Hi = 0.35, H2 = 0.01 and Ho = 0.28 
C = {0.02,0.026,0.028,0.029,0.03} 

The values of C were chosen so as to represent each important region of the 
Szebehely ladder. Recall that Ri, R 2 , R3 and i? 4 , the rungs of the Szebe- 
hely ladder, represent boundaries for Co at which topology of the phase space 
changes. 

In order to analyze the behavior of the motion, for each (/x , fj,i, (J>2, C ) 
case, approximately 3000 integrations were performed for over 1 million time- 
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steps each. A total of 93x3000 = 279000 orbits were investigated numerically. 
Therefore the total integration time was well over 34 million time steps. 

Initially the system is in the 12 (DB) or 24 (SB) hierarchy state when 
r\ > r 2 , see figure (7.3a), and it is in the 12 (DB) or 13 (SB) hierarchy state 
when r 1 < r 2 , see figure (7.3b). Overall the CSFBP system is initially either in 
the 12 (DB), 13 (SB) or 24 (SB) hierarchy states. We do not initially place the 
system in a 14 hierarchy state, since it is equivalent to the 12 hierarchy state 
where only the numbering of the bodies differ. Both the cases 1) initially two 
larger masses in the centre and two smaller masses located outside ( r\ > r 2 ), 
and 2) two smaller masses in the centre and two larger masses located outside 
(ri < r 2 ) are dealt with. 

6.5 Hierarchical Stability of the CS5BP 

We have given an analytical stability criterion for the CS5BP, in chapter 5, 
and we have shown that the CS5BP is hierarchically stable for C >C crit . We 
have also shown that it is possible for the CS5BP to move from one hierarchy 
state to another for C < C crit . We aim to understand the behavior of the 
CS5BP during the transition from the phase space being connected to it being 
disconnected. 

To discuss the hierarchical stability of the CS5BP, we selected five C values 
i.e. C < Ri, Ri < C < R 2 , R 2 < C < R 3 , R 3 < C < R 4 and C > R 4 where 
i?4 is the critical value of C denoted by C crit for each mass ratio. The main 
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Table 6.1: Different Four and Five body systems investigated 



Diagram 


mass ratios 


Section number 


o-o — o-o 

o-O— O-o 

o O O o 

o O O o 


CSFBP 


Section 7.5.1 


Mi = Ai2 = 0.25, Mo = 
Mi = 0.05, ^2 = 0.45, Mo = 
Mi = 0.00495, M2 = 0.49505, Mo = 
Mi = 0.0004995, M2 = 0.4995095, Mo = 


O-O-O-O-O 


Equal mass case of the CS5BP 


Section 7.5.2 


Mi = M2 = Mo = 0.2 


O-O o o-o 
O-O O o-o 

o-o-O-o-o 


Four equal masses with a varying central mass 


Section 7.5.3 


22 475 1 

Mi - M2 - ioo » - ioo 

2 1 

Mi — M2 = 9, Mo = 9 
Mi - M2 - ioo) Mo - ioo 


O-O-O-O-O 
O-O-O-O-O 
O-O-O-O-O 


Three equal masses and two increasing 
symmetrically 


Section 7.5.4 


Mi = Mo = 0.326, M2 = 0.11 
Mi = Mo = 0.326, M2 = 0.11 
^1 = ^0 = 0.01, M2 = 0.485 


O-O-o-O-O 
O-o-o-o-O 

O o O o o 


Non-equal masses 


Section 7.5.5 


Mi = 0.195, M2 = 0.3, Mo = 0.01 

Mi = 0.3, M2 = 0.1, Mo = 0.2 
Mi = 0.35, M2 = 0.01, Mo = 0.28 
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reason for integrating for C >C crit was to validate that the integrator was 
working properly and the numerical investigation agreed with the analytical 
prediction; i.e that the phase space is disconnected for Co >C cr u and therefore 
the system is hierarchically stable. 

For each (/i , n\, (j, 2 , Co) set we integrated numerically approximately 3000 
different orbits for about 1 million time steps of integration time. In order to 
study the hierarchical stability for each (/i , /ii, Co) set we gathered data 
from the numerical integrations outputs on the number of hierarchy changes. 

We have constructed two kinds of tables for each (/x , A*i, ^2, Co) set, namely 

1. The frequency distribution of hierarchy changes: Columns represent fixed 
C values while rows represent the number of each type of hierarchy 
changes. The number of hierarchy changes were calculated by determin- 
ing the total number of hierarchy changes of each type which occurred 
throughout the many different integrations with the same (/i , fJ>i, 1^2, Co) 
set. 

2. Percentage of hierarchy changes: Columns represent fixed C values while 
rows represent percentage of hierarchy changes. 

Now we will discuss the hierarchical stability of the CS5BP for different mass 
ratios. First we will discuss the /x = case which is a special case of the 
CS5BP called the CSFBP in the following section. 
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6.5.1 Hierarchical stability of the CSFBP, = fi 2 and 

fi = 0) 

In this section we will discuss the /io = case of the CS5BP which is the 
CSFBP (Szell, 2003). Szell (2003) discussed a special case of the CSFBP for 
hierarchical stability. He discussed the case when r\ > r 2 or in other words 
when the larger masses were nearer to the centre of mass than the smaller 
masses, see figure (16.3a ). Therefore he covered half of the phase space. He 
also always started in the 12, double binary, or the 24, single binary, hierarchy 
state or using the CSFBP notation the 12 or 23 hierarchy states. Recall that we 
have a difference of notation with Szell (2003), the CSFBP 13 (DB) hierarchy 
state is 14 (DB) and its 23 (SB) is 24 (SB) in our case and vice versa, for 
more details see section 6.1. We will use our notation in the whole thesis 
and whenever a comparison is necessary we will give the CSFBP notation in 
brackets as against to ours. 

Our study of the /io = case will complete Szell's analysis of hierarchical 
stability of the CSFBP, as we will consider both T\ > r 2 and r\ < r 2 cases, see 
figure (16. 3p . 

Equal mass case of the CSFBP //i = /i 2 = 0.25, /io = (or in the 
CSFBP notation fi = 1) 

The critical values in this case are R\ = R 2 = 0.0286267 and R3 = = 
0.0457437. The results of integration are contained in Table (16.20 . 
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a) 




r 2 



Figure 6.3: The two possible initial states in the CSFBP a) r 1 > r 2 b) r 2 > r x 

There are no hierarchy changes for C > C crit which confirm the analytical 
results we obtained in chapter 6. The most likely hierarchy changes for Cq = 
0.01 and C = 0.028 are 12 -> 24 and 12 -> 13 and vice versa. The 12 -> 
24 hierarchy changes die out for higher Co values, while the 12 — > 13 remain 
dominant. For higher C values, the 13 — > 12 hierarchy changes joins the race 
for the most likely hierarchy changes and remains considerably close to the 12 
— > 13 hierarchy changes. All of the above changes are from SB to DB and vice 
versa. 

The most unlikely hierarchy changes are from double binary to double 
binary. These kind of hierarchy changes usually do not happen directly. There 
is usually a third kind of hierarchy change of intermediate hierarchy state which 
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takes the system in to a new configuration. If we look at the average number 
of hierarchy changes from different hierarchy states we come to the conclusion 
that the 12 is the most unstable while 14 is comparatively more stable than all 
others. For example there are nearly twelve hundred hierarchy changes from 
the 12 hierarchy state to another, which is the most, and there are a little 
above five hundred from the 14 hierarchy state to another. 

Hi = 0.05 and /j, 2 = 0.45, /x = or in CSFBP notation /j, = 0.1 

The critical values in this case are Ri = 0.012952, R 2 = 0.0138015, R 3 = 
0.0159447, i?4 = 0.0182837. The total number of hierarchy changes of each for 
different values of Cq are listed in Table (16.41) . The corresponding percentages 
of different hierarchy changes are given in Table (16.51) . 

For C = 0.019 > C crit the phase space is disconnected and therefore there 
are no hierarchy changes. At R 3 < C = 0.018 < C crit = R$ there are no 
hierarchy changes to or from the 24 hierarchy state which is exactly what the 
analytical stability criterion predicts i.e. there are no 12 — > 24, 14^24, 13 — > 
24, 24 — > 12, 24 — > 14 and 24 — > 13 hierarchy changes. 

The most likely hierarchy changes for small Cq values are the 12 — s- 24. And 
the most unlikely hierarchy changes are the 24 — ■> 14 and 14 — > 24. Overall 
hierarchy changes from 14 to any other hierarchy state for all Cq values are 
very small except for small Cq values eg. Cq = 0.018 where the 14 — > 12 and 
12 — >• 14 are the highest hierarchy changes. Hierarchy changes from double 
binary to double binary are very small as one would expect. 
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The main characteristics of this case is similar to those of the equal mass 
case. With small C , we have more hierarchy changes and as C increases the 
number of hierarchy changes decreases. The most unlikely hierarchy changes 
are still from double binary to double binary. 

The main difference in this case and the equal mass case is that we have 
fewer hierarchy changes. Therefore the system is hierarchically more stable. 
Moreover in the equal mass case we were not able to identify any, compara- 
tively, stable hierarchy states as all hierarchy states were equally unstable, but 
in this case we have some hierarchy states which can be considered hierarchi- 
cally stable. For example Hierarchy changes from the 14 (DB) hierarchy state 
is approximately 1 percent, which is very nearly stable. As we decrease the 
value of /ii we will be getting closer to a perturbed two+two body system and 
therefore we should expect a comparatively more stable situation. 

Hx = 0.00495 and /i 2 = 0.49505, /i = (or in CSFBP notation /i = 0.01) 

The critical values in this case are R x = 0.00828195, R 2 = 0.00838264, R 3 = 
0.00850907, i?4 = 0.00869902. The total number of hierarchy changes of each 
for different values of Co are listed in Table f)6.6p . The corresponding percent- 
ages of different hierarchy changes are given in Table (16.71) . 

There are no hierarchy changes with Co = 0.0087, since this value is greater 
than C cr it and therefore the phase space is disconnected. 

At C = 0.0086 the Szebehely constant is a little less than C crit and the 
phase space is partially disconnected. Therefore the hierarchy changes to or 
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from the 24 hierarchy state are forbidden. In table (16.61) it can be seen that 
there are no 12 -> 24, 14 -> 24, 13 -> 24, 24 -> 12, 24 -> 14 and 24 -> 13 
hierarchy changes, confirming numerically the analytical stability criterion. 

As we were expecting from our experience of the previous /ii = 0.05, case 
the number of hierarchy changes decrease further as we decrease the value 
of /ii. But the main characteristics remain the same. With small Co, 12 
— > 24, 24 — > 12, and 13 — > 24 hierarchy changes are dominant. This makes 
sense physically since as \ii is increased, the situation of two large masses 
in the centre forming a single binary with two smaller masses orbiting their 
centre of mass (eg the 24 hierarchy) will begin to dominate the system. We 
can see this growing dominance of 24 hierarchies visually in the projection 
of allowed motion for this case figure (16.41) which shows an increased region 
of allowed motion. When Cq is increased the number of hierarchy changes 
drop dramatically and the system is comparatively much more stable than the 
previous cases discussed so far. Again the most unlikely hierarchy changes are 
from double binary to double binary. 

The main difference between the previous case and this case is that the 
system is hierarchically more stable as the system is now closer to the perturbed 
two+two body system. 
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Figure 6.4: Projection found in the same manner as those in chapter 6 but for 
the case when fii = 0.00495 and fi 2 = 0.49505, /x = 

Hx = 0.0004995 and fi 2 = 0.4995095, /i = (or in CSFBP notation // = 

0.001) 

The critical values in this case are R\ = 0.00785938, R 2 = 0.00786963, R3 = 
0.00788154, i? 4 = 0.00789903. The total number of hierarchy changes of each 
for different values of Cq are listed in Table (I6.8p . The corresponding percent- 
ages of different hierarchy changes are given in Table (16.91) . 

There are only five hierarchy changes for the smallest Cq value, which are 
the only hierarchy changes for the 15,000 orbits integrated for one million time 
steps of integration time. These few hierarchy changes are from 13 — > 24, 24 
— > 13 and 24 — > 12. The very small number of hierarchy changes indicate that 
the system is generally hierarchically stable. This is so because the system is 
close to a slightly perturbed two body system. 
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Comparison with Szell's Analysis 

To compare our analysis of the Caledonian Symmetric Problem (CSFBP) with 
Szell's we have the following similarities and differences. Recall the difference 
in method of research with Szell (2003) which might produce some differences. 

1. Szell (2003) counts for only half of the phase space (ri > r 2 ). This might 
produce a bias in favor of some hierarchy states. 

2. Our hierarchy changing criterion is different from Szell (2003) as our 
criterion is based on the size of r\ and r 2 relative to each other, while 
Szell's criterion is also based on change from libration to rotation or vice 
versa of the angle a between Pi and P 2 . 

Szell (2003) number and percent of hierarchy changes are tabulated in 
Appendix 1. Recall that when making comparisons you must take account of 
Szell's different numbering of the four bodies from our numbering. 

Differences 

1. For //i = 0.25 (or CSFBP notation y, = 1), the 12 -> 24 and 24 -> 
12, (in CSFBP notation 12 — > 23 and 23 — ► 12), hierarchy changes are 
much higher in our case than in Szell's case. They are twenty percent to 
fourteen percent in our case and four percent to two percent in Szell's 
case. See tables ( 16. 2p and (16. 3p and table 1 and 2 in appendix 1. These 
differences are most likely because of the inclusion of the second half of 
the phase space into our research which almost doubles the number of 
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starting conditions for 12 hierarchy state. See appendix for Szell (2003) 
results. When the value of C increases, the hierarchy changes discussed 
above decrease in our case while it increases in Szell's case. 

2. For all values of the Szebehely constant Co in the [L\ = 0.05 (or CSFBP 
notation // = 0.1) case, the 12 — > 24 (CSFBP notation 12 — > 23)hierarchy 
changes are very high in our case while it is very small in Szell's case. The 
24 — > 13 (CSFBP notation 23 — ► 14) hierarchy changes are very small in 
our case and stand at a maximum of nine percent. In Szell's case they 
are very high and stands at thirty six percent maximum and ten percent 
minimum. See tables (16. 4p and (16. 5p and table 3 and 4 in appendix 
1. These differences are also most likely because of the inclusion of the 
second half of the phase space. The size of the region of real motion 
in the non e-qual mass cases for the 13 (CSFBP notation 14) hierarchy 
state is very small therefore we have fewer starting conditions for this 
hierarchy state and thus we will have fewer hierarchy changes to or from 
this hierarchy state. 

3. For /ii = 0.00495 (or CSFBP notation // = 0.01), the 13 -> 24 hierarchy 
changes are very small (maximum fifteen percent) in our case and be- 
comes zero near the critical value of the Szebehely constant Cq which is 
the prediction of the analytical stability criterion. In Szell's case these 
hierarchy changes are comparatively very high (maximum forty percent) 
and do not become zero near the critical value of the Szebehely constant 

193 



Cq which is against the prediction of the analytical stability criterion. 
See tables f !6.6j) and (16.71) . Szell (2003) do not give any reason for this, 
see appendix for their results. 

4. For //i = 0.00495 (Or CSFBP notation fj, = 0.01), the 24 -> 14 (CSFBP 
notation 23 — > 13) hierarchy changes are zero in our case while it goes 
up to fifteen percent in Szell's case. 

5. For /xx = 0.0004995 (Or p = 0.001), the 12 -> 13 (CSFBP notation 12 -> 
14) and 12 — > 24 (CSFBP notation 12 — > 23) hierarchy changes are zero 
in our case while they are non-zero in Szell's case and 12 — > 13 (CSFBP 
notation 12 — > 14) change increases from zero to thirty one percent when 
Co increases. 

Similarities 

1. the main characteristics of both the analysis remain the same. The 
number of hierarchy changes decrease as the value of Co increases and 
the number of hierarchy changes are zero for Co > C cr u 

2. For fix — /j,2 — 0.25, (in CSFBP notation fi = 1), and \i\ = 0.05, (in 
CSFBP notation fi = 0.1), the number of hierarchy changes from double 
binary to double binary are very small. Recall that the double binaries 
in our case are 12 and 14 while they are 12 and 13 in Szell's case. See 
tables |T2]) and fl6\4j) . 

3. For m = 0.05, in CSFBP notation /i = 0.1, the 24 -> 12 (CSFBP: 
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Table 6.2: Total number of hierarchical changes for the equal mass case of the 
CSFBP, n x = [M 2 = 0.25, fio = 0; or in CSFBP notation n=l 





C = 0.01 


C = 0.028 


0.031 


0.038 


0.042 


C = 0.046 


12 13 


179 


143 


150 


125 


104 





12 14 


5 





10 


22 


12 





12 24 


204 


162 


24 


33 


11 





13 12 


91 


9 


141 


119 


90 





13 14 


66 


10 


21 


20 


28 





13 24 


56 


18 


28 


27 


20 





14 12 


5 





11 


13 


19 





14 13 


60 


32 


24 


41 


38 





14 24 


80 


84 


36 


37 


26 





24 12 


146 


130 


32 


33 


21 





24 13 


130 


90 


32 


33 


41 





24 14 


85 


74 


26 


25 


26 





Total 


1022 


678 


509 


503 


410 






23 — > 12) hierarchy changes are very small for all Co- Also the 24 — > 
14 (CSFBP: 23 -> 13) and 14 -> 24 (CSFBP: 13 -> 23) are the most 
unlikely hierarchy changes in both the cases. See tables (I6.4p and (16.51) 

4. For /i! = 0.00495, in CSFBP notation \i = 0.01,the 14 -> 24 (CSFBP: 

13 — > 23) hierarchy changes are the most unlikely in both cases. 

5. For /i! = 0.0004995, in CSFBP notation /i = 0.001, the 13 -> 24 (CSFBP: 

14 -> 23) and 24 -> 13 (CSFBP: 23 -> 14) hierarchy changes are the most 
likely in both the cases. Also the number of hierarchy changes are very 
small with 12 -> 13, 14 -> 12, 14 -> 13 and 13 -»■ 14 (CSFBP: 12 -> 14, 
13 -> 12, 13 -> 14, and 14 -> 13) being zero. 
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Table 6.3: Percentage of hierarchical changes for the equal mass case of the 
CSFBP, t ji 1 = n 2 = 0.25, no = 0; or in CSFBP notation n = 1 





C = 0.01 


C = 0.028 


0.031 


0.038 


0.042 


C = 0.046 


12 13 


17.51 


21.09 


22.12 


18.44 


15.34 





12 14 


0.49 





1.47 


3.24 


1.77 





12 24 


19.96 


24.89 


3.54 


4.87 


1.62 





13 12 


8.9 


1.33 


20.8 


17.55 


13.27 





13 14 


6.46 


1.47 


3.1 


2.95 


4.13 





13 24 


5.48 


2.65 


4.13 


3.98 


2.95 





14 12 


0.49 





1.62 


1.92 


2.8 





14 13 


5.87 


4.72 


3.54 


6.05 


5.6 





14 24 


7.83 


12.39 


5.31 


5.46 


3.83 





24 12 


14.29 


19.17 


4.72 


4.87 


3.1 





24 13 


12.72 


13.27 


4.72 


4.87 


6.05 





24 14 


8.32 


10.91 


3.83 


3.69 


3.83 






Table 6.4: Total number of hierarchical changes for [i\ = 0.05, [ii = 0.45 and 
Ho = 0, in CSFBP notation /i = 0.1 





C = 0.012 


C = 0.013 


Co = 0.015 


C = 0.018 


C = 0.019 


12 13 


132 


98 


121 


9 





12 14 


1 


5 


5 


16 





12 24 


360 


348 


90 








13 12 


55 


61 


73 


2 





13 14 


2 


5 


3 








13 24 


136 


138 


46 








14 12 


5 


11 


3 


16 





14 13 


5 


8 


6 








14 24 





2 


2 








24 12 


78 


57 


22 








24 13 


77 


56 


22 








24 14 





3 


1 








Total 


851 


792 


394 


43 
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Table 6.5: Percentage of hierarchical changes for hi = 0.05, Hi = 0.45 and 
Ho = 0, in CSFBP notation \x = 0.1 





C = 0.012 


Co = 0.013 


C = 0.015 


C = 0.018 


C = 0.019 


12 13 


16 


12 


31 


21 





12 14 





1 


1 


37 





12 24 


42 


44 


24 








13 12 


6 


8 


19 


5 





13 14 





1 


1 








13 24 


16 


17 


12 








14 12 


1 


1 


1 


37 





14 13 


1 


1 


2 








14 24 








1 








24 12 


9 


7 


6 








24 13 


9 


7 


6 








24 14 


















Table 6.6: Total number of hierarchical changes for hi — 0.00495, Hi = 0.49505 
and no = 0, in CSFBP notation h = 0.01 





C = 0.008 


C = 0.0082 


C = 0.0085 


C = 0.0086 


C = 0.0087 


12 13 


41 


17 


4 


8 





12 14 


15 














12 24 


78 


3 


2 








13 12 


19 


5 


4 


4 





13 14 


2 














13 24 


46 


1 


3 








14 12 


9 














14 13 


12 














14 24 

















24 12 


47 





2 








24 13 


37 





4 








24 14 

















Total 


306 


27 


19 


12 
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Table 6.7: Percentage of hierarchical changes for [i\ = 0.00495, jj, 2 = 0.49505 
and yu = 0, in CSFBP notation fi = 0.01 





C = 0.008 


C = 0.0083 


C = 0.0085 


C = 0.0086 


Cq = 0.0087 


1 1 Q 






91 

__, J. 




n 


12 14 


5 














12 24 


25 


11 


11 








13 12 


6 


19 


21 


33 





13 14 


1 














13 24 


15 


4 


16 








14 12 


3 














14 13 


4 














14 24 

















24 12 


15 


4 


11 








24 13 


12 





21 








24 14 


















Table 6.8: Total number of hierarchical changes for fii = 0.0004995, fi 2 = 
0.4995005 and /i = 0, in CSFBP notation // = 0.001 





Co = 0.0078 


C = 0.00786 


Co = 0.00788 


Co = 0.00789 


Co = 0.0079 


12 13 

















12 14 

















12 24 

















13 12 

















13 14 

















13 24 


2 














14 12 

















14 13 

















14 24 

















24 12 


1 














24 13 


2 














24 14 

















Total 


5 
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Table 6.9: Percentage of hierarchical changes for /ii = 0.0004995, y«i = 
0.4995005 and fi = 0, in CSFBP notation fi = 0.001 





Go = U.UU to 


C-o = U.UU^oo 


C-o = U.UU^oo 


Go = U.UlWo9 


Go = U.UU ty 


12 13 

















12 14 














U 


12 24 

















13 12 

















13 14 

















13 24 


40 














14 12 

















14 13 

















14 24 

















24 12 


20 














24 13 


40 














24 14 


















6.5.2 Equal mass case of the CS5BP, Mi = A*2 = A*o = 0.25 

The critical values in this case are R\ — R% — 0.039222 and R3 = R4 = 
0.065551. The results are contained in Table (I6.10p . 

There are no hierarchy changes for Go > C cr u which is also the prediction 
of the analytical criterion. The most likely hierarchy changes for Cq = 0.03 are 
the 14 — > 24 and 24 — > 14. When we increase the value of C to 0.06 the 12 
— > 14 changes become the most likely hierarchy change which is from double 
binary to double binary. The most unlikely hierarchy changes are for all the 
Co values are those involving changes to and from 13. 

The main characteristic of this case is similar to that of the CSFBP i.e. the 
number of hierarchy changes decreases as the value of Cq increases; however 
there are some differences as follows. 
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Table 6.10: Total number of hierarchical changes for equal mass case of the 
CS5BP,//i = fa = /x = 0.25 





C = 0.03 


C = 0.06 


C = 0.07 


12 13 


14 








12 14 


54 


50 





12 24 


19 


10 





13 12 


8 


2 





13 14 


18 


6 





13 24 


10 


1 





14 12 


16 


9 





14 13 


17 


11 





14 24 


83 


15 





24 12 


18 


20 





24 13 


15 


7 





24 14 


67 


10 





Total 


339 


141 






1. There are fewer hierarchy changes overall than for the equal mass case 
of the CSFBP which means that the extra central mass has played a 
stabilizing role in terms of hierarchical stability and therefore the equal 
mass CS5BP is more stable than the equal mass CSFBP. 

2. The 13 — > 12 and 12 — > 13 were among the most likely hierarchy changes 
in the case of the equal mass CSFBP which clearly is the least in this 
case and the double binary to double binary hierarchy changes are most 
likely in the CS5BP case for large Co- See tables (16.21) and (16.101) 



6.5.3 Four equal masses with a varying central mass fiQ 

In this section we introduce a stationary mass to the centre of the equal mass 
CSFBP and will use it to discover the effect of changing its mass on the 
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Table 6.11: Percentage of hierarchical changes for equal mass case of the 
CS5BP,//i = /i 2 = hq = 0.25 





C = 0.03 


C = 0.06 


C = 0.07 


12 13 


4 








12 14 


16 


36 





12 24 


6 


7 





13 12 


2 


1 





13 14 


5 


4 





13 24 


3 


1 





14 12 


5 


6 





14 13 


5 


8 





14 24 


24 


11 





24 12 


5 


14 





24 13 


4 


5 





24 14 


20 


7 






hierarchical stability of the system. First we place a small mass at the centre 
of mass and then we increase it until the mass of the central body is much 
higher than the outer bodies. We discuss the following three different sets of 
mass ratios: 

1. /ii — /i 2 — 2 ^qq 5 and fi — Y5o : Four equal masses with a very small 
central mass; 

2. Hi — H2 — | and /io = |: Four equal masses with a comparatively larger 
central mass but still smaller than the outer bodies; 

3. /ii — ii 2 — an d A*o — jo^o'- Four equal masses with a large central mass 
and small outer bodies. 
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Four equal masses with a very small central mass — // 2 = 2 ^»J 5 
and no = ^) 

The critical values in this case are R 1 = 0.0292835, R 2 = 0.0293423, R 3 = 
0.0468416 and i? 4 = 0.0482373. Results are contained in tables fl6TT2|) and 

(EH. 

There are no hierarchy changes for Co > C CT u which is also the prediction 
of the analytical stability criterion. The most likely hierarchy changes for small 
C are 12 — > 13 and 13 — > 12. We have seen in the equal mass CS5BP that 
the single binary to double binary hierarchy changes and vice versa were small 
FOR large Co- These have increased after reducing the size of the central mass 
with in particular to and from 13 changes become larger. 

The main difference between this case and the equal mass CSFBP is that we 
have fewer hierarchy changes. Compared to the equal mass CS5BP, this case 
has more hierarchy changes. An increasing central mass appears to stabilize 
the system. We shall see this confirmed in the next two sections. 

Four equal masses with a slightly bigger central mass (/i! = fi 2 — § 
and no = §) 

The critical values in this case are R\ = 0.0368188 = R 2 , R3 = 0.062839 = -R4. 
Results are contained in tables (I6.14p and (16.15!) . 

We have increased the mass of the central body further but is still smaller 
than the four outer bodies. The total number of hierarchy changes decreased 
further but are still more than we had in the equal mass case of the CS5BP. 
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The introduction of the central mass is indeed playing a stabilizing role. The 
most likely hierarchy changes in this case are 14 — > 24 which were the most 
unlikely in the previous case with a smaller central mass. The most unlikely 
hierarchy changes are from 13 to any other hierarchy state for Co = 0.036. 

The main characteristic of both the cases remain the same. There are no 
hierarchy changes for Co > C cr u and the number of hierarchies decreases as 
the value of the Szebehely constant increases. 

Four equal masses with a very large central mass (/xi — [Li — ^ an ^ 
^° ~ loo-' 

The critical values in this case are R x = 0.0000300822 = R 2} R 3 = 0.000032294 = 
i?4. Results are contained in tables (16.161) and (16. 17ft . 

We have further increased the mass of the central body so that the central 
body is 96 times bigger than the outer bodies i.e. /i = 96/ii = 96/i 2 . This 
system is close to a planetary system with a star in the centre with four planets. 

The total number of hierarchy changes does not decrease much and is al- 
most the same as we had in the previous case suggesting a similar level of 
instability. This is deceiving if we do not look at the individual members of 
the table. Most of the hierarchy changes are coming from 14 — > 24 and vice 
versa. For small C (C = 0.00003), 14 — > 24 and 24 — > 14 changes con- 
tribute 74 percent to the total number of hierarchy changes and for large Co 
(C = 0.000033), they contribute more than 99.9 percent of the total hierarchy 
changes. See tables (16.161) and (I6.17p . 
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Table 6.12: Total number of hierarchical changes for [i\ — /j, 2 = and 
Ho = 0.01. Four equal masses with a very small central mass 





Co = 0.026 


C = 0.04 


C = 0.05 


12 13 


96 


60 





12 14 


30 


12 





12 24 


22 


20 





13 12 


83 


65 





13 14 


31 


25 





13 24 


26 


10 





14 12 


44 


31 





14 13 


63 


51 





14 24 


26 


5 





24 12 


28 


15 





24 13 


39 


40 





24 14 


19 


8 





Total 


507 


342 






There are no hierarchy changes from 12 — > 13, 12 — > 14, 13 — > 12, 14 — > 
12 and 14 — > 13 throughout the integrations carried out for this case. This 
system is comparatively much more stable than any other case of four equal 
masses with or without a stationary central mass. 

As we have seen from all other cases of four equal masses with a stationary 
central mass and without a stationary central mass, the number of hierarchy 
changes decreases as the central mass increases. In the equal mass case of the 
CSFBP with no central mass the number of hierarchy changes were highest 
then after the introduction of the central mass the number of hierarchy changes 
decreased as we increased the central mass. Therefore now we can conjecture 
that the larger the central mass the more stable the system becomes. 
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Table 6.13: Percentage of hierarchical changes for fx\ = ^2 = 10 q and 
0.01. Four equal masses with a very small central mass 





C = 0.026 


C = 0.04 


C = 0.05 


12 13 


19 


18 





12 14 


6 


4 





12 24 


4 


6 





13 12 


16 


19 





13 14 


6 


7 





13 24 


5 


3 





14 12 


9 


9 





14 13 


12 


15 





14 24 


5 


1 





24 12 


6 


4 





24 13 


8 


12 





24 14 


4 


2 






Table 6.14: Total number of hierarchical changes for /ii — ^2 — § and A*o 
Four equal masses with a slightly bigger central mass 





Co = 0.036 


C = 0.055 


C = 0.063 


12 13 


45 


4 





12 14 


35 


27 





12 24 


29 


19 





13 12 


22 


2 





13 14 


22 


10 





13 24 


22 


15 





14 12 


30 


20 





14 13 


32 


25 





14 24 


93 


37 





24 12 


34 


14 





24 13 


24 


12 





24 14 


75 


29 





Total 


462 


214 
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Table 6.15: Percentage of hierarchical changes for /ii — /i2 — § and fiQ = |. 
Four equal masses with a slightly bigger central mass 





C = 0.036 


C = 0.055 


C = 0.063 


12 13 


10 


2 





12 14 


8 


13 





12 24 


6 


9 





13 12 


5 


1 





13 14 


5 


5 





13 24 


5 


7 





14 12 


6 


9 





14 13 


7 


12 





14 24 


20 


17 





24 12 


7 


7 





24 13 


5 


6 





24 14 


16 


14 






Table 6.16: Total number of hierarchical changes for Hi = /i 2 = ^ and 
f^o = fo%- Four equal masses with a very large central mass 





C = 0.00003 


C = 0.000031 


C = 0.000033 


12 13 











12 14 











12 24 


32 








13 12 











13 14 


1 








13 24 


26 








14 12 











14 13 





1 





14 24 


136 


52 





24 12 


32 








24 13 


27 








24 14 


196 


93 





Total 


450 


146 






206 



Table 6.17: Percentage of hierarchical changes for Hi = ^2 = iM and /i = ^jf,. 
Four equal masses with a very large central mass 





C = 0.00003 


C = 0.000031 


C = 0.000033 


12 13 











12 14 











12 24 


7 








13 12 











13 14 











13 24 


6 








14 12 











14 13 





1 





14 24 


30 


36 





24 12 


7 








24 13 


6 








24 14 


44 


64 






6.5.4 Three equal masses and two increasing symmetri- 
cally 

In this section we will discuss a new case of the Caledonian Symmetric Five 
Body Problem. This is a completely different situation from those discussed 
earlier and little comparison with the previous cases will be made. Here we 
take the masses of three of the bodies to be equal and two of them to be varying 
symmetrically from being very small to very large and effectively approaching 
a situation of a perturbed two body problem. The three bodies with equal 
masses will comprise of two outer bodies and the central body i.e. Hi = Ho- 
H2 will increase to the largest masses from being the smallest. We will discuss 
the following three cases. 

1. Hi = no = 0.326 and /i 2 = 0.11 

2. Hi = Ho = 0.15 and /i 2 = 0.275 
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3. fj,i = /i = 0.01 and /i 2 = 0.485 
Hi = Ho — 0.326 and /i 2 = 0.11 

The critical values in this case are Ri = 0.0269644, R 2 = 0.0270025, R 3 = 
0.0286002,^4 = 0.0294447. Results are contained in tables (IBTTSl) and (1091) . 

There are no hierarchy changes for Cq = 0.03. Since this value is greater 
than the critical value, this fact shows that the numerical integrations follow 
well the analytical stability criterion. 

The most unlikely hierarchy changes are 14 — > 24. The 12 — > 13, 13 — > 
12, 13 — > 24, 24 — > 13 and 24 — > 14 hierarchy changes are also very small for 
small Cq (Co = 0.02). The 14 — > 24 and 13 — > 24 hierarchy changes remains 
the smallest for all values of Cq. 

For small Cq = 0.02, 12 — ► 14 and 14 — * 12 hierarchy changes are the most 
likely and remain amongst the highest throughout the integrations carried 
out for this case of the CS5BP. These hierarchy changes remain the highest 
throughout the integration. 

At Co = 0.029 the Szebehely constant is near C cr u and the phase space 
is partially disconnected. Therefore the hierarchy changes to or from the 13 
hierarchy state are not allowed. This is numerically confirmed in table (I6.16P 
where there are no 12 -> 13, 13 -> 12, 13 -> 14, 13 -> 24, 14 -> 13 and 24 -> 
13 hierarchy changes. 

The main characteristic of this case is similar to all other cases of the 
CS5BP. With small Cq we have more hierarchy changes and as Cq increases 
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the number of hierarchy changes decrease. 
Hi = Ho = 0.15 and pi 2 = 0.275 

The critical values in this case are Ri = 0.0364159, R 2 = 0.0365447, R3 = 

0. 0578382,^4 = 0.0622745. Results are contained in tables (E2QD and <K21\\ . 

In this case we increase the value of H2 to 0.275 and consequently the value 
of Hi decreases and hence we have Hi > Hi — Ho- 

There are no hierarchy changes for Co > C cr u. For Co = 0.06 which is near 
the critical value there are no hierarchy changes from and to the 24 hierarchy 
state which is also the prediction of analytical stability criterion. The most 
likely hierarchy changes for small Co (Co = 0.03) are 14 — > 24 and 12 — > 24. 
When Co increases the 14 — > 24 hierarchy changes decrease and the 12 — > 24 
changes remains the most likely until it dies out for large Co [Co = 0.06). The 
most unlikely hierarchy changes are 12 — > 14 and 13 — > 14. 

The main characteristic of this case and all other cases remain the same 

1. e. the number of hierarchy changes decreases as Cq increases. 
Hi — Ho — 0.01 and /j, 2 = 0.485 

The critical values in this case are R x = 0.00968676, R 2 = 0.00988033, R 3 = 
0.0102131,^4 = 0.0107386. Results are contained in tables fl632|) and (7.23). 

In this case we further decrease the value of Hi an d Ho an d consequently 
the value of Hi increases. As Hi — Ho ^ Hi therefore this can be considered a 
perturbed two body problem. 

There are no hierarchy changes for any value of Co we have investigated. 

209 



Table 6.18: Total number of hierarchical changes for /ii = /i = 0.326, jj2 = 
0.11 





Go = 0.02 


S~1 r\ HOT 

Go = 0.027 


Go = 0.028 


S~1 r\ AHA 

Go = 0.029 


Go = 0.03 


12 13 


15 


7 


11 








12 14 


108 


24 


21 


15 





12 24 


45 


26 


24 


15 





13 12 


22 


6 


6 








1 o -1/1 

lo 14 


7o 


7 


10 








1 ^ 9zL 


1 1 


Q 
d 


Q 
O 


n 

U 


n 

U 


14 12 


103 


24 


28 


21 





14 13 


84 


7 


3 








14 24 


6 


2 


1 


1 





24 12 


35 


14 


8 


4 





24 13 


15 


6 


8 








24 14 


14 


5 


5 


7 





Total 


536 


130 


128 


63 






See table ( 16.221) . We have investigated about 15,000 orbits each for 1 million 
time steps of integration time and there are no hierarchy changes. Therefore 
this is the most stable of the CS5BP systems we have discussed so far in this 
chapter. 

6.5.5 Non-Equal mass cases of the CS5BP 

The cases of the CS5BP we have discussed so far had at least two or more 
of the mass ratios equal which included, fii = /i 2 and fi = 0, /ii = /i 2 = fi , 
fii = fi2 and /io 7^ 0, and /ii = /io and /io ^ 0. In this section we discuss the 
cases of the CS5BP where none of the mass ratios are equal. We will discuss 
the following three sets of mass ratios: 



1. /ii = 0.195, /i2 = 0.3 and no = 0.01 

2. /i! = 0.3, fi 2 = 0.1 and /i = 0.2 
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Table 6.19: Percentage of hierarchical changes for jii — /j, — 0.326, /i 2 = 0.11 





C = 0.02 


C = 0.027 


Co = 0.028 


Co = 0.0029 


Co = 0.03 


12 13 


3 


5 


9 








12 14 


20 


18 


16 


24 





12 24 


8 


20 


19 


24 





13 12 


4 


5 


5 








13 14 


15 


5 


8 








13 24 


2 


2 


2 








14 12 


19 


18 


22 


33 





14 13 


16 


5 


2 








14 24 


1 


2 


1 


2 





24 12 


7 


11 


6 


6 





24 13 


3 


5 


6 








24 14 


3 


4 


4 


11 






Table 6.20: Total number of hierarchical changes for [i\ — // = 0.15, yU 2 = 
0.275 





Co = 0.03 


Co = 0.036 


Co = 0.05 


Co = 0.06 


Co = 0.063 


12 13 


38 


24 


16 








12 14 


26 


13 


74 


14 





12 24 


67 


70 


99 








13 12 


20 


20 


8 








13 14 


14 


7 


38 


3 





13 24 


47 


26 


27 








14 12 


24 


17 


72 


22 





14 13 


30 


18 


36 


4 





14 24 


86 


42 


16 








24 12 


53 


28 


37 








24 13 


44 


17 


34 








24 14 


47 


35 


24 








Total 


496 


317 


480 


43 
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Table 6.21: Total number of hierarchical changes for [i\ — /j, — 0.15, [ii = 
0.275 





C = 0.03 


C = 0.036 


Co = 0.05 


Co = 0.06 


Co = 0.063 


12 13 


8 


8 


3 








12 14 


5 


4 


15 


33 





12 24 


14 


22 


21 








13 12 


4 


6 


2 








13 14 


3 


2 


8 


7 





13 24 


9 


8 


6 








14 12 


5 


5 


15 


51 





14 13 


6 


6 


8 


9 





14 24 


17 


13 


3 








24 12 


11 


9 


8 








24 13 


9 


5 


7 








24 14 


9 


11 


5 









Table 6.22: Total number of hierarchical changes for /jLi — /j, — 0.01, /i 2 = 
0.485 





C = 0.009 


Co = 0.0098 


Co = 0.01 


Co = 0.0106 


C = 0.0108 


12 13 

















12 14 

















12 24 

















13 12 

















13 14 

















13 24 

















14 12 

















14 13 

















14 24 

















24 12 

















24 13 

















24 14 

















Total 
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Table 6.23: Percentage of hierarchical changes for jii = /iq = 0.01, /i2 = 0.485 





C = 0.009 


C = 0.0098 


C = 0.01 


C = 0.0106 


C = 0.0108 


12 13 

















12 14 

















12 24 

















13 12 

















13 14 

















13 24 

















14 12 

















14 13 

















14 24 

















24 12 

















24 13 

















24 14 


















3. /i! = 0.35, n 2 = 0.01 and /i = 0.28 
li x = 0.195, fi 2 = 0.3 and /i = 0.01 

The critical values in this case are Ri = 0.0280812,^2 = 0.0283641,^3 = 
0.0439109,^4 = 0.0469418. Results are contained in tables fl6T24j) and (16T25|) . 

There are no hierarchy changes for Co = 0.055 > C cr i t as was predicted 
by the analytical stability criterion. For Co = 0.044 there are no hierarchy 
changes to or from the 24 hierarchy state which is exactly what the analytical 
stability criterion predicts. 

The most likely hierarchy changes are always to or from the 24 hierarchy 
state until it becomes disconnected from all other hierarchy states. For small 
Co (C = 0.02), the 24 — > 13 hierarchy changes are the highest, for C = 0.0281, 
the 13 — > 24 hierarchy changes is the highest and for large C (C = 0.04), 
the 24 — > 13 hierarchy changes are again the highest which means that the 24 
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hierarchy state is the most unstable for small values of the Szebehely constant, 
C . Also the most number of hierarchy changes are from the 24 hierarchy state 
which is another reason for branding the 24 hierarchy state the most unstable 
for small C . 

The most unlikely hierarchy changes are from the 14 hierarchy state for 
C = 0.02, 0.028 and 0.04. This suggests that the 14 hierarchy state is the 
most stable hierarchy state but perhaps it is the result of not starting in the 
14 hierarchy state. Fewer starts in the 14 hierarchy mean fewer changes from 
it are possible. At C large, C = 0.044, the 14 hierarchy state appears to be 
contributing more than 50 percent to the total hierarchy changes and turns 
out to be unstable. For large values of the Szebehely constant near the critical 
value the 24 hierarchy state is the most stable opposite to the fact that it is 
the most unstable situation for small C . 

The main characteristic of this case and all other cases of the CS5BP dis- 
cussed so far remain the same. There are no hierarchy changes for Co > C cr u 
and the number of hierarchy changes decrease as the value of Co increases. 
There is a sudden decrease in the total number of hierarchy changes near C cr u 
but this is not unexpected as at this point the 24 hierarchy state which was 
the main contributor disconnects and therefore the total number of hierarchy 
changes drop. 
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/ii = 0.3, /x 2 = 0.1 and fi = 0.2 

The critical values in this case are Ri = 0.0334268, R2 = 0.0336685, -R3 = 
0.0526554, R A = 0.0577769. Results are contained in tables (167261) and (167271) . 

We increase the value of fio from 0.01 to 0.2 from the last case and the 
total number of hierarchy changes drop dramatically. As we said earlier, the 
bigger the central mass, the more stable the system would be and this behavior 
compliments our statement. 

There are no hierarchy changes for C > C crit . For C large, C = 0.055 < 
C cr it the number of hierarchy changes are zero which is another reason to 
believe that this system is comparatively more stable than the [i\ = 0.195, 
H2 = 0.3 and /io = 0.01 case discussed earlier. The number of hierarchy 
changes from the 14 hierarchy state are very small throughout the integrations 
performed for this set of mass ratios and becomes zero for Co = 0.05 and 
therefore 14 is the most stable hierarchy state. Recall that 14 was the most 
stable hierarchy state in the previous case too. The 24 — > 12 and 24 — > 13 
hierarchy changes are also very small and is the most unlikely for Cq = 0.03 
and C = 0.0346. See table (167261) . 

The most likely hierarchy changes for Co = 0.03 and Co = 0.0346 are 12 
— > 13. For Co = 0.05 the highest number of hierarchy changes are from 13 to 
24. See table (167261) . 
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fit = 0.35, fi 2 = 0.01 and fi = 0.28 

The critical values in this case are Ri = 0.0243191, R2 = 0.0251055, -R3 = 
0.0361804, R A = 0.0408993. Results are contained in tables (l6\28l) and (l6"29l) . 

The total number of hierarchy changes further dropped as the value of fi 
increased and fi-i decreased. There are no hierarchy changes for Co = 0.029 
and Cq = 0.03 > C cr i t . There are no hierarchy changes from the 14 hierarchy 
state except one 14 — > 24 hierarchy changes for Co = 0.028, see table (16.281) . 
Therefore the 14 hierarchy state continue to be the most stable for all the three 
sets of mass ratios discussed in this section. There are also no 12 — > 14 and 14 
— > 12 hierarchy changes which is the double binary to double binary. For this 
particular set of mass ratios the double binary hierarchy states i.e. 12 and 14 
appear to be more stable as there are very few hierarchy changes from 12 as 
well. 

The most likely hierarchy changes are to or from the 24 hierarchy state as 
the size of the region of real motion is larger than the others for 24 hierarchy 
state. The 24 — > 13 and the 13 — > 24 hierarchy changes are the highest for 
C = 0.02 and for C = 0.026. For C = 0.028 the 24 -> 12 and 12 -> 24 
hierarchy changes are the highest. 

It is very interesting to see that in all the three sets of mass ratios discussed 
in this section, the 14 hierarchy state which is a double binary was the most 
stable and the 24 hierarchy state, which a single binary, was the most unstable. 

The main characteristic of all the cases of the CS5BP discussed in this 
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Table 6.24: Total number of hierarchical changes for [i\ = 0.195, //o = 0.01, 
and ji 2 = 0.3 





C = 0.02 


C = 0.0281 


C = 0.04 


C = 0.044 


C = 0.055 


12 13 


58 


50 


56 


14 





12 14 


14 


7 


22 


12 





12 24 


48 


40 


50 








13 12 


56 


34 


21 


2 





13 14 


9 


7 


6 








13 24 


38 


67 


17 








14 12 


10 


12 


14 


20 





14 13 


16 


12 


9 


10 





14 24 


6 


9 











24 12 


85 


53 


33 








24 13 


93 


60 


49 








24 14 


9 


5 











Total 


442 


356 


277 


58 






chapter remains the same. There are no hierarchy changes for the Szebehely 
constant greater than the critical value and the number of hierarchy changes 
decreases as C is increased. 

6.6 Conclusions 



The numerical investigations shows that the introduction of a small stationary 
mass to the centre of mass of the Caledonian Symmetric Four Body Problem 
(CSFBP) plays a stabilizing role. The number of hierarchy changes decreases 
as we increase the value of /x i.e. the central mass. This can be applied to a 
system of four planets and one star at the centre of mass of the system. We 
have seen that such system is hierarchically stable, see for example the analysis 
and tables provided for ^ — /i 2 — ^ and no — j^. 

The numerical investigations support the analytical predictions i.e. that the 
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Table 6.25: Percentage of hierarchical changes for [i\ = 0.195, no = 0.01, and 
/i 2 = 0.3 





C = 0.02 


C = 0.0281 


C = 0.04 


C = 0.044 


Co = 0.055 


12 13 


13 


14 


20 


24 





12 14 


3 


2 


8 


21 





12 24 


11 


11 


18 








13 12 


13 


10 


8 


3 





13 14 


2 


2 


2 








13 24 


9 


19 


6 








14 12 


2 


3 


5 


34 





14 13 


4 


3 


3 


17 





14 24 


1 


3 











24 12 


19 


15 


12 








24 13 


21 


17 


18 








24 14 


2 


1 












Table 6.26: Total number of hierarchical changes for [i\ = 0.3, /io = 0.2, and 
/i 2 = 0.1 





C = 0.03 


Co = 0.0346 


C = 0.05 


C = 0.055 


Co = 0.06 


12 13 


84 


55 


1 








12 14 


16 


16 


4 








12 24 


12 


18 


18 








13 12 


26 


12 


5 








13 14 


19 


19 











13 24 


4 


10 


30 








14 12 


9 


6 











14 13 


19 


19 











14 24 


19 


28 











24 12 


8 


5 


6 








24 13 


8 


13 


5 








24 14 


26 


35 











Total 


250 


246 


69 
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Table 6.27: Percentage of hierarchical changes for Hi = 0.3, /z = 0.2, and 
H2 = 0.1 





C = 0.03 


C = 0.0346 


C = 0.05 


C = 0.055 


C = 0.06 


12 13 


34 


24 


1 








12 14 


6 


7 


6 








12 24 


5 


8 


26 








13 12 


10 


5 


7 








13 14 


8 


8 











13 24 


2 


4 


43 








14 12 


4 


3 











14 13 


8 


8 











14 24 


8 


12 











24 12 


3 


2 


9 








24 13 


3 


6 


7 








24 14 


10 


15 












Table 6.28: Total number of hierarchical changes for [i\ = 0.35, Ho = 0.28, and 
H2 = 0.01 





C = 0.02 


C = 0.026 


C = 0.028 


C = 0.029 


C = 0.03 


12 13 





3 











12 14 

















12 24 


8 


6 


2 








13 12 





3 











13 14 

















13 24 


17 


3 











14 12 

















14 13 

















14 24 








1 








24 12 


14 


6 


1 








24 13 


17 


4 











24 14 


2 





1 








Total 


58 


25 


5 
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Table 6.29: Percentage of hierarchical changes for [i\ = 0.35, /x = 0.28, and 
n 2 = 0.01 





C = 0.02 


C = 0.026 


C = 0.028 


C = 0.029 


C = 0.03 


12 13 





12 











12 14 

















12 24 


14 


24 


40 








13 12 





12 











13 14 

















13 24 


29 


12 











14 12 

















14 13 

















14 24 








20 








24 12 


24 


24 


20 








24 13 


29 


16 











24 14 


3 





20 









phase space is disconnected for C > C crit . In some cases the analytical curves 
predicted that the phase space is partially disconnected eg. for R 3 < Co < R^- 
This is also supported by the numerical investigations. 

The stability of the four and five body systems increases as the value of 
pii or n 2 decreases. With small /ii or /j 2 the system is close to a three body 
system perturbed by two bodies of negligible mass. In the case of /xi = /io or 
fi 2 = A*o ver Y small as compared to the other mass ratio then, the system is 
close to a perturbed two body system. 

Numerical investigations of all the CS5BP systems discussed in this chapter 
reveal that as the value of the Szebehely constant increases the number of 
hierarchy changes decreases and hence their hierarchical stability improves. 
Therefore it can be stated for the Szebehely constant near zero the system will 
be very unstable, for increasing Co it will become more stable and for greater 
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than the critical value it will be hierarchically stable for all time. 

We have seen that the studies of the CSFBP and the CS5BP have given 
valuable information to understand better the general problems. Therefore 
it is of interest to generalize the CS5BP system in the next chapter for any 
number of bodies. 
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Chapter 7 

The Caledonian Symmetric N 
Body Problem 

In chapter 5 we investigated the Caledonian Symmetric five Body Problem 
(CS5BP), which is a symmetrically reduced five body problem containing all 
possible symmetries in the phase space. We have found an analytical stability 
criterion for the CS5BP. In chapter 5 we numerically verified this stability cri- 
terion and have shown that C , the Szebehely constant, plays a very important 
role in determining the hierarchical stability criterion of the CS5BP. 

In this chapter we generalize the CS5BP to the CSNBP i.e the Caledonian 
Symmetric N Body Problem in such a way that it can be easily used for any 
number of bodies such as three, four (CSFBP) and five (CS5BP) etc. In section 
7.1 we define the CSNBP. The equations of motion and Sundman's inequality, 
the key to the stability criterion, are given in sections 7.2 and 7.3 respectively. 
We then derive in sections 7.4 and 7.5 the analytical stability criterion for the 
CSNBP. The conclusions are given in section 7.6. 
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7.1 Definition of the Caledonian Symmetric N 
Body Problem (CSNBP) 

The CSNBP is formulated by using all possible symmetries. The main feature 
of the model is its use of two types of symmetries. 1. past-future symmetry 
and 2. dynamical symmetry. Past future symmetry exists in an n-body system 
when the dynamical evolution of the system after t = is a mirror image of 
the dynamical evolution of the system before t = 0. It occurs whenever the 
system passes through a mirror configuration, i.e. a configuration in which the 
velocity vectors of all the bodies are perpendicular to all the position vectors 
from the centre of mass of the system (Roy and Ovenden, 1955). 

Let us consider In + 1 bodies P , Pi, P 2 , ■ ■ ■, P 2n , of masses m , mi, m 2 , 
■ ' ^2n respectively existing in three dimensional Euclidean space. The radius 
and velocity vectors of the bodies with respect to the centre of mass of the 
five body system are given by and Tj respectively, % — 0, 1, 2, • • • , n. Let the 
centre of mass of the system be denoted by O. The CSNBP has the following 
conditions: 

1. All 2n + 1 bodies are finite point masses with: 

mi = m n+i , % = 0, 1,2, ...,n. (7.1) 

2. Pq is stationary at O, the centre of mass of the system. Pi and Pi+ n are 
moving symmetrically to each other with respect to the centre of mass 
of the system. Likewise Pj and P n +j are moving symmetrically to each 
other, see figure (JT7TJ) . Thus 
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V» = I-i = -Ti+n, Vj = Tj = -Tj+n, V = r = 0, (7.2) 

This dynamical symmetry is maintained for all time t. 

3. At time t = the bodies are collinear with their velocity vectors perpen- 
dicular to their line of position. This ensures past-future symmetry and 
is described by: 

?i x Vj = 0, Vi Vi = 0, Vj ■ Vj = (7.3) 

Figure (17.11) gives the dynamical configuration of the CSNBP. 

It is useful to define the masses as ratios to the total mass. Let the total 
mass M of the system be 

n 

M = 2 mi + m (7.4) 

Let Hi be the mass ratios defined as /ij = ^ for % = 0, 1, 2, • • ■ , 2n. Equation 
(17.41) then becomes 

n 

M = (2^2 f a + fi )M (7.5) 

i=l 

Thus 

E". + f = 5' (7 ' 6) 

1=1 

and 

< /io < 1, 0</ii<0.5, i = l,2,---,n (7.7) 
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We simplify the problem yet further by studying solely the coplanar CSNBP, 
where the radius and velocity vectors are coplanar. Figure (17. ip gives the dy- 
namical configuration of the coplanar CSNBP at some time t. In the next 
sections we will derive an analytical stability criterion for the CSNBP. This 
criterion was applied to the CS5BP for a wide range of mass ratios in chapter 



Let there be 2n + 1 bodies of masses m , mi, m.2, m 2 „_i, m 2ri ,. Then their 
equations of motion may be written as 



where V« = i^- + j^- + k^-, i, j, k, being unit vectors, along the rectangular 
axes O x , O y , O z respectively, Xi, yi, Zi being the rectangular coordinates of body 
Pi and O being the center of mass of the system. 
The force function U is given by 



5. 



7.2 The Equations of Motions 



rriifi = ViU, 



i = 0, 1, 2, ...k where k = 2n 



(7.8) 



k 



k 



k 





(7.9) 



where 



Then the energy E of the system may be written as 



E = T-U 



(7.10) 
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where T is the kinetic energy given by 

1 N 

T=-Y j m l \r l \\ (7.11) 

i=i 

Note: m remains stationary at the centre of mass so that there is no contri- 
bution to the kinetic energy, angular momentum and moment of inertia from 
m . The angular momentum C is given by 

N 

C=^2m i v i xr i . (7.12) 

%=\ 

We may also write the moment of inertia of the system / as 

N 

1 = ^2^1 (7.13) 

i=\ 

Now we may introduce the following symmetries. 

1. By the Roy-Ovenden mirror theorem the orbital history of the system 
after t = is a mirror image of the history after t = 0, given that their 
velocity vectors are perpendicular to their relative radius vectors at t — 0. 

2. The dynamic symmetry at any time t. Divide the 2n bodies into two sets 
of bodies P a , a = 1, 2, ...n, and Pp, (3 — n + 1, n + 2, ...2n — 1, 2n and let 
P a in the a set have mass m a and position and velocity vectors r a and 
r a at time t. Let the body Pp in the (5 set have mass m a and position 
and velocity vectors — r a and — r a at time t. The (2n + l) th body has 
mass m and position and velocity vectors r = and r = at time t. 

Then the kinetic energy T, the angular momentum C and the moment of 
inertia / may be written as 
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i,n+j 



Figure 7.1: Model of the Caledonian Symmetric N Body Problem (CSNBP) 

n 

T = Y,m l \r l \ 2 J (7.14) 

i=i 

n 

C =2^m i r i x Yi . (7.15) 

n 

1 = 2^^. (7.16) 

i=l 

Now consider the force function given by f)7.9p . The figure (17.1 j) defined by 
any two symmetric pairs and a fifth body at the centre of mass is always a 
parallelogram of changing shape and orientation. Now 



ij ' n+i,n+ j ; ' j,n+i ' i,n+j 



(7.17) 
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also 

r i,n+i — 2 r i] r j,n+j = 2Vj. (7- 18) 

Consider A's PiPjO, PjOP n+i . Let ZPiOPj = 9. 
Then 

rf • = r\ + r 2 — 2rj7j cos 6. (7-19) 



Also 



r-,n + i = r 2 n+i + rj + 2r n+t r 3 cos 9. (7.20) 



But r n+i = r iy so that 



r 2 „ +l = 2(rf + rj) - rj. (7.21) 



as 



Therefore the six mutual radius vectors in the parallelogram may be written 



PiP n+i = 2n; (7.22) 



P 3 P n+l = P t P n+J = ^2 (r? + rj)-rj; (7.23) 

Pn+iPn+j = PiPj = r ij- (7.24) 

Now using the above symmetries the force function can be written as 



i=i ^ ,=i ,=i \ r « ^2 (rf + r 2 ) - r% ) i= i 



m 8 
r,- 



(7.25) 
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7.3 Sundman's Inequality 

Sundman's inequality may be written as (see Orbital Motion by Archie Roy 
(2004)) 

T > i (7.26) 

Now 

E = T-U, (7.27) 

so that 

U + E>^. (7.28) 
Let Eq = —E. Then for real motion, we must have 

U > ^ + E . (7.29) 

Using (17351) . (17391 becomes 

( rl 2 rl n / i i \ n 

c 2 

^ 2 + ^o, where i 7^ j and j < i. (7.30) 

Let M be the total mass of the system, so that 

n 

M = 2^m i + m Q . (7.31) 

i=l 

Let 

Tfl ' 

[Mi = —, so that < fii < 1. (7.32) 

Then ( 17.311) becomes 

M= ^2^/^ + ^ Af, (7.33) 
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thus 



I> + f = 5- (7-34) 

1=1 

Thus the Sundman's inequality becomes 

i 2 n n / i i \ n 

GM2 [J2t + 2 ££** ^ + 2 + 2 /^ 

<=i i=i j=i \ 13 v 2 ^ + r ^ _ 1=1 



^ i rra — — "2 + ^o, where i ^ j and j < i. (7.35) 



Let us now introduce dimensionless variables Pi,pij and a new dimensionless 
constant Co, called the Szebehely constant Steves and Roy (2001) defined by 
the following equations 

Pi = Pij = J^r.. (7.36) 

C ° = GH^' (7 ' 37) 
where E ^ 0. Then Sundman inequality takes the following form 



n 2 n n / i i \ n 



Pi 

c 

- ami V™ 2 + X > wherei ^ j and j < i. (7.38) 

The geometry of triangles places the following constraint on the radius vectors. 

\n -rj\ < r i:j <ri + Tj. (7.39) 

Then from (I7.37P and (I7.39P we get the following relation between the new 
variables, 

\Pi - Pj\ < Pij < Pi + p r (7 AO) 
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7.4 Regions of motion in the CSNBP 

In this section we construct explicit formulae for the boundary surface of real 
motion. 

It is useful to parameterize the surface in terms of variables y^ (i — 1, 2, • • • , n) 
and Xij. We introduce two more variables namely yi and Xij which are defined 
as follows 

yi = —, here p n is the largest value of pi, for all i. (7-41) 

Pn 

Xij = & (7.42) 

Pn 

In the new variables, the kinematic constraints of (17.401) become 

\Vi -Vj\< <yi + Vj. (7.43) 

Then Sundman's inequality takes the following form after introducing the 
above new variables, 



i / 1 n 2 n n 1 \ n 

P» \ 2 U V* UU \ X « J2 (yf + yf)-x% I 



- TTT^ 2 + X ' wnere * ^ 3 and 3 < i- (7-44) 

4p^ E i= i IHVi 

Hence 

—A>\b + 1, (7.45) 

Pn Pt 

where 

^(?/i, 2/2, ■ • • , 2/n, = ^5Z~ + 2 5Z5Z — 



+ 2/^0 V and j<i (7.46) 
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and 

C 

B(yi,y 2 ,---,y n ) = ■: F^n ° 2- ( 7 - 47 ) 

Taking the equality sign in (17.451) which defines a boundary between real and 
imaginary motion. 

—A = \b + 1, (7.48) 

Pn Pi 

this implies that 

p 2 n - A Pn + B = 0. (7.49) 
After solving this quadratic equation and then simplifying we get the following 



pn(Vu V2, ■ • • , Vn, Xn) = ^ ^1 ± y 1 - jp^t ~ y2 ) ■ ( 7 ' 5 °) 

Define C(y u y 2 , - ■ • , y n , x 12 ) by 



n 



C(yi, y 2 , ■ ■ ■ , y n , xi 2 ) = A 2 /i^ 2 , (7.51) 



1=1 

so that 



Pn(yi, 2/2, • • • , 2/n, xia) = ±ll±J\-<£\. {7J)2) 



or 



1/ C(yi,y 2 , ■ ■■,y n ,x 12 ). 

2 l v sr=i^ 2 



Pn{yi,y 2 ,- ■ ■ ,y n ,xi2) = C a / ^ — —2 ) x 



The value of p n depends explicitly on C(yi, y 2 , • • • , y n , x\ 2 ), if 

1. C(yi, y 2 , ■ ■ ■ , y n , X\ 2 ) > C , there are two real roots for p n 

2. C(yi, y 2 ,---, y n , x± 2 ) = Co, there is a double real root for p n 

3. C(yi, y 2 , ■ • ■ , y n , Xi 2 ) < C , there are two imaginary roots for p n 
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7.5 Maximum and Minimum extensions of the 
real motion projected in the phase space 

We define a new function Ky, given by 

(i i \ 

= innj — + = . (7.54) 

We may write 

= fiifij Wij, (7.55) 

r 2 (yf + yf) 



where 

1 1 



Wn = — + . (7.56) 



and 

x 



u l3 = m (7.57) 

V 2 (yf + yf) 



Therefore f)7.5ip becomes 



C(yx,y 2 , • • • , 2/n, ^12) = I X] Mi 

\i=i 

V i = l 1=1 J = l 4 = 1 / 



(7.5* 



where % 7^ j and j < i. 



By the close inspection of W 7 ^ and u;^ we find the following extreme values 



of Wi 



2y/2 at Uij = 4= 



2 v^maajjj/i } [yf+y 



I. In 



when u>ij is either —^l^L= or 



(7.59) 
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We know from (17.551) that 



K ijmln = mm 1 =w ijmia . (7.60) 

As Wa . =1\pl therefore 

b J mm y 

1 

K ijmin = 2 PiPj I =■ (7-61) 

sjvl + v] 

And hence 

c,„„ = cu = E ™ 2 * U E 7 + 4 E E + ^ E 7 W) 

where i ^ j and j < i. 
Therefore the equations giving minimum projections are 



1 f C rn (y 1 ,y 2 ,---,y n,x 12 ) 

2 IV ELi /W 



Pn\VX, V2, ■ ■ ■ , Vn, Xi2) = ~ \ \j ^3 — ~ 2 ) X 



1± ^ 1 C m (yi,y 2 , ■■■,y n ,xu) ' ' < ^ 7 ' 63 ' ) 
Proceeding on the same lines as for i^,-, . we get K,a 

^i m ax = (7.64) 

And thus 

tn \ /i n 2 n n n \ ^ 

E / W >< u E £ + 4 E E 2 ^ + 2 ^o E 7 p- 65 ) 

t=i / V i=i ^ i=i j=i ^ i=i Vi J 

where i ^ j and j < i. 

The equations giving the maximum projections are: 



/ \ 1 Ce{yi,y2, ■ ■ -,y n ,xi 2 ) . 

Pn{yi, 2/2, • ■ ■ , y n , xi 2 ) = - < 5 I x 

z \ V 2^i=i 
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7.6 Summary and Conclusions 

We have generalized the Caledonian Symmetric Five Body Problem (CS5BP) 
to the Caledonian Symmetric N Body Problem (CSNBP) and derived a general 
form of the stability criterion involving the Szebehely constant C . Research 
in chapters 6 and 7 show that the general form of the stability criterion can be 
used for four and five body symmetrical systems. Szell (2003) has shown that 
a guarantee of complete hierarchical stability for all time is not possible for the 
six or more body problems because, unlike the phase space of the CS5BP and 
CSFBP, the phase space of the CSNBP is a connected manifold and transition 
from one hierarchy to another is always possible. However it is hoped that 
the general stability criterion can still be useful in determining restrictions 
placed on the ability of the system to change from one specific hierarchy to 
another in cases of more than five body configurations. It may also be possible 
to find a hierarchical stability criterion that gives partial stability for some 
hierarchies in a similar manner as the criterion R 3 < C < i? 4 found for the 
four and five body problem. Also, although there may be no guarantee of 
hierarchical stability for all time for Co > C cr u, it may be that numerically 
and effectively there is hierarchical stability, since the likelihood of changing 
from one hierarchy to another becomes very small as the connections between 
the hierarchies in the phase space narrow. In other words, a measure of the 
size of C relative to C crit may well give an indication of when a system is very 
unlikely to change its hierarchy, even though theoretically it is still possible. 
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Chapter 8 

Summary and Conclusions 



We have investigated the Caledonian Symmetric N Body Problem with par- 
ticular focus on the Caledonian Symmetric Four Body Problem (CSFBP) and 
Caledonian Symmetric Five Body Problem (CS5BP). Both the CSFBP and the 
CS5BP are symmetrically restricted problems with all possible symmetries uti- 
lized. We have also studied, numerically, a slightly perturbed non-symmetric 
extension of the CSFBP by using the general four body equations. 

In this chapter we provide a brief summary of the results obtained in this 
thesis. For this purpose we divide this chapter into four sections. In Section 8.1 
we summarize the results obtained in Chapter 3 which includes equilibrium so- 
lutions of four body problem and its linear stability analysis. In Section 8.2 we 
summarize the results of Chapter 4 dealing with the slightly perturbed Cale- 
donian Symmetric Four Body Problem (CSFBP). In Section 8.3 we summarize 
the results of both analytical and numerical work undertaken to understand 
the global characteristics of the Caledonian Symmetric Five Body Problem 
(CS5BP). In Section 8.4 we summarize the results of Chapter 7 dealing with 
the Caledonian Symmetric N Body Problem (CSNBP). 
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8.1 Equilibrium configurations of four body 
problem and its linear stability analysis 

The main goal of the first part of the thesis was to review the literature on 
analytical solutions to four and five body problems (chapter 2), to find equi- 
librium solutions for the four body problem (Chapter 3) and to discuss their 
linear stability (Chapter 4). 

In Section 3.1 of chapter 3 we gave a detailed review of the equilibrium 
configurations of equal mass cases of the four body problem (Steves and Roy, 
2001). In Section 3.2 we discussed the non equal mass cases of the four body 
problem. The non equal mass cases have two pairs of equal masses of masses 
m and M with m being the smaller one. We defined the ratio between m 
and Mas/i = m/M < 1. We allowed \i to reduce from 1 to to obtain 
the Lagrange five equilibrium points L 1 ,L 2 ,L 3 , L 4 and L 5 of the Copenhagen 
problem (Steves and Roy, 2001). We completed the analysis of Steves and Roy 
(2001) to include two more examples of equilibrium solutions of four body 
problem: 

1. The Triangular equilibrium configuration of four body problem with the 
two bigger masses making the base of the triangle. 

2. The Triangular equilibrium configuration of four body problem with the 
two smaller masses making the base of the triangle. 

We allowed the two smaller masses to reduce from 1 to 0. Figures (3.9) and 
(3.11) showed the evolution of all the four masses when \x goes to zero from 
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1 in the above two cases of the triangular equilibrium configurations. In the 
Triangular Case-I for /x — 1 we get the two well known equal mass solutions 
i.e. the equilateral triangle solution and the isosceles triangle solution. For the 
triangular case-I as /x is reduced from 1 to there always exists two continuous 
family of solutions: 1) Solution 1 starting at /x = 1 with the isosceles triangle 
solution and ending as /x — > with P2 — > L4; and 2) solution 2 starting at /x — 1 
with the equilateral triangle solution and ending as \i — > with P4 — > Li. In 
the Triangular Case-II for /x = 1 the equilateral triangle solution is obtained. 
In this case there is a continuous family of solutions for /x between 1 and 0.9972. 

8.2 Stability analysis of the CSFBP 

In chapter 4 we analyzed the evolution of the nearly symmetric Caledonian 
Symmetric Four Body systems over a range of initial conditions. We identified 
the stable and unstable regions of orbits by integrating over 1 million time 
steps. The possible endpoints of the numerical integration were: 1. close 
encounter between any two of the bodies occurs, 2. the dynamical Symmetry 
is broken or 3. 1 million time steps of integration. Two types of systems were 
analyzed, the symmetric and nearly symmetric CSFBP by using a general 
four body integrator. The general four body integrator was developed using 
the Microsoft Visual C++ software. The results obtained were processed using 
Matlab 6.5 to produce the graphs given in Figures (14.31) to (I4.18P 

We analyzed the phase space of the symmetric and nearly symmetric CSFBP's 
in detail for the /x = 1 and /x = 0.1 cases. Graphs of the phase space of the 
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CSFBP's with different Co values were determined. Each grid point r 1? r 2 on 
the graphs denoting initial conditions of an orbit, Figures H4.3[) to (14.181) . was 
color coded according to the different final outcomes of the evolution of the 
orbits. 

We found that the size of the stable regions in the phase space of the 
CSFBP is dependent on the value of the Szebehely constant. The larger the 
value of the Szebehely constant the more stable overall is the CSFBP system 
across the phase space. For the equal mass case of the CSFBP the phase space 
was very chaotic as most of the orbits were collision orbits with the single 
binary regions being the most chaotic. In the single binary cases, almost all 
of the orbits were collision orbits. In the /i = 0.1 case, the SB2 region is 
the most chaotic and was comparatively very very small. The SB1 region 
is the most stable as overall there were few collisions. When comparing the 
two cases, /i = 1 and /x = 0.1 with each other, the \x = 0.1 case is the most 
stable. In the \i = 1 large number of orbits fail the symmetry breaking 

criterion. We know that those orbits are not necessarily unstable. Therefore it 
will be interesting to investigate, in future, the long term stability of the equal 
mass case of the symmetric and nearly symmetric CSFBP past the symmetry 
breaking point. 
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8.3 Topological stability of the CS5BP and nu- 
merical verification of the analytical re- 
sults 

The main goal of this part of the thesis (Chapter 5, 6) was to derive an analyt- 
ical stability criterion to discuss the global hierarchical stability of the CS5BP 
and then to confirm the prediction of these theoretical results using numerical 
integrations. In order to do the first part we introduced a stationary mass 
to the centre of mass of the CSFBP. The resulting system was named The 
Caledonian Symmetric Five Body problem (CS5BP). 

In Section 1 of Chapter 5 we defined the CS5BP in such a way that the 
CSFBP became a special case of the CS5BP. We derived the Sundman's in- 
equality in its simplest possible form which is the key to the derivation of the 
stability criterion. In Section 3 to 5 we derived the analytical stability criterion. 
We showed that the hierarchical stability of the CS5BP for a range of different 
mass ratios solely depends on the Szebehely constant, Co which is a function 
of the total energy and angular momentum. We also showed that Co > 0.0659 
will guarantee that the CS5BP for all values of the mass ratios = 0, 1,2 
will be hierarchically stable. It is also shown that when /x > A*i or /i 2 , the 
14 hierarchy state is hierarchically stable for R 3 < Cq < while in all other 
cases the 24 hierarchy state is hierarchically stable at R 3 < C < R 4 . In all 
cases hierarchical stability for all states occurred when Co > -R4. Of course 
when /ii — fj, 2 , hierarchical stability for all time occurs when Co > -R3 = Ra 

In Chapter 6 we confirmed the theoretical results of Chapter 5 with the 
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help of numerical integrations, i.e. we showed that the Szebehely constant 
can be used to predict hierarchical stability. We also confirmed that when the 
value of Co was increased, the system became more hierarchically stable. In 
order to do this an appropriate form of equations of motion were derived and a 
numerical integrator was developed using Fortran and Microsoft Visual C++ 
software. Initial conditions were set in such a way that the whole phase space 
was covered. We used the Matlab 6.5 Software package to process the results. 
The results of thousands of integrations were analyzed in Section 5 and 6 and 
the results were collected in different kinds of tables. 

The table of hierarchy changes contained the number of hierarchy changes 
throughout the whole range of different integrations. It can easily be seen 
from these tables that as we increase the value of the Szebehely constant, the 
number of hierarchy changes drop until it becomes zero for C > C crit = i? 4 . 
Therefore it confirms the analytical predictions of Chapter 6. For values of C 
, i? 3 < C < i? 4 , the analytically predicted hierarchically stable states were 
verified as well. These tables also showed that the introduction of a stationary 
mass into the centre of mass of the CSFBP played a stabilizing role. The 
number of hierarchy changes decreased as we increased the value of the central 
mass. Also when one of the fii or ji 2 was very small we had fewer number of 
hierarchy changes. 

Our analysis in Chapter 6, also extended and completed the analysis of 
Szell (2003) for the CSFBP, since he checked the hierarchical stability of the 
CSFBP for only half of the phase space. We have discussed the CSFBP as a 
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special case of the CS5BP, where /x = 0, and covered the whole phase space. 

8.4 The Caledonian Symmetric N Body Prob- 
lem (CSNBP) 

In Chapter 7 we generalized the Caledonian Symmetric Five Body Problem 
(CS5BP) to the Caledonian Symmetric N Body Problem (CSNBP) which is the 
most general form of the Caledonian problem. Szell (2003) showed that there 
cannot be any hierarchical stability criterion for the six or more body problem 
which is valid for all time. Our generalization showed that the Szebehely 
constant still exists and still defines the topological behavior of the phase 
space. We have shown in Chapter 5 and 6 that there is a relation between the 
Szebehely constant and the hierarchical stability. Therefore it is very likely 
that this constant will also play some role in determining the hierarchical 
stability of the N Body Problem where n > 6. It has been proven by Szell 
(2003) that it cannot give complete hierarchical stability criterion for all time, 
but perhaps it can give a criterion that stipulates the hierarchical stability of 
certain states of hierarchy as R3 did in the CSFBP and CS5BP. Or perhaps 
even though there is no guarantee of hierarchical stability for all time using 
Co as a stability criterion, it may be that the criterion can provide a measure 
that in effect provide stability, i.e. analytically it is possible for the system to 
change from one hierarchy to another but effectively it does not because the 
chances of this occurring are so slim. We will address these question in the 
future. 
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Chapter 9 
Future Work 



We investigated Caledonian Symmetric N Body Problems. The main focus was 
on the Caledonian Symmetric Five Body Problem (CS5BP) with the Caledo- 
nian Symmetric Four Body Problem (CSFBP) studied as a special case. Re- 
search on the CS5BP and hence the CSFBP is far from complete. In every 
chapter we answered some very interesting questions but with every answer 
more questions arose which we would like to answer in the future. We will list 
some of these questions in this chapter with some explanations on how these 
questions might be addressed. 

Section 9.1, 9.2 and 9.3 discusses the questions still unanswered, with some 
proposed methods of tackling the solutions for the following areas respectively, 
CSFBP, CS5BP and CSNBP. 

9.1 Future lines of research on the CSFBP 

In Chapter 4 we investigated the phase space of the CSFBP for stable and 
unstable regions using the general four body integrator. We investigated the 
four body problem with two types of initial conditions i.e. symmetric initial 
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conditions and nearly symmetric initial condition where symmetry was broken 
slightly in the ^-component of Pi. 

In the equal mass case there are a large number of orbits which fails the 
symmetry criterion. We know that because an orbit fails to maintain its sym- 
metry does not necessarily mean that it is unstable. Therefore it will be inter- 
esting to repeat the analysis of Chapter 4 for those regions of the phase space 
of the CSFBP, in the equal mass case , that failed the symmetry, without the 
symmetry restrictions. This research can be done with our existing numerical 
tools except for a slight modification in the integrator. We could use the same 
equations of motion and initial conditions given in Chapter 4. 

The long-term integration of the general four body equations is a very 
important exercise but at the same time very time consuming too. For example 
only long-term integrations can tell you which of the perturbed system will 
remain close to the symmetric CSFBP systems. There is no alternate method 
for this particular exercise except long-term integration. 

There are, however some fast chaos detection methods which can be used to 
study the chaotic behavior of the CSFBP. Szell, erdi, Sandor and Steves (2003) 
have investigated the chaotic behavior of the CSFBP using the fast chaos de- 
tection methods of Relative Lyapunov Indicators (RLI) and Smaller Alignment 
Indices (SALI). By using the general four body equations, the CSFBP system 
becomes closer to the general four body system. It would be interesting to 
investigate the phase space of the CSFBP by using the fast chaos detection 
methods like RLI and SALI with the general four body equations. 
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This exercise will provide us with the unique opportunity of comparing 
these results with our results obtained in Chapter 4 and Szell, erdi, Sandor and 
Steves (2003) results. We expect to get similar results to Szell, erdi, Sandor and 
Steves (2003) about the connection between the chaotic and regular behavior 
of the phase space of the CSFBP with the Szebehely constant i.e. the phase 
space of the CSFBP is more regular for Co near C cr u and particularly for 
Co > C cr it. We also think that by using the general four body integrator 
the phase space of the CSFBP will be more chaotic which is indicated by the 
results in Chapter 4. To be able to do this exercise we would have to modify 
the integrator or perhaps write a new integrator. 

Another aspect of the CSFBP which is completely unexplored is an inves- 
tigation of the non-coplanar case. The degrees of freedom will increase which 
will make it comparatively more challenging than the coplanar case. The an- 
alytical results obtained by Steves and Roy (2001) will not change, but we 
will have to modify the equations of motion and hence the integrator to carry 
out this research. The three dimensional CSFBP model will be a very good 
approximation of the general four body model near a symmetrical system. 
The results obtained in this case could be used to model three dimensional 
symmetrical quadruple stellar and planetary systems. 

We believe that apart from being an approximation to the general four 
body problem, the CSFBP has theoretical as well as practical applications. 
The phase space of the CSFBP is a subspace of the phase space of the general 
four body problem. Therefore some typical behavior of the CSFBP may occur 
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in the general four body problem. We would like to examine how the dis- 
connectedness of the restricted phase space influences the global phase space 
nearby. 

9.2 Future lines of research on the CS5BP 

In Chapter 5 we developed a symmetrically restricted five body problem called 
the Caledonian Symmetric Five Body problem (CS5BP), for which we derived 
an analytical stability criterion valid for all time. Following in the footsteps of 
Steves and Roy (1998, 2000, 2001) we have shown that the hierarchical stability 
of the CS5BP depends solely on a parameter called the Szebehely constant, C . 
We verified this stability criterion numerically in Chapter 7. The CS5BP is 
modelled by introducing a stationary mass to the centre of mass of the CSFBP 
and we have seen that the bigger the central mass, the more stable the CS5BP 
system becomes. This generates some very interesting questions. 

The characteristics of the CS5BP studied in Chapter 5 and 6 are its global 
characteristics. In future we would like to look at the finer details of the 
CS5BP and study it locally. We will repeat the exercise of Chapter 4 for the 
CS5BP i.e. we will perform long term integrations for different mass ratios 
of the CS5BP with the general five body integrator for both symmetric and 
nearly symmetric cases. This will help us to discover the effect on the stability 
of the general system near the symmetric system, of adding a central body to 
the CSFBP. We will also be able to find out if the system remains close to 
symmetry for a longer time than it was found in the CSFBP. 
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To be able to complete this research we would need to extend the general 
four body integrator to include five bodies. We would also need to find suitable 
initial conditions and derive the equations of motion for the general CS5BP. 
We should also be prepared to dedicate a lot of CPU time. 

Another method of studying the CS5BP locally is to use the fast chaos 
detection methods like RLI and SALI. Proceeding along the same lines as we 
suggested for the CSFBP, we plan to study the chaotic and regular behavior 
of the CS5BP. In this research we will have the flexibility to vary the value 
of the central mass. Therefore we can approximate this CS5BP system by 
CSFBP by approximating the central mass to zero. This exercise will provide 
us with the unique opportunity of comparing these results with the results 
gained by investigating the long-term evolution of CS5BP orbits that are both 
symmetric and nearly symmetric. We expect to get results similar to the long 
term integrations. 

Like the CSFBP the non-coplanar case of the CS5BP is completely unex- 
plored. Recall that in this case the degrees of freedom increases in the non- 
coplanar case and therefore will be more challenging to explore. The method 
of research will be the same as was for the CSFBP. In fact we can analyze the 
CS5BP first and then let the mass of the central body be equal to zero in order 
to analyze the CSFBP. The three dimensional CS5BP model could be a good 
approximation of the general five body model. These results can be used to 
model planetary systems. 

After completing the proposed research on the CS5BP in this section we 
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should be able to answer the following very important question. 

Does the central mass of the CS5BP forces the system to remain close to 
symmetry for longer than the CSFBP and what role does the central mass play 
in the overall stability/instability of the CS5BP in the more general problem. 

Our current analysis of the CS5BP suggests that if the mass of the central 
body is bigger than the other bodies, then the system will be more stable. 

9.3 Future lines of research on the CSNBP 

In chapter 7 we have generalized the CS5BP to any number of bodies to give 
the CSNBP. In the CSNBP there exist critical values of the Szebehely constant 
but we know that they do not predict hierarchical stability for all time because 
unlike the phase space of the CS5BP the phase space of the CSNBP is a 
connected manifold and transition from one hierarchy state to another is always 
possible. Although complete hierarchical stability is not possible, we think that 
the Szebehely constant could still have some role to play in the hierarchical 
stability of more than five body systems. For example for a very large value 
of C , greater than C crit , the six or seven body system should be more stable 
than it will be for a smaller C . We plan to investigate this numerically and 
confirm or contradict our expectations. 

A very interesting application of the CSNBP can be in Galactic Dynamics 
or stellar cluster dynamics. For example by taking N to be a very large number 
or in other words by taking N — > oo it may be possible to model globular 
clusters or galaxies and their stability at least in an empirical sense. 
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Appendix A 

Results of ? on the numerical investigation of 
hierarchical stability of the CSFBP 

In this appendix results on the numerical investigations of the hierarchical 
stability of the CSFBP are given (?), which include tables on total number of 
hierarchy changes and percentage of hierarchy changes for different values of 
fi. 
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Table 1: Hierarchical changes for // = 1 
Total number of hierarchy changes 





Co = 10 


Co = 29 


Co = 30 


Co = 40 


Co = 46 


Co = 47 


Co = 60 


12 13 





1 





3 











12 14 


230 


136 


120 


79 


7 








12 23 


143 


111 


74 


87 


10 








13 12 





2 


2 


4 











13 14 


133 


79 


88 


68 











13 23 


64 


97 


75 


69 


9 








14 12 


137 


106 


97 


56 


2 








14 13 


95 


112 


108 


65 


4 








14 23 


1119 


185 


158 


42 











23 12 


80 


164 


151 


170 


56 








23 13 


217 


198 


191 


184 


57 








23 14 


1316 


267 


236 


42 


1 












Percentage of hierarchical changes 








Co = 10 


C = 29 


C = 30 


C = 40 


C = 46 


Co = 47 


C = 60 


12 13 


0.0 


0.1 


0.0 


0.3 


0.0 


0.0 


0.0 


12 14 


6.5 


9.3 


9.2 


9.1 


4.8 


0.0 


0.0 


12 23 


4.0 


7.6 


5.7 


10.0 


6.8 


0.0 


0.0 


13 12 


0.0 


0.1 


0.2 


0.5 


0.0 


0.0 


0.0 


13 14 


3.8 


5.4 


6.8 


7.8 


0.0 


0.0 


0.0 


13 23 


1.8 


6.7 


5.8 


7.9 


6.2 


0.0 


0.0 


14 12 


3.9 


7.3 


7.5 


6.4 


1.4 


0.0 


0.0 


14 13 


2.7 


7.7 


8.3 


7.5 


2.7 


0.0 


0.0 


14 23 


31.7 


12.7 


12.2 


4.8 


0.0 


0.0 


0.0 


23 12 


2.3 


11.2 


11.6 


19.6 


38.4 


0.0 


0.0 


23 13 


6.1 


13.6 


14.7 


21.2 


39.0 


0.0 


0.0 


23 14 


37.2 


18.3 


18.2 


4.8 


0.7 


0.0 


0.0 
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Table 2: Hierarchical changes for /i = 0.1 
Total number of hierarchical changes 





C = 0.4 


C = 0.6 


C = 0.65 


Co = 0.7 


C = 0.8 


C = 0.9 


12 13 





1 


3 


2 








1 O 1 A 

12 14 


loo 


a n 


41 


51 


Q 
O 


U 


1 1 OQ 
1Z Zo 


P.A 
04 


/in 
4a 


32 


38 


U 


n 
U 


lO 1Z 


z 


n 
U 








n 
U 


n 
U 


lO 14 


08 
Zo 


QD 

ou 


37 


50 


i 


n 
U 


1 Q OQ 
lo Zo 


QQ 
OO 


QQ 
OO 


38 


18 


n 
u 


n 
u 


1/1 10 
14 1Z 


QS 
OO 


07 
Z I 


31 


37 


i 


n 
u 


1/1 1 Q 
14 IO 


Q1 
Ol 


07 
Z I 


32 


30 


i 


n 
u 


1 /I OQ 
14 ZO 


A 81 
4ol 


oon 
zzu 


147 


90 


Q 
O 


n 
u 


23 12 


69 


65 


57 


62 








OQ 1 Q 
ZO 10 


A 7 
4 I 


fiO 

oz 


64 


61 


n 
u 


n 
u 


OQ 1 /I 
Zo 14 


c;/i i 

041 


OA A 
Z44 


177 


84 


i 


n 
u 








5e of hierarchical chang 








Co = 0.4 


C = 0.6 


C = 0.65 


Co = 0.7 


C = 0.8 


C = 0.9 


12 13 


0.0 


0.1 


0.5 


0.4 


0.0 


0.0 


12 14 


11.1 


6.1 


6.2 


9.8 


30.0 


0.0 


12 23 


4.3 


6.1 


4.9 


7.3 


0.0 


0.0 


13 12 


0.1 


0.0 


0.0 


0.0 


0.0 


0.0 


13 14 


1.9 


3.7 


5.6 


9.6 


10.0 


0.0 


13 23 


2.2 


4.1 


5.8 


3.4 


0.0 


0.0 


14 12 


2.5 


3.3 


4.7 


7.1 


10.0 


0.0 


14 13 


2.1 


3.3 


4.9 


5.7 


10.0 


0.0 


14 23 


32.1 


27.3 


22.3 


17.2 


30.0 


0.0 


23 12 


4.6 


8.1 


8.6 


11.9 


0.0 


0.0 


23 13 


3.1 


7.7 


9.7 


11.7 


0.0 


0.0 


23 14 


36.1 


30.2 


26.9 


16.1 


10.0 


0.0 
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Table 3: Hierarchical changes for ji = 0.01 



Total number of hierarchical changes 





C = 0.2 


Co = 0.27 


Co = 0.283 


Co = 0.29 


Co = 0.3 


12 13 











1 





12 14 


17 


16 


7 


6 





12 23 


14 


7 











13 12 





1 





1 





13 14 


3 


1 


6 


5 





13 23 


3 


3 


4 








14 12 


3 


4 


1 


6 





14 13 


1 





3 


7 





14 23 


105 


38 


13 


2 





23 12 


2 


4 


4 








23 13 


5 


6 


8 








23 14 


113 


41 


7 


2 






Percentage of hierarchical changes 





Co = 0.2 


C = 0.27 


Co = 0.283 


C = 0.29 


Co = 0.3 


12 13 


0.0 


0.0 


0.0 


3.3 


0.0 


12 14 


6.4 


13.2 


13.2 


20.0 


0.0 


12 23 


5.3 


5.8 


0.0 


0.0 


0.0 


13 12 


0.0 


0.8 


0.0 


3.3 


0.0 


13 14 


1.1 


0.8 


11.3 


16.7 


0.0 


13 23 


1.1 


2.5 


7.5 


0.0 


0.0 


14 12 


1.1 


3.3 


1.9 


20.0 


0.0 


14 13 


0.4 


0.0 


5.7 


23.3 


0.0 


14 23 


39.5 


31.4 


24.5 


6.7 


0.0 


23 12 


0.8 


3.3 


7.5 


0.0 


0.0 


23 13 


1.9 


5.0 


15.1 


0.0 


0.0 


23 14 


42.5 


33.9 


13.2 


6.7 


0.0 
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Table 4: Hierarchical changes for /i = 0.001 
Total number of hierarchical changes 





Co = 0.25 


C = 0.2529 




Co 


= 0.2533 


Co 


= 0.2538 


Co 


= 0.2545 


12 13 























12 14 





1 






5 












12 23 


2 


2 



















13 12 























13 14 












2 












13 23 


1 


2 



















14 12 












1 












14 13 





1 






3 












14 23 


56 


4 






2 












23 12 





2 






1 












23 13 


1 


1 






1 












23 14 


75 


9 






1 
















Percentage 


of hierarchical changes 








C = 0.25 


Co = 0.2529 




Co 


= 0.2533 


Co 


= 0.2538 


Co 


= 0.2545 


12 13 


0.0 


0.0 




0.0 




0.0 




0.0 


12 14 


0.0 


4.5 






31.3 




0.0 




0.0 


12 23 


1.5 


9.1 






0.0 




0.0 




0.0 


13 12 


0.0 


0.0 




0.0 




0.0 




0.0 


13 14 


0.0 


0.0 






12.5 




0.0 




0.0 


13 23 


0.7 


9.1 






0.0 




0.0 




0.0 


14 12 


0.0 


0.0 






6.3 




0.0 




0.0 


14 13 


0.0 


4.5 






18.8 




0.0 




0.0 


14 23 


41.5 


18.2 






12.5 




0.0 




0.0 


23 12 


0.0 


9.1 






6.3 




0.0 




0.0 


23 13 


0.7 


4.5 






6.3 




0.0 




0.0 


23 14 


55.6 


40.9 






6.3 




0.0 




0.0 
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